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Abstract 

We consider general physical characteristics of doubly heavy baryons: the spectroscopy 
in the framework of potential approach and sum rules of QCD, mechanisms of production 
in various interactions on the basis of fragmentation model with account of preasymptotic 
corrections caused by higher twists over the transverse momentum of baryon, inclusive decays 
and lifetimes in the operator product expansion over the inverse powers of heavy quark masses 
as well as the exclusive decays in the sum rules of QCD. We generalize the methods developed 
in the effective theory of heavy quarks towards the description of systems with two heavy 
quarks and a single light quark. The calculations are presented for the masses, decay widths 
and yields of baryons with two heavy quarks in the running and planned experimental facilities. 
We discuss the prospects of search for the baryons and possibilities of experimental observation. 
The most bright physical effects concerning these baryons as well as their position in the system 
for the theoretical description of heavy quark dynamics are considered. 
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Introduction. 



After the precision investigations of neutral intermediate Z boson at LEP (CERN) and observation 
of t quark at FNAL were done, the study of electroweak interactions in the sector of heavy quarks is a 
task among the most actual problems in the physics of fundamental particles, since in the framework 
of standard model |I| such the searches can result in a complete picture of effects responsible for the 
unreversability of time at energies below the scale of electroweak symmetry breaking. A comparative 
analysis of processes with violation of combined CP parity inverting both the charges (C) and space 
orientation (P) in decays of hadrons containing the heavy quarks becomes accessible for quite precise 
measurements due to putting into action specialized facilities Belle (KEK) and BaBar (SLAG) along 
with the upgraded detectors of CDF and DO at FNAL. The results of such experiments probably 
allow us to add an essential missing chain-cell in the standard model of interactions, i.e. a complete 
description of charged currents for three quark generations [0. This study is the most important 
problem along with both the observation of scalar Higgs particles providing the mechanism of 
electroweak symmetry breaking, and the investigation of neutrino currents. 

The pricision study of electroweak properties of heavy quarks make rise a deep theoretical prob- 
lem, that is the description of dynamics for the strong interactions of quarks, since the strong forces 
cause the formation of bound states, hadrons: mesons and baryons, and the observable quantities 
like characteristics of rare decays with the CP violation are assigned to the bound states, so that 
we need definite and reliable representations about the relation of measured properties with the 
parameters of heavy quark interactions themselves. In this way, fine effects of electroweak physics 
can be extracted under a high precision description for the dominant contributions of quantum chro- 
modynamics (QCD). In this respect, we deal with a general problem on quantitative understanding 
the quark confinement in QCD, that can be efficiently investigated not only in the spectroscopy 
and processes with the production and decays of exotic hybrid and glueball states P], but also in 
the study of hadrons with heavy quarks. In practice, the measured quantities of asymmetries in the 
decays of heavy hadrons, for instance, are expressed in termes of functions parametrically depend- 
ing on both the primary characteristics assigned to charged and flavor changing neutral currents of 
quarks, and the hadronic matrix elements of quark operators, which usually cannot be straightfor- 
wardly determined from a wide set of various experimental data. So, we need a sound theoretical 
anlaysis of such the matrix elements in QCD. 

A more complete list of heavy quark bound states under study provides a more wide region 
for the variation of conditions, wherein the forces of QCD act on the heavy quarks, so that the 
theoretical methods for the description of hadrons containing the heavy quarks should be more 
accomplished in order to get a consistent understanding of various quark systems. In this way, we 
see a new field of activity is baryons, containing two heavy quarks. So, the theoretical predictions for 
the properties of such the baryons are of interest and actual. The doubly heavy baryons naturally 
continue the list of long-lived heavy hadrons with both a single heavy quark (D, B mesons and A^, 
Ec, Sc, fie, Afe baryons) and two heavy quarks (the meson). With respect to the character of 
strong interactions, these baryons could have common features with the heavy quarkonia cc and hh. 
In practice, we could expect experimental observation of doubly heavy baryons in the searches at the 
modern hadron colliders with high lumonosities (Tevatron, LHC), since the yield of such baryons 
should be of the same order of magnitude as for the doubly heavy meson containing the quarks of 
different flavors, the B^. mesonQ, while the methods for the registration of rare decays with heavy 
particles recently get a high efficiency due to a development in the technique of vertex detectors as 

^See the review on the physivs of i?c in Q. 
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was successfully exhibited in the first experimental observation of Be by the CDF collaboration 

Constructing the theoretical methods for the description of QCD dynamics with the heavy quarks 
is based on a clear physical definition: the quark Q is ascribed heavy, if its mass ttlq is much greater 
than the scale of heavy quark confinement in the bound state, so niQ 3> Aqcd- Thus, considering 
the problem of strong interactions with the heavy quarks, i.e. calculating hadronic matrix elements 
of quark operators, we could involve a small parameter AqcD/mQ 1, which could serve for the 
development of formal approximate methods. 

So, first of all, in hard processes with virtualities of the order of heavy quark masses, say, in 
the production of heavy quarks, the QCD coupling constant is small, ~ inmQ/AQca ' 
usual technique of perturbative theory in powers of is justified. Another productive method is the 
operator product expansion (OPE) in inverse powers of heavy quark mass. In such the approach a 
calculation of hadronic matrix element for a quark operator leads to summing up the matrix elements 
of operators, whose properties assume a hierarchy in terms of small parameter Kqcd/^q ^ 1, i-e. 
a suppression of some contributions by powers of Kqcd/i^q, since the interaction in the bound state 
containing the heavy quark is characterized by the energies close to Aqcd determining the inverse 
size of the hadron. This expansion gives a complete description of hadronic systems containing 
a single heavy quark. If there are two heavy quarks in the hadron, then along with the scale 
of nonperturbative interactions another energetic scale is a momentum transfer in the coulomb 
interaction with a virtuality of /i ~ a^mg, so that a relative velocity for two heavy quarks moving 
inside the hadron, v, is determined by a comparatively small value of coupling constant in QCD, and 
V ^ as, where the coupling constant is taken at the scale of virtualities prescribed to the coulomb 
interactions. Thus, in the heavy quarkonium QQ', say, there is an additional small parameter, 
that is the relative velocity of nonrelativistic quarks v, which can be used in the OPE in order to 
calculate the hadronic matrix elements. Following such the approach, we underline three methods 
for the evaluation of quantities characteristic for the bound states with heavy quarks: 

- the operator product expansion in the inverse powers of heavy quark mass in QCD for calcula- 

tions of inclusive width and lifetimes of heavy hadrons, where terms correcting the leading 
approximation are given by some external parameters |Q; 

- sum rules of QCD and nonrelativistic QCD (NRQCD) for both two-point correlators of quark 

currents in spectroscopic calculations and three-point correlators in estimations of exclusive 
decay modes |]^; 

- potential models used for the evaluation of exclusive characteristics of hadrons containing the 

heavy quarks [§]. 

Let us stress, first, that the sum rules are based on the operator product expansion, too. However, 
the only external parameters of sum rules are fundamental quantities such as the masses of heavy 
quarks, the normalization of QCD coupling constant and the quark-gluon condensates in contrast 
to inclusive estimates in the general OPE, wherein we should use, for instance, the value of heavy 
quark binding energy in the hadron, the average square of heavy quark momentum etc. Second, in 
the expansion over the inverse powers of heavy quark mass as well as over the relative velocity of 
heavy quarks inside the hadron an important role is played by both the perturbation theory and 
the renormalization group relations, which are necessary for the calculation of Wilson coefficients 
in the operator product expansion, since these coefficients enter the factors in front of operators or 
matrix elements under consideration. 
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It is worth to emphasize that in the method of OPE for the heavy quarks we could consider an 
actual operator, say, a current of weak decay or a product of currents as it does in sum rules of 
QCD, with the consequent expansion, while the approach of effective field theory could be useful, 
too, wherein the starting point of construction is a formal expansion of QCD lagrangian itself for the 
heavy quarks. Then we can isolate a leading term in the effective lagrangian and treat the rest of 
terms as perturbations in the expansion. In this way the definition of leading term is determined by 
the character of the problem, i.e. by the actual convergency in the estimation of physical quantities 
calculated in the effective lagrangian. So, for the hadrons with a single heavy quark the effective 
theory of heavy quarks (HQET) |Q was developed. In this theory we can neglect the binding 
energy of heavy quark inside the hadron to the leading order, therefore, the kinetic energy of heavy 
quark is the perturbation. It is important to note that, first, in the leading order of HQET the 
effective lagragian possesses the following symmetries: a) the heavy quarks with identical velocities 
equal to the velocity of hadron, where the quarks are bound, can be permuted, that implies the 
heavy flavor symmetry, b) the heavy quark spin is decoupled from the interaction with low-virtual 
gluons, since the current is given by the quark velocity Vu, that implies the spin symmetry. Second, 
the leading term of effective lagrangian provides the renormalization group behaviour different 
from the full QCD. Particularly, the currents being conserved in full QCD and, hence, having 
zero anomalous dimensions become divergent after the transition to the fields of effective theory. 
The same note concerns to the correction terms of effective lagrangian: the corresponding Wilson 
coefficients in the effective theory have nonzero anomalous dimensions, too. Thus, we deal with the 
situation, wherein we need an infinite number of normalization conditions for the anomalous Wilson 
coefficients in the renormalization group in order to correctly define the theory. This problem has 
a clear physical reason, since in the effective theory constructed for the fields with low virtualities 
one should introduce a cut off in the ultraviolet region at the scale of the order of heavy quark mass 
because at high virtualities the assumptions made in the derivation of the theory are not correct. 
Since the effective theory was derived from the full QCD, we have to constructively put the effective 
lagrangian equal to the lagrangian of full QCD at a scale fihard about the heavy quark mass mg. 
This fact implies that in the given order over the coupling constant of QCD we have to calculate 
the effective action of QCD with account of corresponding loop corrections and to transform it by 
the expansion in inverse powers of heavy quark mass, so that the expanded lagrangian should be 
equalized to the effective lagrangian calculated in the effective theory in the same order of coupling 
constant with the anomalous Wilson coefficients at the scale fihard, that results in the matching 
conditions for the unknown constans of integration in the renormalization group equations for the 
Wilson coefficients in the effective theory. After the matching is done, we remove an ambiguity in 
the choice of finite contraterms in the renormalization of effective lagrangian, and it becomes definite 
at /i below the scale of matching with full QCD fitard- Thus, for the hadrons containing single heavy 
quark along with a light one the scheme of effective theory HQET is cosistently constructed. 

For the heavy quarkonia composed of heavy quark and heavy anti-quark another physical situ- 
ation takes place. Indeed, the coulomb interaction of nonrelativistic heavy quarks leads to that the 
kimetic energy is of the same order of magnitude as the potential energy, while one could naively 
expect that the kinetic term p^/2mQ should be suppressed by the heavy quark mass. However, 
since in the coulomb exchange we have p ~ Osmq, such the suppression of kinetic energy does not 
take place, and the wave functions of heavy quarkonia depend on the quark masses, i.e. they are 
flavor-dependent. In the formal approach of effective theory for the nonrelativistic quarks in the 
heavy quarkonium the leading term of lagrangian is defined with account of the kinetic energy, so 
that we deal with the nonrelativistic QCD (NRQCD) |ig. In NRQCD compared with HQET, the 
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Figure 1. The character of strong interactions in the doubly heavy baryon Ei^c- the compton lengths of quarks 
Aq = l/mq, the size of heavy diquark r^c ~ l/(f • mg) and the scale of nonperturbative confinement of light 
quark rqco = 1/Aqcd are arranged \b ^ \ ^H^^ ^ tqcd- 

spin symmetry of leading order effective lagragian survives, while there is no heavy flavor sym- 
metry because the contribution of kinetic energy explicitely depends on the quark masses. In the 
same manner as in HQET, the Wilson coefficients in the effective lagrangian of NQRCD have to be 
matched with the full QCD at a scale about the heavy quark mass, so that, in general, these coef- 
ficients have anomalous dimensions different from those of HQET, since the kinetic energy results 
in the ultraviolet behaviour of quark propagators, which is different from the behaviour in HQET. 

The potential approach is also based on the operator product expansion. So, the static potential 
means the expansion of effective action of QCD for two infinitely heavy sources j posed at a fixed 
distance r, so that the effective action has the form of r(j) = —V{r) ■ T, where T ^ oo is a time 
of the sources are switched on. For actual problems the long time interval means that virtualities 
of external gluon fields interacting with the heavy quarks are much less than the inverse distance, 
i.e. /i ~ |; -C i ~ rriQV. This constraint for the consistency of potential approach can be expressed 
in terms of effective theory, which is called the potential nonrelativistic QCD (pNRQCD) ||1 1|| . 
The construction of pNRQCD follows the matching of NQRCD action with the effective action of 
supersoft fields treated in the multipole expansion of QCD with nonrelativistic heavy quarks at 
/i ~ f ■ TTiQ. The leading order includes the kinetic energy as well as the Wilson coefficient meaning 
the static potential depending on the distance between the quarks, r. Thus, the potential approach 
also has the status of OPE, and the static approximation for the potential is determined by the 
covergency of this expansion in pNRQCD. 

In this respect, the baryons with two heavy quarks especially are of interest for the theoretical 
consideration, since for their description one should develop and use a combined approach involving 
the features of HQET, NRQCD and pNRQCD because in this systems the interaction of light quark 
with the heavy quarks is essential as well as the interaction between the heavy quarks (see Fig. |l]). 

In this review we solve the problem on the description of bound baryonic states QQ'q = '^qq' 
with two heavy quarks Q, Q' and a light quark q on the basis of factorizing the interactions with 
various virtualities determined by the following: 

- the confinement scale Aqqd for the nonperturbative interactions of both the heavy quarks with 
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the light one and the quarks with the quark-lguon condensates, 

- the size of heavy diquark tqq/ ~ 1/ {itlq^qi ■ v) , which is composed of two nonrelativistic heavy 

quarks moving with a smaU relative velocity -C 1, for the interactions between the heavy 
quarks, 

- the scale of hard gluon corrections at the energies about the heavy quark masses mg, 

so that we assume that the leading approximation is valid if there is the hierarchy of QCD interaction 
scales in Hqq/ 

Aqcd < ■ < niQ. (1) 

In this approach the doubly heavy diquark acts as a local heavy source of gluon field for the light 
quark. This source is charged as the anti-triplet under the color group of QCD, while the diquark 
itself is the system of two nonrelativistic quarks in the low-frequency field of light quark. Thus, for 
the motion of light quark and diquark we can use the effective theory of HQET, to the moment for 
the motion of heavy quarks inside the diquark we should modify NRQCD and pNRQCD in order 
to consider the nonrelativistic fields of heavy quarks in the anti-triplet state, but the singlet one. 

On the basis of quark-diquark representation for the interactions we consider various physical 
aspects of baryons with two heavy quarks. In Chapter 1 the mass spectrum of baryons Sqq/ 
is constructed in the potential approach. We calculate the characteristics of groud state and its 
excitations in the system of heavy diquark as well as in the system of light quark and diquark. We 
show that there is a family of Sqq/ levels with the quasi-stable states for the heavy diquark composed 
of identical heavy quarks. For these states the Pauli principle dictates quite definite values for the 
sum of quark spins, so, the total spin equals unit for the P-even wave functions in the configuration 
space, while it is equal to zero for the P-odd wave functions, since the anti-triplet color state of 
diquark is anti-symmetric under the permutation of color indices. Taking into account the small size 
of diquark, the nonrelativistic motion of heavy quarks and a small ratio of Aqcd/itt-q, the operators 
for the transitions of excited P-wave diquark into the ground S-wave level with the emission of n 
meson are suppressed because in this transition both the spin state of diquark and orbital state 
should change. We determine the region of consistency for the quark-diquark approximation in the 
calculations of Sqq/ mass spectra in the framework of potential approach. 

In Chapter 2 the two-point sum rules of NRQCD are considered for the baryonic currents with 
two heavy quarks. We discuss a criterium of stability for the results of such the sum rules in estimates 
of masses and coupling constants of Sqq/ baryons. We show that reliable results can be obtained 
after account of the quark and gluon condensates as well as their product and mixed consdensate, 
i.e. after the introduction of combined condensates of higher dimensions. We calculate the masses 
and coupling constants for the ground states of baryons Sqq/ and the doubly heavy baryons with 
strangeness floo', too. Reliable predictions for the mass splitting of — Mq are presented. 

We get anomalous dimensions for the baryonic currents in NQRCD up to two-loop approximation, 
that allows us to estimate the coupling constants of baryons with the baryonic currents not only in 
NRQCD but also in full QCD. 

Mechanisms for the production of Sqq/ baryons arc considered in Chapter 3. Following the 
factorization of interactions in the quark-diquark approach, in e+e~ annihilation at high energies 
the inclusive production of doubly heavy baryons can be described in the form of subsequent frag- 
mentation for both the heavy quark into the heavy diquark and the diquark into the baryon. In 
this way the virtualities in the first process are determined by the masses of heavy quarks, so that 
factorizing the soft movement of heavy quarks inside the diquark we can use the perurbative QCD 
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and get an analytic form of fragmentation functions for the states with various spins and orbital 
quantum numbers of diquark. Keeping in mind the difference between the color structures of di- 
quark and heavy quarkonium, these calculations repeat the consideration of fragmentation into the 
doubly heavy mesons up to a color factor. Supposing the representation of heavy diquark by the 
local field in its interactions with the light quark and approximating the field components of baryon 
by its dominant fast valence quarks, we develop a QCD-motivated perturbative model for the frag- 
mentation of diquark into the baryon and calculate an analytic form of fragmentation function for 
the vector and scalar diquarks entering the baryon Hqq/ with the spin i. The estimates for the pair 
production of baryons in the e~^e~ annihilation close to the threshold are presented, too. The anal- 
ysis turns to be more complicated in the consideration of mechanism for the production of baryons 
Sqq/ in hadron collisions due to subprocesses of quark-antiquark annihilation (low energies, fixed 
target experiments) and gluon-gluon fusion (high energies, hadron colliders). This complication is 
connected to a greater number of diagrams in the leading approximation, i.e. in the fourth order 
of QCD coupling constant. Making use of numerical method, we show that at high energies of 
partonic subprocess and at transverse momenta much greater than the baryon mass the complete 
set of diagrams in the given order of perturbation theory in QCD leads to the factorization for the 
production of heavy quark and its fragmentation into the doubly heavy diquark, expressed in the 
form of univerasl fragmetation function analytically derived in the perturbative QCD. This fact 
implies the consistency of approach used. The advantage of such considereation with the complete 
set of diagrams dictated by the gauge invariance is a possibility to calculate not only the leading 
term with respect to the transverse momentum p±, that gives the fragmentation dropping as ~ 

p± 

but also the correction terms, the higher twists over the transverse momentum. In this way we get 
a definite estimate for the transverse momentum determining the boundary between the regions 
of fragmentation and recombination (the higher twists). We point out that the statistics of events 
with the production of doubly heavy baryons Sqq/ is dominantly integrated out at small transverse 
momenta lying in the region of recombination. We represent estimates for the total and differential 
cross sections of Hqq/ baryon production in hadronic experiments at various energies and for the 
pair production of doubly heavy baryons in the quark-antiquark annihilation. 

The method of operator product expansion in the inverse powers of heavy quark masses is 
exploited in Chapter 4 for the analysis of lifetimes and inclusive width in decays of Sqq/ baryons. 
In this consideration of mechanisms for the decays of heavy quarks entering the doubly heavy 
baryons, the following three physical effects are essentially important: 

- valuable corrections to the spectator width in decays of heavy quarks appear because of movement 

of quarks in the heavy diquark, whose size is small and, hence, the relative momenta of quarks 
p ~ f ■ TTiQ are greater than the momentum of heavy quark in the hadron with the single heavy 
quark, when k ~ Aqcd, since we have v ■ uiq ^ Aqcd] 

- nonspectator contributions due to the Pauli interference between the products of heavy quark 

decay and the valence quarks in the initial state can give a fraction 30-50% of the total width, 
so that the feature of anti-symmetric color wave function for the baryons is a possibility of 
both a positive or negative overall sign for the term with the interference[| 

- along with the Pauli interference the weak scattering of quarks in the initial state appears in the 

OPE in the form of operator of higher dimension and it is enhanced by a two-particle phase 

^The overall sign is determined by the anti-symmetric permutation of fermions multiplied by the color factor. 
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space in comparison with other operators with the same dimension, since they have a three- 
particle phase space in the final state, so that the three-particle phase space is suppressed in 
units of heavy quark mass; the weak scattering gives about 30% of total width for the baryons 
with charmed quark. 

In this way, we show how the arrangement of lifetimes is produced for the baryons with two heavy 
quarks. Then we determine the parametric dependence of estimates for the total and inclusive 
widths on the physical quantities of hadron system. So, the small size of heavy diquark determines 
its wave function yielding the factor in the evaluation of nonspectator decays. The masses of heavy 
quarks in OPE are essentially correlated for the hadrons with different quark contents, so that 
current experimental data on the semileptonic, nonleptonic and total widths decrease uncertainties 
of estimates. Furthermore, the experimental data are able to essentially improve the qualitative and 
quantitative knowledges on the dynamics of heavy hadrons with the single or two heavy quarks. 

We analyze exclusive semileptonic decays and nonleptonic decays in the assumption of factor- 
ization in the framework of three-point sum rules of NRQCD, which allow us to derive relations 
for the formfactors of transitions given by hadronic matrix elements from the spin symmetry of 
effective lagrangian. We discuss uncertainties of calculations and compare the sum rule results with 
the predictions of potential models for the exclusive decays. 

In Conclusion we summarize our results on the physics of baryons containing two heavy quarks. 
The obtained predictions do not only point to a way for a goal-recognized search of such baryons, 
but also make a basis for a more accomplished and detailed theoretical analysis of physical effects 
in the hadronic systems with two heavy quarks, which are, no doubts, of interest, particularly, 
for the reliable predictions of total and exclusive widths. Finally, we consider possibilities for an 
experimental observation of doubly heavy baryons. 
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Chapter 1. Spectroscopy of doubly heavy baryons: 
the potential approach 



In this chapter we analyze basic spectroscopic characteristics for the famihes of doubly heavy baryons 
Sqq' = {QQ'q), where q = u, d and Qqq' = [QQ's). 

A general approach of potential models to calculate the masses of baryons containing two heavy 
quarks was considered in refs. [0. The physical motivation used the pair interactions between the 
quarks composing the baryon, that was explored in the three-body problem. Clear implications for 
the mass spectra of doubly heavy baryons were derived. So, for the given masses of heavy charmed 
and beauty quarks, the approximation of factorization for the motion of doubly heavy diquark and 
light quark is not accurate. It results in the ground state mass and excitation levels, essentially 
deviating from the estimates in the framework of appropriate three-body problem. For example, we 
can easily find that in the oscillator potential of pair interactions an evident introduction of Jacobi 

variables leads to the change of vibration energy uj "^/f^ comparison with naive expectations 
of diquark factorization. 



Figure 1.1. The representation of doubly heavy baryon QQq with the colored fields forming the strings between 
the heavy and light quarks, that destroys the picture of pair interactions and involves the additional 'centre-of- 
mass' point close to the centre of mass for the heavy-heavy system. 



There is another point of view to the problem of three-quark bound states in QCD, i.e. the 
representation of quark-gluon string. In the string-like picture of doubly heavy baryon shown in 
Fig.|TTT|, the above conclusions on the structure of mass spectra for the doubly heavy baryons derived 
on the basis of pair interactions should be essentially modified. Indeed, to the moment we have to 
introduce the centre of string, which is very close to the centre of mass for the doubly heavy diquark. 
Furthermore, the light quark interacts with the doubly heavy diquark as a whole, i.e. with the string 
tension identical to that in the heavy-light mesons Qq. Therefore, two different assumptions on the 
nature of interactions inside the doubly heavy baryons: pair interactions or string-like picture, result 
in a distinct variation of predictions on the mass spectra of these baryons for both the ground states 
and excitation levels. The only criterion testing the assumptions is provided by an experimental 
observation and measurements. 

In this review we follow the approximation of doubly heavy diquark, which is quite reasonable as 
we have clarified in the discussion given above. To enforce this point we refer to the consideration 
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of doubly heavy baryon masses in the framework of QCD sum rules (see Chapter 2), which result 
in the masses of ground states in a good agreement with the estimates obtained in the potential 
approach with the factorization of doubly heavy diquark. 

The qualitative picture for the forming of bound states in the system of (QQ'q) is determined 
by the presence of two scales of distances, which are given by the size of QQ'-diquark subsystem, 
tqqi, in the anti-triplet color state as well as by the confinement scale, Aqcd, for the light quark q, 
so that 

rQQ' ■ Aqcd < 1, Aqcd < mg. 

Under such conditions, the compact diquark QQ' looks like a static source approximated by the 
local colored QCD field interacting with the light quark. Therefore, we can use a set of reliable 
results in models of mesons with a single heavy quark, i.e. with a local static source belonging 
to the anti-triplet representation of SU(3)c group. The successful approaches are the potential 
models ^ and the Heavy Quark Effective Theory (HQET) |^ in the framework of expansion in 
the inverse heavy quark mass. We apply the nonrelativistic quark model with the potential by 
Buchmiiller-Tye |13|. Then theoretically we can talk on the rough approximation for the light 
quark. Indeed, since m^^^ <^ Aqcd the light quark is relativistic. Nevertheless, we introduce the 
system with a finite number of degrees of freedom and an instantaneous interaction V^(r). This fact 
is a disadvantage because the confinement supposes the following: a) the generation of sea around 
the light quark, i.e. the presence of infinite number of gluons and quark-antiquark pairs, and b) 
the nonperturbative effects with the correlation time tqcd ~ I/Aqcd, that is beyond the potential 
approach. However, phenomenologically the introduction of constituent mass m^^ ~ Aqcd as a 
basic parameter determining the interaction with the QCD condensates, allows us to successfully 
adjust the nonrelativistic potential model with a high accuracy {SM 30 — 40 MeV) by fitting the 
existing experimental data, that makes the approach to be quite a reliable tool for the prediction 
of masses for the hadrons, containing the heavy and light quarks. 

As for the diquark QQ', it is completely analogous to the heavy quarkonium QQ' except the 
very essential peculiarities. 

1. {QQ')% is a system with the nonzero color charge. 

2. For the quarks of the same flavor Q = Q' ii is necessary to take into account the Pauli principle 
for the identical fermions. 

The second item turns out to forbid the sum of quark spins S=0 for the symmetric, spatial parity P- 
even wave functions of diquark, \l'd(r) (the orbital momentum equals = 2n, where n = 0, 1, 2 . . . 
), as well as S=l is forbidden for the anti-symmetric, -odd functions ^E'(i(r) (i.e. = 2r;, + 1). The 
nonzero color charge leads to two problems. 

First, we cannot generally apply the confinement hypothesis on the form of potential (an infinite 
growth of energy with the increase of the system size) for the object under consideration. However, 
it is unpossible to imagine a situation, when a big colored object with a size r > 1/Aqcd has a 
finite energy of self-action, and, to the same moment, it is confined inside a white hadron (the 
singlet over SU(3)c) with r ~ 1/Aqcd due to the interaction with another colored source. In the 
framework of well-justified picture of the hadronic string, the tension of such string in the diquark 
with the external leg inside the baryons is only two times less than in the quark-antiquark pair 
inside the meson qq', and, hence, the energy of diquark linearly grows with the increase of its size. 
So, the effect analogous to the confinement of quarks takes place in the similar way. In the potential 
models we can suppose that the quark binding appears due to the effective single exchange by a 
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colored object in the adjoint representation of SU(3)c (the sum of scalar and vector exchanges is 
usually taken). Then, the potentials in the singlet {qq') and anti-triplet {qq') states differ by the 
factor of 1/2, that means the confining potential with the linear term in the QCD-motivated models 
for the heavy diquark (QQOsc- ■'■^ present chapter we use the nonrelativistic model with the 
Buchmiiller-Tye potential for the diquark, too. 

Second, in the singlet color state {QQ') there are separate conservations of the summed spin 
S and the orbital momentum L, since the QCD operators for the transitions between the levels 
determined by these quantum numbers, are suppressed. Indeed, in the framework of multipole 
expansion in QCD |T^, the amplitudes of chromo- magnetic and chromo-electric dipole transitions 
are suppressed by the inverse heavy quark mass, but in addition, the major reason is provided by 
the following: a) the necessity to emit a white object, i.e. at least two gluons, which results in the 
higher order in I/itlq, and b) the projection to a real phase space in a physical spectrum of massive 
hadrons in contrast to the case of massless gluon. Furthermore, the probability of a hybrid state, 
say, the octet sybsystem {QQ') and the additional gluon, i.e. the Fock state IQQ'^^g), is suppressed 
due to both the small size of system and the nonrelativistic motion of quarks (for a more strict 
consideration see ref. |1^). In the anti-triplet color state, the emission of a soft nonperturbative 
gluon between the levels determined by the spin Sd and the orbital momentum Ld of diquark, is 
not forbidden, if there are no some other no-go rules or small order-parameters. For the quarks 
of identical flavors inside the diquark, the Pauli principle leads to that the transitions are possible 
only between the levels, which either differ by the spin {ASd = 1) and the orbital momentum 
{ALd = 2?T, -|- 1), instantaniously, or belong to the same set of radial and orbital excitations with 
ALd = 2n. Therefore, the transition amplitudes are suppressed by a small recoil momentum of 
diquark in comparison with its mass. The transition operator changing the diquark spin as well as 
its orbital momentum, has the higher order of smallness because of either the additional factor of 
l/mg or the small size of diquark. These suppressions lead to the existence of quasi-stable states 
with the quantum numbers of Sd and Ld. In the diquark composed by the quarks of different 
flavors, be, the QCD operators of dipole transitions with the single emission of soft gluon are not 
forbidden, so that the lifetimes of levels can be about the times for the forming of bound states or 
with the inverse distances between the levels themselves. Then, we cannot insist on the appearance 
of excitation system for such the diquark with definite quantum numbers of the spin and orbital 
momentum^. 

Thus, in the present review we explore the presence of two physical scales in the form of factoriza- 
tion for the wave functions of the heavy diquark and light constituent quark. So, in the framework 
of nonrelativistic quark model the problem on the calculation of mass spectrum and characteristics 
of bound states in the system of doubly heavy baryon is reduced to two standard problems on the 
study of stationary levels of energy in the system of two bodies. After that, we take into account the 
relativistic corrections dependent of the quark spins in two subsystems under consideration. The 
natural boundary for the region of stable states in the doubly heavy system can be assigned to the 



threshold energy for the decay into a heavy baryon and a heavy meson. As was shown in [|T5|, the 
appearance of such threshold in different systems can be provided by the existence of an universal 
characteristics in QCD, a critical distance between the quarks. At distances greater than the critical 
separation, the quark-gluon fields become unstable, i.e. the generation of valence quark-antiquark 
pairs from the sea takes place. In other words, the hadronic string having a length greater than 
the critical one, decay into the strings of smaller sizes with a high probability close to unit. In the 



•^In other words, the presence of gluon field inside the baryon She leads to the transitions between the states with 
the different excitations of diquark, like \bc) \bcg) with ASd = 1 or ALd = 1, which are not suppressed. 
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framework of potential approach this effect can be taken into account by that we will restrict the 
consideration of excited diquark levels by the region, wherein the size of diquark is less than the 
critical distance, vqq/ < ~ 1.4 — 1.5 fm. Furthermore, the model with the isolated structure of 
diquark looks to be reliable, just if the size of diquark is less than the distance to the light quark 
rQQ' < ri. 

The peculiarity of quark-diquark picture for the doubly heavy baryon is the possibility of mixing 
between the states of higher diquark excitations, possessing the different quantum numbers, because 
of the interaction with the light quark. Then it is difficult to assign some definite quantum numbers 
to such excitations. We will discuss the mechanism of this effect. 

In Section 1.1 we describe a general procedure for the calculation of masses for the doubly heavy 
baryons in the framework of assumptions drawn above. We take into account the spin-dependent 
corrections to the potential motivated in QCD. The results of numerical estimates are presented in 
Section 1.2, and, finally, our conclusions are discussed in the end of this chapter. 



1.1. Nonrelativistic potential model 

As we have mentioned in the Introduction, we solve the problem on the calculation of mass spectra of 
baryons containing two heavy quarks, in two steps. First, we compute the energy levels of diquark. 
Second, we consider the two-body problem for the light quark interacting with the point-like diquark 
having the mass obtained in the first step. In accordance with the effective expansion of QCD in the 
inverse heavy quark mass, we separate two stages of such the calculations. So, the nonrelativistic 
Schrodinger equation with the model potential motivated by QCD, is solved numerically. After that, 
the spin-dependent corrections are introduced as perturbations suppressed by the quark masses. 



1.1.1. Potential 



The potential of static heavy quarks illuminates the most important features of QCD dynamics: the 
asymptotic freedom and confinement. In the leading order of perturbative QCD at short distances 
and with a linear confining term in the infrared region, the potential of static heavy quarks was 
considered in the Cornell model [1^, incorporating the simple superposition of both asymptotic 
limits (the effective coulomb and string- like interactions). The observed heavy quarkonia posed in 
the intermediate distances, where both terms are important for the determination of mass spectra 
(see Fig. |1.2| ). So, the phenomenological approximations of potential (logarithmic one [l^ and 
power law [1^), taking into account the regularities of such the spectra, were quite successful [Q. 

The quantities more sensitive to the global properties of potential are the wave functions at the 
origin as related to the leptonic constants and production rates. So, the potentials consistent with 
the asymptotic freedom to one and two loops as well as the linear confinement were proposed by 
Richardson |T9[, Buchmiiller and Tye |]13[, respectively. 

In QCD the static potential is defined in a manifestly gauge invariant way by means of the 
vacuum expectation value of a Wilson loop P0| , 

V{r) = - Urn ^ In(Wr) , 



T^oo 



iT 



trP exp 



Here, F is taken as a rectangular loop with time extension T and spatial extension r. The gauge 
fields An are path-ordered along the loop, while the color trace is normalized according to tr(..) = 
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V{r), GeV 
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Figure 1.2. The Corenell model of static potential and sizes of observed heavy quarkonia with charmed quarks 
(the family ■0) and bottom quarks (the family T). 

tr(..)/trl. This definition corresponds to the calculation of effective action for the case of two 
external sources fixed at a distance r during an infinitely long time period T, so that the time- 
ordering coincides with the path-ordering. Moreover, the contribution into the effective action by 
the path parts, where the charges have been separated to the finite distance during a finite time, 
can be neglected in comparison with the infinitely growing term of V{r) ■ T. Let us emphasize that 
the defined static potential is, by construction, the renormalization invariant quantity, since the 
action, by definition, doe not depend on the normalization point. 

Generally, one introduces the V scheme of QCD coupling constant by the definition of QCD 
potential of static quarks in momentum space as follows: 

V(c-) = -C,'-^^, (1.2) 

q 

so that for the such- way introduced value of one can derive some results at large virtualities in 
the perturbative QCD as well as at low transfer momenta in the approximation of linear term in 
the potential confining the quarks. 

In this section, first, we discuss two regimes for the QCD forces between the static heavy quarks: 
the asymptotic freedom and confinement. Then we follow the method by Buchmiiller and Tye and 
formulate how these regimes can be combined in a unified /? function for cty obeyed both limits of 
small and large QCD couplings. 




14 



1.1.2. Perturbative results at short distances 



Technically, using a given scheme of regularization, say, MS, one has to calculate the perturbative 
expansion for the potential of static quarks. This potential can be written down as the coulomb one 
with the running coupling constant in the so-called V scheme. Thus, the perturbative calculations 
provide us with the matching of MS scheme with V-one. The calculations with the n loop running 
of requires the n — 1 loop matching to a^. Note, that initial two coefficients of corresponding 
j3 functions are scheme and gauge independent, while others generally depend. The V scheme is 
defined for the observed quantity, that implies its j3 function to be gauge invariant. 
In the perturbative QCD the quantity can be matched with a-Ms 



a. 



oo ^ 
n=0 ^ 



47r 



n=0 



ftMsiq 



1.3) 



Two loop results for the j3 function and the one loop matching condition for the potential were 
available to the moment of Buchmiiller-Tye publication. Recently, the progress in calculations has 
provided us with the two loop matching of V and MS schemes |3T,^, that can be combined with 
the three loop running of af^. At present, in expansion ( |1.3| ) the coefficients of tree approximation 
ao, the one loop contribution ai and new results for the two-loop term 02 (see are known. 

After the introduction of a = the (3 function is actually defined by 



;i-4) 



n=0 



SO that = p^^l and /3J = Pf"" - ai/^^' + (aa - ai)(3^^ 

The Fourier transform results in the position-space potential ||21 



V{r) 



+ 



'^r/,2u 2 



4tt 



aMs(/i 
47r 



P^{W{nr') + —) + 2(/3i + 2poai) ln(/ir') + as 



;i.5) 



with r' = rexp(7£;). Defining the new running coupling constant, depending on the distance. 



(1.6) 



we can calculate its /3 function from ( |1.5|) , so that |2T 



TT 



P2 = P2 + y/5o^ 



;i-7) 



and the minor coefficients (3q-^ are equal to the values independent of the scheme. Note that 
the perturbative potential ( |1.5| ), by construction, is independent of normalization point, i.e. it is 
the renormalization group invariant. However, in the problem under consideration the truncation 
of perturbative expansion, wherein the coefficients do not decrease^, leads to a strong custodial 

"^Moreover, according to the investigations of renormalon, the coefficients in the series of perturbation theory for 
the potential increase in the factorial power, so that the series has a meaning of asymptotic one. 
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dependence on the normalization point. So, putting the normahzation point fi in the region of 
charned quark mass, we find that the two-loop potential with the three-loop running coupling 
constant af^ has an unremovable additive shift depending on /i. This shift has variation in wide 
limits. This fact illuminates the presence of infrared singularity in the coupling constant of QCD, 
so that the //-dependent shift in the potential energy has the form of a pole posed at Aqcd • 



Thus, in order to avoid the ambiguity of static potential in QCD we have to deal with infrared 
stable quantities. The motivation by Buchmiiller and Tye was to write down the (3 function of 
ay consistent with two known asymptotic regimes at short and long distances. They proposed 
the function, which results in the effective charge determined by two parameters, only: the per- 
turbative parameter is the scale in the running of coupling constant at large virtualities and the 
nonperturbative parameter is the string tension. The necessary inputs are the coefficients of (3 func- 
tion. The parameters of potential by Buchmiiller and Tye were fixed by fitting the mass spectra of 
charmonium and bottomonium Particularly, in such the phenomenological approach the scale 
A^=^ ^510 MeV was determined. It determines the asymptotic behaviour of coupling constant at 
latge virtualities in QCD. This value is in a deep contradiction with the current data on the QCD 
couphng constant af^ [^. In addition, one can easily find that the three- loop coefficient /9J for 
the (3 function suggested by Buchmiiller and Tye is not correct even by its sign and absolute value 
in comparison with the exact coefficient recently calculated in . 



Thus, the modification of Buchmiiller-Tye (BT) potential of static quarks as dictated by the 
current status of perturbative calculations is of great interest. 



To normalize the couplings in deep perturbative region, we use (|1.3| ) at = m 



1.1.3. The quark confinement 

The nonperturbative behaviour of QCD forces between the static heavy quarks at long distances r 
is usually represented by the linear potential (see discussion in ref. p5| ) 



l^™^f(r) = k-r, (1.8) 

which corresponds to the square-law limit for the Wilson loop. 

We can represent this potential in terms of constant chromo-electric field between the sources 
posed in the fundamental representation of SU(Ai'c). So, in the Fock-Schwinger gauge of fixed point 

■ A^^lx) = 0, we can represent the gluon field by means of strength tensor A^{x) ~ ^x'^Gfj,^{0), 
so that for the static quarks separated by the distance r we have Qi{0) Gmo(O) Qji.^) = 
where the heavy quark fields are normalized to unit. Then, the confining potential is written down 
as 

V'^"-'{r) = ^gsCFE-r. 

Supposing, that the same strength of the field is responsible for the forming of gluon condensate 
(see Fig. |L3|) , and introducing the colored sources rii, which have to be averaged in the vacuum, 
we can easily find 

{Gl)=ACj.E'{nn). 



For the linear term in the potential, the consideration in [24] leads to 



■■CfJ{-GI^). (1.9) 
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Figure 1.3. Under the action of charged sources, the vacuum chromo-electric field is aligned along the axis 
connecting the heavy quarks. 



The k term is usually represented through a parameter a'p as 

1 



27ra'p 



Following Buchmiiller and Tye, we put a'p = 1.04 GeV This value of tension, that is related 
with a slope of Regge trajectories, can be compared with the estimate following from (|1.9|) . At 
i^GlJ = (1.6 ± 0.1) ■ 10-2 GeV^ we have found 

a'p = 1.04 ±0.03 GeV"^ 

which is in a good agreement with the well known value of Regge trajectory slope. 



The form of (|1.8|) corresponds to the limit, when at low virtualities the coupling ay 

tends to 



f 2^ K 



so that 



In 

which gives the confinement asymptotics for the jSy function. 



^ -«v(q ), (1-10) 



1.1.4. Unified /3 function and potential 

Buchmiiller and Tye proposed the procedure for the reconstruction of (3 function in the whole region 
of charge variation by the known limits of asymptotic freedom to a given order in as and confine- 
ment regime. Generalizing their method, the /5pt function found in the framework of asymptotic 
perturbative theory (PT) to three loops, is transformed to the /? function of effective charge as 
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follows 



1 



/?pt(o) 

1 



+ 



/5J 



/3o 



/3(a) 



exp 



1 

/3oo 



+ 



/3o 



exp 



(1.11) 



where the exponential factor in the second term contributes to the next-to-next-to-leading order 
at a — 0. This function has the essential peculiarity at a — 0, so that the expansion is the 
asymptotic series in a. At a ^ oo the j3 function tends to the confinement limit represented in 
LlOp . Remember that the one and two- loop static potentials matched with the linear term of 



confinement lead to the contradiction with the value of QCD coupling constant extracted at the 
scale of Z boson mass if we fit the mass spectra of heavy quarkonia in such potentials. We will show 
that the static potential in the three-loop approximation results in the consistent value of QCD 
couphng constant at large virtualities. 

In the perturbative limit the usual solution for the running coupling constant 



a(/i^ 



/Soln^ 



Pi 



In ln-^+ 



PI In ^ 



A2 



PI 



/^oln^^ 



In In -— ■ — In In -— — 1 + 



A2 



A2 



PI 



is valid. Using the asymptotic limit of ( [L.12[ ), one can get the equation 

A2 



InC 



/3oa(/i2 



Pi /""^^^^ 
-2 ln/5oa(/i^) + / dx 
Po Jo 



Pi ^ 1 



PqX"^ PqX P{x 



[1.12) 



[1.13] 



which can be easily integrated out, so that we get an implicit solution for the charge depending on the 
scale. The implicit equation can be inverted by the iteration procedure, so that well approximated 
solution has the form 

a(^^) = ^ TT, (1-14) 



Poln (l+r^(^2)£ij 



where //(/x^) is expressed through the coefficients of perturbative P function and parameter /, which 
is related to the slope of Regge trajectories and the integration constant, the scale A, by the relation 



In Att^ C pa'p = In Po + 



Pi_ 
2Pi 



IE + 



/2 



2/5, 



P^Po-P! 

PI 



I 



:i.i5) 



which completely fixes the parameters of P function and the charge in terms of scale A and slope 



a'p. 



1.1.5. Setting the scales 

As we have already mentioned the slope of Regge trajectories, determining the linear part of poten- 
tial, is fixed as a'p 
and pose 



1.04 GeV~^ We use also the measured value of QCD coupling constant 



0.123, 
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as the basic input of the potential. Then, we evaluate q;v(^|) ~ 0.1306, and put it as the normal- 
ization point for a(m|) = a;v(^|)/(47r). Further, we find the values of A for the effective charge. 



depending on the number of active flavors ||2^. After determining the momentum space dependence 



of the charge, we perform the Fourier transform to get 



8C 

V{r) = k ■ r u{r), (1-16) 



with the function 



which is calculated numerically at r > 0.01 fm, while at short distances the behaviour of potential 
is purely perturbative, so that at r < 0.01 fm we put 

Vir) = -C,^^^^, (1.17) 



where the running q;v(1/t^) is given by eq. (|1.12|) with the appropriate value of /3J at n/ = 5, and with 
the matching with the potential ( |1.16| ) at = 0.01 fm, where we have found aYi^/f^s) — 0.22213, 
which implies A^^^g = 617.42 MeV. 

Thus, we have completely determined the potential of static heavy quarks in QCD with the three- 
loop running of coupling constant. In Fig. |L^we present the potential versus the distance between 
the quarks. As we can see the potential is very close to what was obtained in the Cornell model 
in the phenomenological manner by fltting the mass spectra of heavy quarkonia. We can draw the 
conclusion that accepting the normalization by the value of QCD coupling constant at the virtuality 
= m| and using the three-loop evolution for the effective charge incorporating the conflnement 
with the linear term, we have got the static potential consistent with the phenomenological models 
and, hence, with the calculations of mass spectra for the heavy quarkonia in the nonrelativistic 
approximation. 

Such the consistency of potential with the parameters of QCD has become possible due to the 
fact that in the two-loop approximation for the coulomb potential the calculations have led to the 
essential corrections to the /3 function of effective charge, so that A(3/f3 ~ 10%. This correction 
is important for the determination of critical values of charge, i.e. the value in the intermediate 
region between the perturbative and nonperturbative regimes. Moreover, the two-loop matching 
condition and the three-loop running of coupling constant normalized by the data at the high energy 
of mz determine the region of energetic scale for changing the regimes mentioned above. This scale 
strongly correlates with the data on the mass spectra of heavy quarkonia. Indeed, it is connected 
with the splitting of masses between the IS* and 25* levels. We stress that the two-loop improvement 
gives the correct normalization of effective coulomb exchange at the distances characteristic for the 
average separation between the heavy quarks inside the heavy quarkonium and determines the 
evolution at short distances r < 0.08 fm, that is important in the calculations of leptonic constants 
related with the wave functions at the origin. 

The analysis of potential for the static quarks in the calculations of both the mass spectra for 
the heavy quarkonia and the leptonic constants for the vector states is presented in |2^, where the 
heavy quark masses are determined in the potential approach. 

We emphasize that the potential by Buchmiiller and Tye was obtained under fltting the experi- 
mental mass spectra for the heavy quarkonia, and it is numerically very close to the static potential 
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Figure 1.4. The potential of static heavy quarks in QCD (solid line) in comparison with the Cornell model 
(dashed line) (up to an additive shift of energy scale). 

under consideration consistent with the normahzation of QCD couphng constant at large viruah- 
ties. Therefore, the Buchmiiller-Tye potential continues to keep its phenomenological usage for the 
calculations of mass levels for the hadrons with c and b quarks with the accuracy about 40-70 MeV, 
that is a systematic uncertainty of potential approach. 



1.1.6. System of levels 



Following ref. in the potential model we use the Buchmiiller-Tye ansatz, which takes into 
account the coulomb corrections at short distances with the running coupling constant in two loops, 
while at large distances the interaction energy linearly increases, that provides the confinement. In 
the anti-triplet quark state we introduce the factor of 1/2 because of the color structure of bound 
quark-quark state. For the interaction of diquark with the light constituent quark, the corresponding 
factor is equal to unit. 

As was shown in , the nonperturbative constituent term introduced into the mass of nonrela- 
tivistic quark, exactly coincides with the additive constant, subtracted from the coulomb potential. 

Thus, we extract the masses of heavy quarks by fitting the real spectra of charmonium and 
bottomonium, 

nic = 1.486 GeV, = 4.88 GeV, (1.18) 

so that the mass of the level in the heavy quarkonium has been calculated as, say, M(c) = 2mc -|- E, 
where E is the energy of stationary Schrodinger equation with the model potential V. Then, we have 
supposed that the mass of meson with a single heavy quark is equal to M{Qq) = ttiq + rriq + E, and 
E = (T) -|- (V — SV) , whereas the additive term in the potential is introduced because the constituent 



mass of light quark is determined as a part of interaction energy SV 



m, 



In accordance with 
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fitting the masses of lieavy mesons, we get ruq = 0.385 GeV. 

Table 1.1. The spectrum of 66-diquark levels without spin-dependent splittings: masses and mean-squared 
radii. 



diquark level 


, GeV 


{r'y/'', fm 


diquark level 


, GeV 


(r^)i/^ fm 


IS 


9.74 


0.33 


2P 


9.95 


0.54 


2S 


10.02 


0.69 


3P 


10.15 


0.86 


3S 


10.22 


1.06 


4P 


10.31 


1.14 


4S 


10.37 


1.26 


5P 


10.45 


1.39 


5S 


10.50 


1.50 


6P 


10.58 


1.61 


3D 


10.08 


0.72 


4D 


10.25 


1.01 


5D 


10.39 


1.28 


6D 


10.53 


1.51 


4F 


10.19 


0.87 


5F 


10.34 


1.15 


6F 


10.47 


1.40 


5G 


10.28 


1.01 


6G 


10.42 


1.28 


6M 


10.37 


1.15 



Table 1.2. The spectrum of 6c-diquark levels without spin-dependent splittings: masses and mean-squared 
radii. 



diquark level 


, GeV 


(r2)V2^ fm 


diquark level 


, GeV 


(r2)i/2^ fm 


IS 


6.48 


0.48 


3P 


6.93 


1.16 


2S 


6.79 


0.95 


4P 


7.13 


1.51 


3S 


7.01 


1.33 


3D 


6.85 


0.96 


2P 


6.69 


0.74 


4D 


7.05 


1.35 


4F 


6.97 


1.16 


5F 


7.16 


1.52 


5G 


7.09 


1.34 


6H 


7.19 


1.50 



The results of calculations for the energy levels in the Schrodinger equation with the Buch- 
miiller-Tye potential for the various diquark systems are presented in Tables p..l| - p..3| , while the 
characteristics of corresponding wave functions are shown in Tables |1.4H L6. 



We have checked that with a good accuracy the binding energy and the wave function of light 
quark practically do not depend on the flavors of heavy quarks. Indeed, large values of diquark 
masses give small contributions into the reduced masses. This fact leads to small corrections to the 
wave functions in the Schrodinger equation. 

So, for the states lying below the threshold of doubly heavy baryon decay into the heavy baryon 
and heavy meson, the energies of levels of light constituent quark are equal to 

E{ls) = 0.38 GeV, E{2s) = 1.09 GeV, E{2p) = 0.83 GeV, 

where the energy has been defined as the sum of light quark constituent mass and eigen-value of 
Schrodinger equation. In HQET the value of A = E[ls) is generally introduced. Then we can draw 
a conclusion that our estimate of A is in a good agreement with calculations in other approaches. 
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Table 1.3. The spectrum of cc-diquark levels without spin-dependent splittings: nnasses and mean-squared 
radii. 



diquark level 


, GeV 


(r2)V2^ fm 


diquark level 


, GeV 


(r2)V2^ fm 


IS 


3.16 


0.58 


3P 


3.66 


1.36 


2S 


3.50 


1.12 


4P 


3.90 


1.86 


3S 


3.76 


1.58 


3D 


3.56 


1.13 


2P 


3.39 


0.88 


4D 


3.80 


1.59 



Table 1.4. The characteristics of radial wave function for the fe6-diquark: Rd{ns){'^) (GeV^/^), R'^^^^^iO) (GeV 

5/2^ 



iiL 




nL 




IS 


1.346 


2P 


0.479 


2S 


1.027 


3P 


0.539 


3S 


0.782 


4P 


0.585 


4S 


0.681 


5P 


0.343 



Table 1.5. The characteristics of radial wave function for the 6c-diquark: Rd{ns){'^) (GeV^/^), -R^(„p)(0) (GeV 

5/2y 



iiL 




nL 




IS 


0.726 


2P 


0.202 


2S 


0.601 


3P 


0.240 



Table 1.6. The characteristics of radial wave function for the cc-diquark: Rd(ns)i'^) (GeV^/^), -R^(„p)(0) (GeV 

5/2 y 



iiL 


(0) 


iiL 


-^:/(„.)(o) 


IS 


0.530 


2P 


0.128 


2S 


0.452 


3P 


0.158 



22 



This fact confirms the rehabihty of such the phenomenological predictions. For the hght quark 
radial wave functions at the origin we find 

R^s{0) = 0.527 GeV^/^ R2s{0) = 0.278 GeV^/^ ^2p(0) = 0.127 GeV^/^ 

The analogous characteristics of bound states of the c-quark interacting with the 66-diquark, are 
equal to 

E{ls) = 1.42 GeV, E{2s) = 1.99 GeV, E{2p) = 1.84 GeV, 
with the wave functions 

Ris{0) = 1.41 GeV^/^ i?2s(0) = 1.07 GeV^/^ ^2p(0) = 0.511 GeV^^l 

For the binding energy of strange constituent quark we add the current mass ^ 100 — 150 
MeV. 

1.1.7. Spin-dependent corrections 



According to we introduce the spin-dependent corrections causing the splitting of nL-levels 
of diquark as well as in the system of light constituent quark and diquark (ra = + L + 1 is the 
principal number, is the number of radial excitation, L is the orbital momentum). For the heavy 
diquark containing the identical quarks we have 

^ 2 I 2m% \ rdr 3 



+ + Asqi ■ SQ,[inS{r)] (1.19) 

3 rriQ r-^ 3 orriQ 

-l^s^ ,^ I J 6(Ld ■ Sd)' + 3(Ld ■ Sd) - 2Ld'Sd']4, 
3 m^4Ld -3 

where Ld, Sd are the orbital momentum in the diquark system and the summed spin of quarks 
composing the diquark, respectively. Taking into account the interaction with the light constituent 
quark gives (S = Sd + Si) 



^^^^ 4 I 2m2 ^ 2mf \ rdr ^ 3 



1 1 (L ■ Sd + 2L ■ Si) , 4 ^1 ^ QM ^1 nom 

-as z h -as- (Sd + Ld) ■ Si 47r(5(r)J (1.20) 

3 rnqmi r-^ 3 onigmi 



-\c^s ^ J 6(L ■ Sr + 3(L ■ S) - 2L'S' 

3 mQiTLi AL^ — 3 

-6(L-Sd)'-3(L-Sd) + 2L2Sd2]^, 

where the first term corresponds to the relativistic correction to the effective scalar exchange, and 
other terms appear because of corrections to the effective single-gluon exchange with the coupling 
constant as- 
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The value of effective parameter ag can be determined in the following way. The splitting in the 
yS-wave heavy quarkonium {Q\Q2) is given by the expression 

AMM = las-^\Rnsm\ (1.21) 

9 171x7712 

where Rns{'>^) is the radial wave function of quarkonium. From the experimental data on the system 
of cc 

AM(15, cc) - 117 ±2 MeV, (1.22) 

and i?is(0) calculated in the model, we can determine as(^)- 

Let us take into account the dependence of this parameter on the reduced mass of the system, 
II. In the framework of one-loop approximation for the running coupling constant of QCD we have 

4:77 

whereas 6=11 — 2n//3 and n/ = 3 at < m^. From the phenomenology of potential models we 
know that the average kinetic energy of quarks in the bound state practically does not depend on 
the flavors of quarks, and it is given by the values 

(Ta) ^ 0.2 GeV, (1.24) 

(Ti) ^ 0.4 GeV, (1.25) 

for the anti-triplet and singlet color states, correspondingly. Substituting the definition of the 
nonrelativistic kinetic energy 

(T) = (1.26) 

we get 

47r 

= , |^2 y (1-27) 

b ■ \n{2{T}ii/AgcD) 

whereas numerically Aqcd ~ 113 MeV. 

For the identical quarks inside the diquark, the scheme of L5'-coupling well known for the 
corrections in the heavy quarkonium, is applicable. Otherwise, for the interaction with the light 
quark we use the scheme of jj-coupling (here, LSi is diagonal at the given Ji, (Ji = L+Si, J = Ji-|-J), 
where J denotes the total spin of baryon, and J is the total spin of diquark, J = Sd + Ld)- 

Then, to estimate various terms and mixings of states, we use the transformations of bases (in 
what follows S = Si -I- J) 

I J; = x:(-i)^'^^'^^"''V(2^ + mji + 1) { ^ J i^}'-^'^^ (^-28) 

and 

I J; Ji) = J](-l)(^"+^'+^+-^) V(2J<i + 1)(2J, + 1) { J t}'-^;-^'^)- (1-29) 

Thus, we have defined the procedure of calculations for the mass spectra of doubly heavy baryons. 
This procedure leads to results presented in the next section. 
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1.2. Numerical results 



In this Section we present the resuhs on the mass spectra with account for the spin-dependent 
sphtting of levels. As we have clarified in the Introduction, the doubly heavy baryons with identical 
heavy quarks allow quite a reliable interpretation in terms of diquark quantum numbers (the summed 
spin and the orbital momentum). Dealing with the excitations of 6c-diquark, we show the results 
on the spin-dependent splitting of the ground IS-state, since the emission of soft gluon breaks the 
simple classification of levels for the higher excitations of such diquark. 

For the doubly heavy baryons, the quark-diquark model of bound states obviously leads to the 
most reliable results for the system with the larger mass of heavy quark, i.e. for S^b. 

For the quantum numbers of levels, we use the notations ridLdUili, i.e. we show the value of 
principal quantum number of diquark, its orbital momentum by a capital letter and the principal 
quantum number for the excitations of light quark and its orbital momentum by a lower-case letter. 



The splitting of Ehi, baryon levels was in detail considered in [26|. The states with the total spin 



J = I (or i), can have different values of J/, and, hence, they have a nonzero mixing, when we 
perform the calculations in the perturbation theory built over the states with the definite total 
momentum Ji of the light constituent quark. For ^ = | the mixing matrix can be approximated by 
the diaginal matrix with a high accuracy. For J = i the mixing of states with the different values 
of total spin-orbital momentum of light quark is strong. The analysis for the lS2p and 2S2p levels 
in the system Ehf, was done in ref. [^, where one can see that the difference between the wave 



functions because of the slow change of diquark subsystem mass is not essential within the accuracy 
of the method. 

The splittings of D and G levels in the diquark are less than 11 MeV, so that these corrections 
are small for the diquark excitations having the sizes less than the distance to the light quark, i.e. 
for the states with low value of principal quantum number, under the systematic uncertainty about 
^ 30 - 40 MeV. 

For the hyper-fine spin-spin splitting in the system of quark-diquark, we have 



J(J+ 1) - J(J+ 1) - - as{2fiT)^—\Ri{0)\'\ (1.30) 

rribmi 



where -R;(0) is the radial wave function at the origin for the light constituent quark, and for the 
analogous shift of diquark level, we find 



A^;^y = -as{2fxT) — \R,m\ (1.31) 

The mass spectrum of and baryons is shown in Fig. |1.5| , wherein we restrict ourselves 
by the presentation of S-, P- and D-wave levels, while the table containing the numerical values of 
masses for the Ebb baryons is presented in [P^ . 

We can see in Fig. [L^ that the most reliable predictions are the masses of baryons ISls {J^ = 
3/2+, 1/2+), 2Pls (J^ = 3/2-, 1/2-) and 3Dls (J^ = 7/2+, . . .1/2+). The 2Pls-level is quasi- 
stable, because the transition into the ground state requires the instantaneous change of both 
the orbital momentum and the summed spin of quarks inside the diquark. The analogous kind 
of transitions seems to be the transition between the states of ortho- and para-hydrogen in the 
molecule of This transition take place in a non- homogeneous external field due to the magnetic 
moments of other molecules. For the transition of 2Pls 15*15, the role of such the external field 
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11.5 



3S2p 

^^li At B threshold 

1^2. 5P1. ^ Jm^ 

— Sm. 2P2p 

iPls 



, ADls 

3^1^ lS2p l^/l- 



10.0 



3P1. 3/2- 1/2^ 

10.5 ^ ^12- ? ^ 3m. ^ Si 

2S'ls ^ — — 1/2+ 

?J'2+ 2P1S 3/2- S/2+ 



1/2- 



^■^^■^ 3/2+ 

1/2+ 



Figure 1.5. The spectrum of baryons, containing two fe-quarks: H^^ and H^^, with account for the spin-dependent 
splittings of low-lying excitations. The masses are given in GeV. 

is played by the non-homogeneous chromo-magnetic field of the light quark. The corresponding 
perturbation has the form 

5V ~ — [Si-Hi + S2-H2-(Si + S2)-(H)] 

= :r^(V • ra) {S^ - S^) • H ^ — ^ (Si - S^) • Ji /(n), 

where f{ri) is a dimensionless nonperturbative function depending on the distance between of 
the light quark and diquark. The 5V operator changes the orbital momentum of light quark, 
too. It results in the mixing between the states with the same values of J^. If the splitting 
is not small (for instance, 2Pls — lS2p, where AE ~ Aqcd), then the mixing is suppressed as 
SV/AE ~ mgrn V^'Ke ^ ^incc thc admixture of 15*2^ in the 2Pls-statc is low, the 2Pls-lcvcls 
are quasi-stable, i.e. their hadronic transitions into the ground state with the emission of 7r-mesons 
are suppressed as we have derived, though an additional suppression is given by a small value of 
phase space. Therefore, we have to expect the presence of narrow resonances in the mass spectra 
of pairs Sf,f,7r, as they are produced in the decays of quasi-stable states with = 3/2~, l/2~. The 
experimental observation of such levels could straightforwardly confirm the existence of diquark 
excitations and provide the information on the character of dependency in f{ri), i.e. on the non- 
homogeneous chromo-magnetic field in the nonperturbative region. 

Sure, the 3Dls — 7/2+, 5/2+ states are also quasi-stable, since in the framework of multipole 
expansion in QCD they transform into the ground state due to the quadrupole emission of gluon 
(the E2-transition with the hadronization gq q''rr). 
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As for the higher excitations, the 3Pls-states are close to the lS'2p-levels with = 3/2^, 1/2^, 
so that the operators changing both the orbital momentum of diquark and its spin, can lead to the 
essential mixing with an amplitude 5y„„//Ai?„„/ ~ 1, despite of suppression by the inverse heavy 
quark mass and small size of diquark. We are sure that the mixing slightly shifts the masses of 
states. The most important effect is a large admixture of lS2p in 3Pls. It makes the state to 
be unstable because of the transition into the ground IS'ls-state with the emission of gluon (the 
El-transition). This transition leads to decays with the emission of vr-mesonsj^. 

The level 15*2^ = 5/2~ has the definite quantum numbers of diquark and light quark motion, 
because there are no levels with the same values of in its vicinity. However, its width of transition 
into the ground state and vr-meson is not suppressed and seems to be large, F ~ 100 MeV. 

For the transitions we have 

3- 3+ 3- 1+ 

- ^ - TT in S — wave, - — vr in D — wave, 
2 2 2 2 

1~ 3"'" 1~ l"*" 

- — i> - TT in D — wave, - — - vr in S — wave. 
2 2 '22 

The D-wave transitions are suppressed by the ratio of low recoil momentum to the mass of baryon. 

The width of state = 3/2"*" is completely determined by the radiative electromagnetic Ml- 
transition into the ground = 1/2+ state. 

The calculation procedure described above leads to the results for the doubly charmed baryons 
as presented in Table [L 



Table 1.7. The mass spectrum of H++ and baryons. 



{ndLdUiLi), 


mass, GeV 


{ndLdUiLi), 


mass, GeV 


(IS ls)l/2+ 


3.478 


(3P ls)l/2- 


3.972 


(IS ls)3/2+ 


3.61 


(3D ls)3/2'+ 


4.007 


(2P ls)l/2- 


3.702 


(IS 2p)3/2'- 


4.034 


(3D ls)5/2+ 


3.781 


(IS 2p)3/2- 


4.039 


(2S ls)l/2+ 


3.812 


(IS 2p)5/2- 


4.047 


(3D ls)3/2+ 


3.83 


(3D ls)5/2'+ 


4.05 


(2P ls)3/2- 


3.834 


(IS 2p)l/2'- 


4.052 


(3D ls)l/2+ 


3.875 


(3S ls)l/2+ 


4.072 


(IS 2p)l/2- 


3.927 


(3D ls)7/2+ 


4.089 


(2S ls)3/2+ 


3.944 


(3P ls)3/2- 


4.104 



As we have already mentioned, the heavy diquark composed of the quarks of different flavors, 
turns out to be unstable under the emission of soft gluons. So, in the Fock state of doubly heavy 
baryon, there is a sizable nonperturbative admixture of configurations including the gluons and 
diquark with the various values of its spin 5*^ and orbital momentum 

\Bbc,) = 0B|6cJ'^^ q) + i/i|6cff (7, q) + H2\h4'j'''^\ g,q) + ..., 
^Remember, that the Sgg'-baryons are the iso-dublets. 
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whereas the amphtudes of Hi, H2 are not suppressed with respect to O^. In the heavy quarkonium, 
the analogous operators for the octet-color states are suppressed by the probability of emission by 
the nonrelativistic quarks inside a small volume determined by the size of singlet-color system of 
heavy quark and anti-quark. In the baryonic system under consideration, a soft gluon is restricted 
only by the ordinary scale of confinement, and, hence, there is no suppression. 

We suppose that the calculations of masses for the excited S^c baryons are not so justified in 
the given scheme. Therefore, we present only the result for the ground state with = 1/2+ 

M=/ =6.85GeV, = 6.82 GeV, 

whereas for the vector diquark we have assumed that the spin-dependent splitting due to the 
interaction with the light quark is determined by the standard contact coupling of magnetic moments 
for the point-like systems. The picture for the baryon levels with no account for the spin-dependent 



perturbations suppressed by the heavy quark masses is shown in [EH] 



1.2.1. The doubly heavy baryons with the strangeness ^qq' 

In the leading approximation, we suppose that the wave functions and the excitation energies 
of strange quark in the field of doubly heavy diquark repeat the characteristics for the analogous 
baryons containing the ordinary quarks u, d. Therefore, the level system of baryons VLqqi reproduces 
that of Sqq/ up to an additive shift of the masses by the value of current mass of strange quark, 
nis ^ M{D,) - M{D) ^ M{Bs) - M{B) ^ 0.1 GeV. 

Further, we suppose that the spin-spin splitting of 2Pls and 3Dls levels of ^qq' is 20-30% 
less than in Eqq> (the factor of mu^/ms). As for the lS'2p-level, the procedure described above 
can be applied. So, for the matrix of mixing for the states with the different values of total 
momentum Ji practically can be assigned to be diagonal. This fact means that the following term 
of perturbation is dominant: 

1 / 2L- SA / dV{r) ^ 8^ l\ 
4 \ 2mf J \ rdr 3 ^ J 

Therefore, we can think that the splitting of 15*2^ is determined by the factor of Tn^^/rn^ with 
respect to the splitting of corresponding Sftb, i.e. it is 40% less than in S^b. Hence, the splitting is 
very small. 

For the baryon f2, the factor of ms/mc is not small. Hence, for 15*2^, the mixing matrix is not 
diagonal, so that the arrangement of 15*2^ states of VL can be slightly different from that of S. 

The following peculiarity of VLqqi is of great interest: the low- lying S- and P-excitations of 
diquark are stable. Indeed, even after taking into account the mixing of levels, a gluon emission 
makes a hadronization into the i^'-meson (the transitions of ^qqi '^qq' + K), while a single 
emission of vr-meson is forbidden because of the conservation of iso-spin and strangeness. The 
hadronic transitions with kaons are forbidden because of insufficient splitting between the masses 
of VLqqi and Sgg/. The decays with the emission of pion pairs belonging to the iso-singlet state, are 
suppressed by a small phase space or even forbidden. Thus, the radiative electromagnetic transitions 
into the ground state are the dominant modes of decays for the low- lying excitations of ^qq'- 
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1.2.2. Qbbc baryons 



In the framework of quark-diquark picture, we can build the model for the baryons containing 
three heavy quarks, bbc. However, as we estimate, the size of diquark turns out to be compa- 
rable with the average distance to the charmed quark. So, the model assumption on the com- 
pact heavy diquark cannot be quite accurate for the calculations of mass levels in this case. The 
spin-dependent forces are negligibly small inside the diquark, as we have already pointed out 
above. The spin-spin splitting of vector diquark interacting with the charmed quark, is given 
by A(ls) = 33 MeV, A(2s) = 18 MeV. For lS2p, the level shifts are small. So, for the state 
= 1/2^ we have to add the correction of —33 MeV. For the 3Dls-state the splitting is deter- 
mined by the spin-spin interaction. The characteristics of excitations for the charmed quark in the 
model with the potential by Buchmiiller and Tye have been presented above. Finally, we obtain 
the picture of flbbc levels presented in Table pT8| . 



Table 1.8. The mass spectrum of il^j^^ baryons. 



{ndLdUiLi), 


mass, GeV 


{ndLdUiLi), 


mass, GeV 


(IS ls)l/2+ 


11.12 


(3D ls)3/2'+ 


11.52 


(IS ls)3/2+ 


11.18 


(3D ls)5/2'+ 


11.54 


(2P Is) 1/2- 


11.33 


(IS 2p)l/2- 


11.55 


(2P ls)3/2- 


11.39 


(3D ls)7/2+ 


11.56 


(2S ls)l/2+ 


11.40 


(IS 2p)3/2'- 


11.58 


(3D ls)5/2+ 


11.42 


(IS 2p)3/2- 


11.58 


(3D ls)3/2+ 


11.44 


(IS 2p)l/2'- 


11.59 


(3D ls)l/2+ 


11.46 


(IS 2p)5/2- 


11.59 


(2S ls)3/2+ 


11.46 


(3P ls)3/2- 


11.59 


(3P ls)l/2- 


11.52 


(3S ls)l/2+ 


11.62 



Further, the excitations of ground fi^^^ state can strongly mix with large amplitudes because of 
small splittings between the levels, but they have small shifts of masses. This effect takes place for 
3Pls - lS2p with = 1/2", 3/2", and for 2^1s - 3Dls with = 1/2+, 3/2+. We suppose the 
prediction to be quite rehable for the states of ISls with = 1/2+, 3/2+, 15*2^ with = 5/2^ 
and 3Dls with = 5/2+, 7/2+. For these excitations, we might definitely predict the widths 
of their radiative electromagnetic transitions into the ground state in the framework of multipole 
expansion in QCD. The widths for the transitions will be essentially determined by the amplitudes 
of admixtures, which have a strong model dependence. Therefore, the experimental study of elec- 
tromagnetic transitions in the family of baryons could provide a significant information on the 
mechanism of mixing between the different levels in the baryonic systems. The electromagnetic 
transitions combined with the emission of pion pairs, if not forbidden by the phase space, saturate 
the total widths of excited levels. The characteristic value of total width is about F ~ 10 — 100 
keV, in the order of magnitude. 

Thus, the system of fi^^^ can be characterized by a large number of narrow quasi-stable states. 
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1.3. Discussion 



In this paper we have calculated the spectroscopic characteristics of baryons containing two heavy 
quarks, in the model with the quark-diquark factorization of wave functions. We have explored the 
nonrelativistic model of constituent quarks with the potential by Buchmiiller and Tye. The region 
of applicability of such the approximations has been pointed out. 

We have taken into account the spin-dependent relativistic corrections to the potential in the 
subsystems of diquark and light quark-diquark. Below the threshold of decay into the heavy baryon 
and heavy meson, we have found the system of excited bound states, which are quasi-stable under 
the hadronic transitions into the ground state. We have considered the physical reasons for the 
quasi-stability taking place for the baryons with two identical quarks. In accordance with the Pauli 
principle, the operators responsible for the hadronic decays and the mixing between the levels, are 
suppressed by the inverse heavy quark mass and the small size of diquark. This suppression is caused 
by the necessity of instantaneous change in both the spin and the orbital momentum of compact 
diquark. In the baryonic systems with two heavy quarks and the strange quark, the quasi-stability 
of diquark excitations is provided by the absense of transitions with the emission of both a single 
kaon and a single pion. These transitions are forbidden because of small splitting between the levels 
and the conservation of iso-spin and strangeness. 

The characteristics of wave functions can be used in calculations of cross sections for the doubly 
heavy baryons in the framework of quark-diquark approximation. 

The quark-diquark factorization in calculating the masses of ground states for the baryon systems 
with two heavy quarks was also considered in ref . , where the quasi-potential approach ||30| was 
explored. There is a numerical difference in the choice of heavy quark masses, that leads to that 
in 1^ the mass of doubly charmed diquark, say, about 100 MeV greater that the mass used in 
the above calculations. This difference determines the discrepacy of estimates for the masses of 
ground states presented in this paper and in |^. We beheve that this deviation between the 
quark masses is caused by the use of Cornell potential with the constant value of effective coulomb 
exchange coupling in contrast to the above consideration with the running coupling constant, that 
cancels the uncertainty in the arbitrary additive shift of energy. Furthermore, in the potential 
approach the masses of heavy quarks depend on the mentioned additive shift, which adjusted in the 
phenomemological models by comparing, say, the leptonic constants of heavy quarkonium calculated 
in the model with the values known from experiments. In the QCD motivated potential such the 
ambiguity of potential because of the additive shift is absent, so that the estimates of heavy quark 
masses have less uncertainities. Let us stress that in the Cornell model the leptonic constants 
were calculated by taking into account the one-loop corrections caused by the hard gluons. This 
correction is quite essential, in part, for the charmed quarks. The two-loop corrections are also 
important for the consideration of leptonic constants in the potential approach |^|. Moreover, 
in the constituent mass of light quark is posed with no correlation with the normalization of 



potential, while we put the constituent mass to be a part of nonperturbative energy in the potential. 
This can lead to an additional deviation between the estimates of baryon masses about 50 MeV. 
Taking into account the above notes on the systematic differences, we can claim that the estimates 



of ground states masses for the baryons with two heavy quarks in ||29| agree with the values obtained 
in the presented approach (see Table |1.9|) . 

In ref. [^], following in the framework of quasi-potential approach, the analysis of spin- 
dependent relativistic corrections was performed so that the overestimated, to our opinion, value of 
heavy diquark from |29] was used. Unfortunately, there is an evident mistake in the description of 
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Table 1.9. The masses of ground states M (in GeV) for the baryons with two heavy quarks calculated in various 
approaches (* denotes the results of authors in this review). The accuracy of predictions under the variation of 
model parameters is about 30-50 MeV. The systematic unceratinties are discussed in the text. 



baryon 


* 


|31| 


^9| 






H 


[p5l 




3.48 


3.74 


3.66 


3.66 


3.61 


3.65 


3.71 


^cc 


3.61 


3.86 


3.81 


3.74 


3.68 


3.73 


3.79 




3.59 


3.76 


3.76 


3.74 


3.71 


3.75 


3.89 




3.69 


3.90 


3.89 


3.82 


3.76 


3.83 


3.91 


'^bb 


10.09 


10.30 


10.23 


10.34 


- 


- 


10.43 




10.13 


10.34 


10.28 


10.37 


- 


- 


10.48 


^^66 


10.18 


10.34 


10.32 


10.37 






10.59 


"66 


10.20 


10.38 


10.36 


10.40 






10.62 


^cb 


6.82 


7.01 


6.95 


7.04 






7.08 


^cb 


6.85 


7.07 


7.00 


6.99 






7.10 


^cb 


6.90 


7.10 


7.02 


7.06 






7.13 


^cb 


6.91 


7.05 


7.05 


7.09 






7.23 


Kb 


6.93 


7.11 


7.09 


7.06 






7.24 


Kb 


6.99 


7.13 


7.11 


7.12 






7.27 



calculations in |^ , because both the parameter giving the relative contribution of scalar and vector 
part in the potential and the anomalous chromo-magnetic moment of heavy quark are denoted by 
the same symbol, that leads to numerical errors, since in |^ it was shown that these quantities have 
different values. This mistake enlarges the unceratinty about 100 MeV into th estimates of |^ 



so 



that we can consider that the results of ||3T| do not contradict with the presented description (see 
Table [LI). 

The estimates based on the hypothesis of pair interactions were presented in ref. p2[) so that 



in the light of discussion given in the beggining of this chapter the difference about 200-300 MeV, 
that follows from values in Table L9, is not amazing. This deviation is, in general, related with the 
different character of interquark forces in the doubly heavy baryon, though the uncertainty in the 
heavy quark masses is also important. 

In ref. simple speculations based on the HQET with the heavy diquark were explored, so 
that the estimates depend on the supposed mass of diquark composed of two heavy quarks. In this 
way, if we neglected the binding energy in the diquark, that is evidently related with the choice of 
heavy quark masses, then we got the estimates of ground state masses shown in Table 0|. 



Finally, in |35[ the analysis given in ^6[ was modified on the basis of interpolation formulae for 
the mass of ground state with account for the dependence of spin forces on both the wave functions 
and the effective coupling constant, which were changed under the quark contents of hadrons. In 
this way, the parameter of energy shift enters the fitting function, so that this parameter essentially 
changes under the transition from the description of mesons to baryons: 6m ~ 80 MeV — > 6b ~ 210 
MeV. This shift of energy provides a good agreement of fitting with the mass values for the mesons 
and baryons observed experimentally. However, if we suggest that the doubly heavy baryon is 
similar with the meson containing the local heavy source in the picture of strong interactions, then 
we should use the energy shift prescribed to the heavy mesons but the heavy baryons, wherein the 
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presense of system with two light quarks leads to the essential difference in the calculation of bound 
state masses, hence, to the energy shift different from the mesonic one. Such the substitution of 
parameters would lead to a more good agreement between the results of (see Table [L^) and 
the values obtained in this review. 

Summarizing, we can claim that, first of all, in the framework of potential approach in the 
calculations of masses for the doubly heavy baryons the dominant uncertainty is caused by the 
choice of heavy quark masses, so that due to the adjustment on the systems with heavy quarks, the 
analysis presented in the QCD-motivated model of potential with the running coupling constant at 
short distances and the linear nonperturbative term confining quarks at large distances, gives the 
most reliable predictions. 

A new field of interest for the investigations is radiative, electromagnetic or hadronic, transitions 
between the quasi-stable states in families of baryons with two heavy quarks. A first step in the 
study of this problem was recently done in wherein some preliminary results were obtained on 
the electromagnetic transitions between the levels of Sf,c. 
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Chapter 2. Sum rules of nonrelativistic QCD: 
two-point correlators 



In the framework of potential models in Chapter 1 we have described the families of baryons with 
two heavy quarks, which form a set of narrow excited states in addition to the ground states, so 
that the mass spectrum is very similar with the system of levels in the heavy quarkonium. In the 
method of QCD sum rules [0 for the two-point correlators of baryonic currents, the masses and 



coupling constants for the baryons with two heavy quarks were calculated in |3^. However, the 
analysis in has some disadvantages connected to an unstable divergency of sum rules in the 
range of parameters determining the baryonic currents. This leads to quite large uncertainties in 
calculations. 

In this chapter we investigate the sum rules of NRQCD for the two-point correlators of currents 
corresponding to the baryons with two heavy quarks. The main physical argument of such the 
consideration is the nonrelativistic motion of heavy quarks inside the diquark of small size, that 
interacts with the light quark. This fact leads to some definite expressions for the structure of 
baryonic currents written down in terms of nonrelativistic heavy quark fields. In the leading order 
of inverse heavy quark mass and relative velocity of heavy quarks inside the diquark, we have to 
take into account the hard gluon corrections in order to derive the relations between the correlators 
of nonrelativistic quarks in NRQCD and the correlators in full QCD. The corresponding anomalous 



dimensions for the baryonic currents were calculated up to two loops in ref. The structure of 
currents in NRQCD corresponds to some choice of parameters in expressions of full QCD. These 
values of parameters are posed in the range of significant uncertainty observed in the analysis per- 



formed in . We find a simple physical reason for failure of stability in this case: the behaviour of 
quantities versus the scheme parameters of sum rules (the Borel variable or the number of spectral 
density moment) is determined by forming the doubly heavy diquark inside the baryon and, hence, 
the difference between the masses of baryon and diquark. This difference of masses determines the 
basic characteristics of correlators if we do not take into account the corrections connected to the 
nonperturbative interaction of doubly heavy diquark with the light quark in the baryon. In the 
NRQCD sum rules, the introduction of such the interaction is related with the nonperturbative 
condensates caused by operators of higher dimensions. We show that the better stability and lower 
uncertainty of sum rules can be achieved by taking into account the product of quark and gluon 
condensates in addition to quark, gluon and mixed condensates. Moreover, we accurately introduce 
the coulomb a^/f -corrections inside the heavy diquark. These corrections enforce the relative con- 
tribution of perturbative part of sum rules in comparison with the contribution of condensates into 
the correlators under consideration. 

Further, we comparatively analize the sum rules for the doubly heavy baryons with the both 
strange and light massless quarks. 

In Section 2.1 we define the currents and represent the spectral densities in the NRQCD sum 
rules for various operators included into the consideration. Section 2.2 is devoted to numerical 
estimates. We calculate the masses of ground states, which are in a good agreement with the values 
obtained in potential models. Finally, we briefly summarize the results. 
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2.1. Sum rules for doubly heavy baryons 
2.1.1. Baryonic currents 

The currents for the baryons with two heavy quarks S^^, S^^ and SJ,*^, where o denotes the electric 
charge of baryon depending on the flavor of light quark, corespond to the quantum numbers of 
spin and parity = 1+ and = O"*" for the systems of heavy diquark with the symmetric and 
anti-symmteric structure of flavor matrix, correspondingly (if the identical heavy quarks form the 
diquark then the scalar ground state = 0+ is forbidden). Adding the light quark to the system 
of heavy quarks gives = | for the S^* baryons and the couple of degenerate states = \ 
and = I"*" for the baryons S*^, S^^, and H*^, S^*, S^^. Generally, the structure of baryonic 
currents with two heavy quarks is written down in the form 

J=[Q'''CTrQYq'e,,k. (2.1) 

Here T denotes the transposition, C is the charge conjugation matrix with the properties C'~fJ^C~^ = 
—7^ and C'-f'[C~^ = 75, i,j, k are color indices and r is a matrix in the flavor space. The effective 
static field of the heavy quark is denoted hj Q. To the leading order over both the relative velocity 
of heavy quarks and their inverse masses, this field contains the "large" component only in the 
hadron rest frame. 

Here, unlike the case of baryons with a single heavy quark , there is the only independent 



current component J for each of the ground state baryon currents. They equal 



QQ' 



where JS*o satisfies the spin-3/2 condition 'y„JS*o = 0. The flavor matrix r is anti-symmetric for 
SJ,* and symmetric for Sqq and Sqq. The currents written down in Eq. (|2.2|) are taken in the rest 
frame of hadrons. The corresponding expressions in a general frame moving with a velocity can 
be obtained by the substitution of 7™ — 7^ = 7^" — f)v'^. 

Similar expressions can be written down for the doubly heavy baryons with the strange quark. 

To compare with the full QCD analysis we represent the expression for the J^/o current given 
in [||] 

Jh- = {n[u'^C^,c^]b' + r2K^Cc']75&' + rs[u'^C^57mUci]rb'}e^,k, (2.3) 

so that the NRQCD structure can be obtained by the choice of ri = r2 = 1 and = and the anti- 
symmetric permutation of c and b flavors. This connection can be achieved by the nonrelativistic 
limit of full QCD spinors of heavy quarks, so that in the leading order of 1/mQ-expansion the "large" 
components of spinors contribute only. Therefore, for the ground states of douby heavy baryons 
containing the heavy quarks with the identical flavors, the leading approximation of NRQCD leads 
to the only structure of baryonic current expressed in terms of nonrelativistic spinors of heavy 
quark, since the total spin of heavy diquark is fixed by S* = 1 because of the Pauli principle. The 
corrections of the first l/mg-order can contribute with the other Lorentz structures, of course. 
However, we deal with the leading approximation of NRQCD in the present paper. For the ground 
states of doubly heavy baryons containing the heavy quarks with the different flavors, the 1/2-spin 
state of baryon can contain the mixture of diquark states with S = 1 and S* = 0, as it does in 
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full QCD. We perform the separate consideration of these two currents in NRQCD, and such the 
approach can generally be not optimal in full QCD currents. Nevertheless, as we show, in the leading 
order of NRQCD there are relations between the masses and coupling constants of baryons because 
of the spin symmetry, so that the NRQCD does not distinguish these spin states untill the spin 
dependent l/mg-corrections are taken into account. In addition, as we have already mentioned, 
the analysis in the sum rules of full QCD was done with a large uncertainty beacuse of difference 
in the evaluations of coupling constants from two correlation functions |3^, and the authors noted 
that this uncertainty became large in a region of parameters ri 2 defined above, so that this region 
of bad accuracy is placed in the vicinity of point gining the NRQCD choice in { ^■'S\) . To get rid this 
disadvantage and to find its reason, we analyze the NRQCD sum rules in details. 



2.1.2. Description of the method 

In this section we describe steps required for the evaluation of two-point correlation functions in 
the NRQCD approximation and the connection to physical characteristics of doubly heavy baryons. 
We start from the correlator of two baryonic currents with the half spin 

U{w) = i j (fxe'P%Q\TJ{x), J(0)|0) = i)Fi{w) + F2{w), (2.4) 

where w is defined by = {Ai + w)^, while = mg + mg/ + nis, itiq^q' are the heavy quark 
masses, and is the strange quark mass. Sure, the correlators for the baryonic currents with the 
light quarks instead of the strange one can be obtained, if we put — > rriu^a ~ in expressions 
below. The appropriate definitions of scalar formfactors for the 3/2-spin baryon are given by the 
following: 

ll^uiw) = I j dWP^(0|T{ J^(x), J.(0)}|0) = ~g^,[i,F,{w) + F^iw)] + . . . , (2.5) 

where we do not concern for distinct Lorentz structures. The scalar correlators F can be evaluated in 
a deep euclidean region by employing the Operator Product Expansion (OPE) for the chronological 
product of baryonic currents in the framework of NRQCD, for instance, in Eqs. 



F,,,{w) = Y,C^^'"\^)Od. (2.6) 

d 

where Od denotes the local operator with a given dimension d: Oq = 1, O3 = (gg), O4 = (^G^), 
Od>i = ■ ■ ■, and the functions Cd{w) are the corresponding Wilson coefficients of OPE. 

In this review we take into accouont the nonperturbative terms connected to the quark and 
gluon condensates, their product and mixed condensates. For the contribution of quark condensate 
operator we explore the following OPE for the correlator of two quark fields: pT| : 



(0|TK(x)g;(0)|0) = + ^ + + •••]> (2-7) 

where the value of mixed condensate is parametrized by introducing the variable ml, which is 
numerically determined as ~ 0.8 GeV^. Taking into account the nonzero mass of strange quark 
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in the framework of OPE we get the following expression for the quark condensate up to terms of 
fourth order in x ||42|| : 



^ ^ {ml -2ml) ^ x\^^^G^) - fm^m^ - m^Y 



16 288 
1 x^ f 3mo 2 
48 + 24^1^""'^ 



(2.8) 



Note that at 7^ the expansion of quark condensate (|2.8| ) gives contributions in both correlators 
in contrast with the sum rules for Sqq/ |^3|, where putting m^ = and neglecting the higher con- 



densates, the authors found the factorization of diquark correlator in F2 and full baryonic correlator 
in Fi. This factorization led to the systematic instability for the estimates in the sum rules. 

We write down the Wilson coefficients in front of unity and quark-gluon operators by making 
use of the dispersion relation over w, 

C,(^) = irEl^, (2.9) 

TT Jo UJ -W 

where p denotes the imaginary part of corresponding Wilson coefficient in the physical region of 
NRQCD. Thus, the calculation of Wilson coefficients for the operators under consideration is reduced 
to the problem on the derivation of corresponding spectral densities. 

To relate the NRQCD correlators to the real hadrons, we use the dispersion representation for 
the two-point function with the physical spectral density given by the appropriate resonance and 
continuum part. The coupling constants of baryons are defined by the following expressions: 

(0|J(x)|S(f^)^Q(p)) = tZ^^nr^^uiv,M^^n))e'^^, 
{0\J-{x)\E{QrQ^Q{p,X)) = ^ZH(^,).o^«"^(^;,MH(^))e^^^ 

where the spinor field with the four-velocity v and mass M satisfies the equation 

i}u{v,M) = u{v,M), 

and u^{v, M) denotes the transversal spinor, so that (7"^ — v'^'^)u"^{v, M) = 0. 

We suppose that the continuum densities starting from the threshold tOcont, is modelled by the 
NRQCD expressions. Then, in the sum rules equalizing the correlators in NRQCD and those of 
given by the physical states, we assume the model of continuum given by the calculated pertur- 
bative term. This model cannot be exact because of binding effects as well as the truncation of 
perturbative expansion in the given order of a^. Therefore, the integration above ujcont cannot be 
strictly cancelled, and the model introduces the implicit dependence of masses and couplings on the 
choice of value Ucont- This dependece causes an uncertainty, which is not essential in comparison 
with uncertainties following from another methodics and the variation of quark masses. 

Then we use the nonrelativistic expressions for the physical spectral functions 

P??M = ^l^r5(A-^), (2.10) 
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where we have performed the substitution — M^) — ^ "^(^ ^^^^ ^ is the binding energy 
of baryon and M = A4 + A. The nonrelativistic dispersion relation for the hadronic part of sum 
rules has the form 

Llif = J \ 2 11) 

u-w 2MA-w' ^ ■ ^ 

Further, we write down the correlators in the deep underthreshold point of w = —A4 + 1 at t — 0, 
which corresponds to the limit of — ^ 0. The approximation of hadronic part by the only bound 
state leads to the expression, which can be expanded in series of t. Thus, the sum rules result in 
the equality of coefficients at the same powers of t 

^ ' ' ' (2.12) 



vr 7o {yj^ MY 2M M'^ ' 

where pj contains the contributions given by various operators in OPE for the corresponding scalar 
functions Fj. Introducing the following notation for the n-th moment of two-point correlation 
function: 

for the baryon mass M= we have the following estimate: 

M[n] = (2.14) 
and the coupling constants are determined by the expression 

\Z[n]\' = '^MnM-^\ (2.15) 

where we see the dependence of sum rule results on the scheme parameter. Therefore, we have to 
find the region of parameters, wherein, first, the results are stable with respect to variation of n, 
and, second, the both correlation functions Fi and F2 repeats the same values of physical quantities: 
the masses and coupling constants. The problem of analysis given in the full QCD was the existance 
of significant difference between the masses and coupling constants calculated under different F. 

2.1.3. Calculating the spectral densities 

In this subsection we present analytical expression for the perturbative spectral functions in the 
NRQCD approximation. The evaluation of spectral densities involves the standard use of Cutkosky 
rules with some modifications motivated by NRQCD. We use the rule under which the jump of 



two-point function under study is calculated by the following substitutions for the propogators of 
heavy and fight quarks, respectively: 

1 p2 

heavy quark: ^ 2m ■ S{po - {m + —)), 

Po-{m + ^) 2m 

light quark : — 2Tii ■ 5{p'^ — m^). 

We derive the spin symmetry relations for all the spectral densities due to the fact that in the 
leading order of the heavy quark effective theory the spins of heavy quarks are decoupled, so 

Pi,n(H)«^ = 3Pi,f,(H)'^-^, = 3pi,n(H)-,, (2.16) 
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p2,mrQQ = ^P2,n(s)'^Q, = 3P2,n(H)-,, (2.17) 
and we have the following relation for the baryon couplings in NRQCD: 

= 3|Zf7'(s)|^ = 3|Zn(H).|^. (2.18) 

For the perturbative spectral densities pi^uil^) and P2,h{^^) in front of unit operator for Fi and 
F2, respectively, we explore the smallness of strange quark current mass with respect to the masses 
of heavy quarks and make the expansion in series of m^. So, we get the following expressions: 

Pin'- , (^) = ip^nip; 3^^T \ ^(^l.oH + ^sVl,l{^) + "^sm,2(^)), (2.19) 



so that mQQi = mQmqi / {mq + mqi) is the reduced diquark mass,, M.diq = f^Q + "^q', and the 
coefficients of spectral densities 77 are presented in Appendix The first term of this expansion 
reproduces the result obtained in for the zero mass of light quark. For the strange baryons the 
perturbative density Pgf^'* prportional to m^, and it is not equal to zero 



QQ' 



P2,n''> , (^) = 3(\^ — —^(^2,0 + msV2,i + "^3^2,2)- (2.20) 

' QQ' lUOTT'^^Alrfjg + CUj^ 

In the leading order of perturbative NRQCD the correlators F2 are equal to zero for the massless 
light quark. This fact is caused by the absense of interaction between the light quark and the heavy 
diquark in this order, and, therefore, there is no massive term in this correlator. 

The coulomb interaction inside the diquark can be taken into account by the introduction of 
Sommerfeld factor C for the spectral density of diquark before the integration over the invariant 
mass of diquark in order to get the baryonic spectral densities, so that 

P diquark P diquark ^ (2'21) 

whereas 



27ra, 



3v, 



QQ' 



1 — exp 



27ias 



-1 



(2.22) 



where we have taken into account the anti-triplet color structure of diquark, and vqq> denotes the 
relative velocity of heavy quarks inside the diquark: 



V Q- - mQ>y 

where is the square of heavy diquark four- momentum. In NRQCD we take the limit of low 
velocities, so that denoting the diquark invariant mass squared by = {M.diq + e)^, we find 

— , Ur 



at e <^ rriQQ'. The modified spectral densities are equal to 



mlQ,asUj{2Mdiq + uj) ^ c , 2 c 



-"1 = a2(u fao + ^sVi,i + mivi,2)- (2.24) 



^In what follows the coefficients of spectral densities, which are not given explicitly, are also given in Appendix I. 
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The leading term gives the result for the zero mass of light quark . For pz we have 

= (KA 1 ^2 ("12,0 + msvt + ^s^s'^a)- (2.25) 

The truncated expansion in the mass of light quark leads to a small deviation about 0.5% from 
the exact integral representation, so that this approximation is quite justified with the initial three 
terms in the explicit analytic form. 

Further, the spectral functions, connected to the condensates of light quarks and gluons, can be 
derived in the analogous way. For the moments of coefficients in front of quark condensate terms 
we have the following expressions: 

n\ nl 

Mfhn) = VoM''"'{n) - ^,'^^' V2M''"^{n + 2)+ ^'^ ^ V^M''"' [n + 4), (2.26) 

nl nl 



where we have used the coefficients of expansion in x, Vi, in Eq. (|2.8|) , while the n-th moment of 
two-point correlation function for the diquark is denoted by A1*^(n), and it is derived under the 
integration of spectral density 

Pci^q = — (2.27) 

TT 

which should be multiplied by the Sommerfeld factor C, where the variable e is substituted by u, 
since in this case there is no integration over the invariant mass of diquark. The modified density 
is equal to 

P% = (2.28) 

and it does not depend on uj. 

It is interesting to stress that in NRQCD the light quark condensate contributes to the F2 
correlators, only. This fact has a simple physical explanation: to the leading order the light quark 
operator can be factorized in the expression for the correlator of baryonic currents. Indeed, we can 
write down for the condensate contribution 

(0|T{ J(x), J(0)}|0) ^ (0|r{g«(a;)gt(0)|0) ■ 1 ■ (0|T{ Ji(x), Ji(0)}|0) + . . . , 

where Jj(a;) denotes the appropriate diquark current with the color index j, as it is defined by 
the baryon structure in Eqs. (|2.2| ). So, we see that the restriction by the first term independent 
of X in the expansion for the quark correlator in ( |2.7D results in the independent contribution of 
diquark correlator to the baryonic one. Then, since the diquark correlator is isolated in F2 from 
the baryonic formfactor Fi, the NRQCD sum rules lead to the evaluation of diquark masses and 
couplings from F2, and estimation of baryon masses and couplings from Fi. These masses and 
couplings are different. The positive point is the possibility to calculate the binding energy for the 
doubly heavy baryons A = Ms — M.diq- The disadvantage is the instability of NRQCD sum rules at 
this stage, since the various formfactors or correlators lead to the different results. In sum rules of 
full QCD various choices of parameters in the definitions of baryonic currents result in an admixture 
of pure diquark correlator in various formfactors, so that the estimations acquire huge uncertainties. 
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Say, the characteristic ambiguity in the evaluation of baryon mass in full QCD is about 300 MeV, 
i.e. the value close to the expected estimate of A. The analysis in the framework of NRQCD 
makes this result to be not unexpectable. Moreover, it is quite evident that the introduction of 
interactions between the light quark and the doubly heavy diquark destroys the factorization of 
diquark correlator. Indeed, we see that due to the higher terms in expansion ( |2.7| ), the diquark 
factorization is explicitly broken, which has to result in the convergency of estimates obtained from 
Fi and F2. Below we show numerically that this fact is valid. Technically, we point out that the 
contribution to the moments of spectral density, determined by the light quark condensate including 
the mixed condensate and the product of quark and gluon condensates, can be calculated after the 
exploration of (p7| ), so that 

For the corrections determined by the gluon condensate and written down for the operator O4 = 
i^G"^), we have 

21 ■ 2^oV2TTmlml,{Mdiq + cu)^ 
For the nonzero mass of light quark the density p^^ (u) is proportional to and equal to 



pf (-) = : 1 ; 7 r "'" + -.c + (2-30) 



so that 



Pfi^) = , 2^J\_V + ^svt)^ (2.31) 

3 ■ 2V27rm^m^,(A<dig + u) 

V% = -{9M,., + u), vS = ^-^^- (2.32) 

For the product of condensates (qq) {^G"^) , wherein the gluon fields are connected to the heavy 
quarks in contrast to the light quark, we have computed the contribution to the two-point correlation 
function itself. It has the following form: 

F^iG'(^) ^ ^Q^'i^Q + K' + ^^'^Q^Q') .2 33) 

and F^'^'^ (cu) = 0, so that we restrict the consideration by the operators with the dimension not 
greater than 7, while the nonzero contribution to Fi appears in the fifth order of expansion ( p.8|) . The 
obtained result is presented in the form, which allows the analytic continuation over u = —M. + w. 

Thus, we provide the NRQCD sum rules, where we take into account the perturbative terms 
and the vacuum expectations of quark-gluon operators up to the contributions by the light quark 
condensate, gluon condensate, their product and the mixed condensate. Note, that the product of 
condensates is essential for the doubly heavy baryons, and we present the full NRQCD expression 
for this term, including the interaction of nonperturbative gluons with both the light and heavy 
quarks. The correct introduction of coulomb-like interactions is done for the perturbative spectral 
densities of heavy diquark, which is important for the nonrelativistic heavy quarks. Finally, we find 
the spin-symmetry relation for the baryon couplings in NRQCD 

= 3|Z=/P = 3|Z=.|l 
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2.1.4. Anomalous dimensions for the baryonic currents 

To connect the NRQCD sum rules to the quantities in full QCD we have to take into account the 
anomalous dimensions of effective baryonic currents with the nonrelativistic quarks. They determine 
the factors, which have to multiply the NRQCD correlators to obtain the values in full QCD. Indeed, 
to the leading order of NRQCD we have the relation 

where the coefficient /Cj(as, /isoft, /Uhard) depends on the normalization scale //soft and obeys the 
matching condition at the starting point of /Xhard = -Mdiq- The anomalous dimensions of NRQCD 
currents are independent of the diquark spin structure in the leading order. They are equal to 



d In C J {as, jJ') / a 



^' ' m=l 



(^-2Cs(3a-3) + 3Cj.(a-2)), (2.34) 
^(2) = i(-48(-2 + 6C(2))C2 +C^((104- 240C(2))CB-101Ci.) 



-GACBUfTF + CFi-9CF + 52nfTF)), 

where Cf = {N^ - l)/2Nc, Ca = N^, Cb = (A^c + l)/2N^, and Tp = 1/2 for = 3, n/ being 
the number of light quarks. In Eq. (|2.34|) we give the one- loop result with the arbitrary gauge 
parameter a, and the two-loop anomalous dimension is represented in the Feynman gauge a = 1. 
So, numerically at = 3 and a = 1 we find 

7^^) = -4, 7(2) ^ -188.24. (2.35) 

In the leading logarithmic approximation and to the one-loop accuracy, the coefficient JCj is given 
by the expression 

r- f \ ( "^(/^hard) ^ 

/Cj(Q;s, /isoft,Athard) = — 7 T , (2.36) 

V «s(/isoft) / 

where /3o = lliVc/3 — 2n//3 = 9. To evaluate the two- loop expression for Cj we have to know sub- 
leading corrections in the first order in addition to the anomalous dimensions. These corrections 
are not available yet, so we restrict ourselves by the one- loop accuracy. 

Further, we have to determine the normalization point for the NRQCD estimates /isoft- We put 
it to the average momentum transfer inside the doubly heavy diquark, so that nl^^^ = 2mQQ/Tdiq, 
where Tdiq denotes the kinetic energy in the system of two heavy quarks, which is phenomenologically 
independent of the quark fiavors and approximately equal to 0.2 GeV. Then, the coefficients JCj are 
equal to 

/Cf,(H).. ^ 1.95, /Cn(H),. ^ 1.52, /Co(H),, ^ 1.30, (2.37) 
with the characteristic uncertainty about 10% because of the variation of initial and final points 

/^hard,soft. 

Finally, we emphasize that the values of /Cj do not change the estimates of baryon masses 
calculated in the sum rules of NRQCD. However, they are essential in the evaluation of baryon 
couplings, which acquire these multiplicative factors. 
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2.2. Numerical estimates 



Evaluating the two-point sum rules, we explore the scheme of moments. We point out the well- 
known fact that an essential part of uncertainties is caused by the variation of heavy quark masses. 
In the analysis we chose the following region of mass values: 

mb = 4.6 - 4.7 GeV, = 1.35 - 1.40 GeV, (2.38) 

which is ordinary used in the sum rule estimates for the heavy quarkonia. Next critical point is 
the value of QCD coupling constant determining the coulomb-like interactions inside the doubly 
heavy diquark. Indeed, it stands linearly in front of the perturbative functions of diquark contri- 
butions. Thus, the introduction of a^/f-corrections is essential for both the baryon couplings and 
the relative contributions of perturbative terms and condensates to the baryon masses. To decrease 
the uncertainty we impose the same approach to the heavy quarkonia, where it is well justified, 
and then, we extract the characteristic values for the heavy-heavy systems from the comparison of 
calculations with the current data on the leptonic constants of heavy quarkonia, which are known 
experimentally for cc and bb or evaluated in various approaches for be. So, our calculations give the 
following couplings of coulomb interactions 

as{bb) = 0.37, a,(c6) = 0.45, a,(cc) = 0.60. (2.39) 

Since the squared size of diquark is two times larger than that of heavy quarkonium composed of 
the same heavy quarks (see the dependence of average square of relative momentum on the kinetic 
energy of heavy quarks), the effective coulumb constants have to be rescaled according to the 
equation of evolution in QCD. Since we use the one- loop approximation, we explore the evolution 
equation 

'^^^^ l-^«.(QQ01n2- 

So, 

a,(bb) = 0.45, a,(6c) = 0.58, a,(cc) = 0.85. (2.40) 

As for the dependence of results on inputs for the quark masses, we have to remark that so called 
pole masses are not well defined due to infrared problems, usually mentioned as the renormalon 



ambiguity Thus, it is important to fix the definition of mass 

To the given order in for the NRQCD sum rules, we use the leading quark loop approximation 
with account for the coulomb exchange between the heavy quarks. At this stage the heavy quark 
masses and coulomb coupling constants are strictly fixed by the data on the charmonium and bot- 
tomonium leptonic constants and masses as described by the QCD sum rules to the same accuracy. 
The stability or convergency of sum rule method applied to these heavy quarkoniaQ results in the 
following masses of quarks: 

= 1.40 ± 0.03 GeV, m^, = 4.60 ± 0.02 GeV, 

which well agree with the values of heavy quark masses defined as free of infrared contributions: 
the potential subtracted mass m^^ = 4.60 ± 0.11 GeV and the kinetic mass m^™ = 4.56 ± 0.06 GeV 
both obtained in the QCD sum rules for the bottomonium within the two- loop accuracy |]46|,^ . 



^We have required that the ratio of initial moments for the spectral densities calculated over the data and in QCD 
sum rules was stable. 
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The corresponding IS-mass defined in has a shghtly larger value. We think that in the leading 
order over the PS and kinetic masses above determine the threshold of quark contribution and 
can be taken as the appropriately defined heavy quark masses in the calculations of characteristics 
for the doubly heavy baryons. The mass values are dependent of normalization point, which was 
chosen in the range of 1 — 2 GeV. Nevertheless, we slightly enlarge the region of mass variation. 

The QCD sum rules for the bottomonium and charmonium fix the values of coulomb couplings, 
too, since the momentum stability yields the heavy quark masses, while the leptonic constants 
linearly determine the corresponding values of as shown above (see Fig. |2.1| ). Note that the 
dependence of coulomb coupling constant on the quark contents of quarkonia well agree with the 
renormalization group evolution with the change of size for the system composed of two heavy 
quarks. The uncertainty of further estimates on the supposed values of coulomb couplings is about 
5%. 
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Figure 2.1. The leptonic constant of T in the two-point sum rules formulated in the scheme of spectral density 
moments. The dashed line represents the result at = 4.63 GeV, the solid curve gives the estimate at 
mfe = 4.59 GeV. 



As we have already mentioned, in the coupling constants of baryons the uncertainty connected 
to the coefficients for the matching the NRQCD with full QCD, is about 10%. 

The dependence of estimates on the value of thershold for the continuum contribution is not so 
significant as on the quark masses. We fix the region of ujcont as 



cont 



1.3 - 1.4 GeV. 

For the condensates of quarks and gluons the following regions are under consideration: 



(2.41) 



ml 



-(250 - 270 MeV)3, 
0.75 - 0.85 GeV^ 
(1.5 - 2) ■ 10-2 GeV^ 



(2.42) 



43 



The main source of uncertainties in the ratios of the baryonic couphngs is the ratio of the condensates 
with the strange and hght quarks. We use (ss)/ (qq) = 0.8 ± 0.2 that corresponds to the variations 
of sum {rriu + md)[l GeV] = 12 ^ 14 MeV |9|. 

We suppose the strange quark mass equal to = 150 ± 30 MeV, that is the wide-accepted 
estimate consistent with both the sum rules and the quark current algebra, wherein the current 
mass of quark operates. 

So, we have described the set of parameters entering the scheme of calculations. 

Fig. |2]^ represents the calculated difference of masses extracted from the Fi and F2 correlators^ 
for the baryon (we have not shown the similar figures for the Sec and Sf,;, baryons, since they 
qualitatively and quantitatively repeat the picture clearly given by Fig. |2.2|). 




Figure 2.2. The difference between the H^c-baryon masses calculated in the NRQCD sum rules for the form- 
factors Fi and F2 in the scheme of moments for the spectral densities. 

We certainly see that at low numbers of moments for the spectral densities, the baryon-diquark 
mass difference can be evaluated as 

A = 0.40 ± 0.03 GeV, (2.43) 

which is quite a reasonable value, being in a good agreement with the estimates in the heavy-light 
mesons. In the region of mass difference stability we can fix the number of moment for the spectral 
density, say, n = 27 ± 1 for Ef,c, and calculate the corresponding masses of baryons, which are equal 
to 

Ms,, = 3.47 ± 0.05 GeV, Ms,, = 6.80 ± 0.05 GeV, Me,, = 10.07 ± 0.09 GeV, (2.44) 

^In these figures we have fixed the value of gluon condensate (^G^) = 1.7 • 10"'^ GeV* and arranged mg in the 
above region to reach zero differences between the masses, though the variation of parameters leads to errors in the 
estimates quoted below. 
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where we do not take into account the spin-dependent sphtting caused by the a^-corrections to 
the heavy-hght interactions, which are not available yet. The uncertainties in the mass values 
are basically given by the variation of heavy quark masses. The convergency of NRQCD sum 
rules allows one to improve the accuracy of estimates in comparison with the previous analysis in 
full QCD The obtained values are in agreement with the calculations in the framework of 



nonrelativistic potential models (see Chapter 1). 



GeV 



6.9 



6.8 



6.7 



6.6 




5 10 15 20 25 30 35 



n 



Figure 2.3. The Ebc (lower curve) and Qbc (upper curve) masses obtained in the NRQCD sum rules by averaging 



the results from two correlators Fi 2 



In the two-point sum rules for the mass of Qf,c (the conclusions for other doubly heavy baryons 
are similar) we can observe the stability of estimates with respect to changing the moment numbers 
in both correlators Fi and F2. We suppose that this fact is connected to the destroying of mentioned 
factorization for the correlators of diquark and baryon in the perturbative limit in contrast to the 
case of S-baryons. The stability regions for Fi and F2 do not coincide because the contributions 
of higher dimension operators become valuable at the different numbers of moments. However, the 
quantity 

^(MiM+MsM) 

has the larger stability region, and we explore this fact to determine the Q baryons masses as well as 
that of S baryons (see Fig. p.3|) . Thus, in the present review we consider two criteria for the stability 
of baryon mases: the first is based on the study of mass difference obtained from two correlators 
F12, the second investigates the half sum of masses extracted from two correlators. The second way 
is especially reliable for the doubly heavy baryons with the strangeness, because the both correlators 
have the stability regions at various numbers of moments. In this way, the difference between the 
masses in two correlators at the stable points determines the accuracy of estimates in the framework 
of NRQCD sum rules. 
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Figure 2.4. The mass difference AM = Mq^^ — M^^^ obtained from the results shown in Fig. 



The second method results in the following values of baryon masses: 



Mn 



3.65 ±0.05 
6.89 ±0.05 
10.09 ±0.05 



GeV, 
GeV, 
GeV, 



3.55 ±0.06 GeV, 
6.79 ±0.06 GeV, 
10.00 ±0.06 GeV, 



(2.45) 



We can see that the estimates of Hqq/ masses in both methods of ( p.45|) and ( p.44| ) well agree with 
each other. 

Then, we investigate the difference between the masses of doubly heayy baryons with strangeness 
and wothout it: 1/2((Mi_q ± M2^n) — (Mi = ± M2^e.)), shown in Fig. |2.4 In our scheme of baryon 
masses determination this quantity has the meaning of average difference between the masses, for 
which we observe a wide interval of stability indicating a good systematic accuracy of estimate. We 
have obtained 



AM 



Mr, 



M= 



100 ± 30 MeV. 



Figs. |2]^, ^]6| show the dependence of baryon couplings calculated in the moment scheme of 
NRQCD sum rules for the doubly heavy baryons with the strangeness and without it, respectively. 
Numerically, we find 



(2.46) 



\Zn P = 


(10.0 ± 1.2) 




GeV^ 


|Z= P = 

1 — 'CC 1 


(7.2 ±0.8) 


10-^ 


GeV 




(15.6 ± 1.6) 




GeV^ 


I^hJ^ = 


(11.6 ± 1.0) 




GeV 


l^nttl = 


(6.0 ±0.8) 


10-2 


GeV^ 


I^hJ^ = 


(4.2 ±0.6) 


10-2 


GeV 



In Fig. we present the sum rules results for the ratio of baryonic constants \Zq^J 
We have also found 



1.3 ±0.2. 
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Figure 2.5. The couplings jZ^^J^^p calculated in the NRQCD sum rules for the formfactors Fi and F2 (solid 
and dashed lines, correspondingly) in the scheme of moments for the spectral densities. 

The uncertainty of this resuh as was mentioned above is mainly connected with the pourly known 
ratio of {ss)/{qq) = 0.8 ± 0.2. 

We can see in figures that the region of stabihty for the baryonic constants coincides with the 
region of stabihty for the average mass described above. 

For the sake of comparison, we derive the relation between the baryon coupling and the wave 
function of doubly heavy baryon evaluated in the framework of potential model, where we have 
used the approximation of quark-diquark factorization. So, we find 



7PM I 



2v^|^rf(0)-^^,,(0)|, (2.47) 



where \l'd(0) and \l'/,s(0) denote the wave functions at the origin for the doubly heavy diquark and 
hght (strange) quark-diquark systems, respectively. In the approximation used, the values of ^^(0) 
were calculated in the potential by Buchmiiller and Tye |T^, so that 



4^1^,(0)1 = 0.53 GeV^/^^ V47r [^^(O)! = 0.64 GeV3/^ 
4^|^^^(0)| = 0.53 GeV^/^ v^|^fee(0)| = 0.73 GeV^/^ 
4^|$b^(0)| = 1.35GeV^/l 



In the static limit of potential models, these parameters result in the estimates 

\Z^^\' = 8.8 ■ 10-=^ GeV^ = 6.0 • 10-^ GeV^ 

PM|2 _ i« in-2n„^r6 i'7PM|2 _ 11 in-2n„AA6 



|Z™f = 1.6■10-'GeV^ \Z^y = l.l•10-'GeV^ (2.48) 
\zi^\' = 5.6- 10-2 GeV^ \zi^\^ = 3.9 ■ 10"^ GeV^ 



The estimates in the potential model ( p.48|) are close to the values obtained in the sum rules of 

NRQCD (|2.46| ). We also see that the SU(3)-flavor splitting for the baryonic constants 

is determined by the ratio |\E's(0)p/|\I';(0)p = 1.45 which is in agreement with the sum rules result. 
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Figure 2.6. The couplings calculated in the NRQCD sum rules for the formfactors Fi and F2 (solid 

and dashed lines, correspondingly) in the scheme of moments for the spectral densities. 
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Figure 2.7. The ratio , ^^'"',■2 calculated in the NRQCD sum rules in the scheme of moments for the spectral 



densities at {ss)/{qq) = 0.8. 
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The values obtained in the NRQCD sum rules have to be multiphed by the Wilson coefficients 
coming from the expansion of full QCD operators in terms of NRQCD fields, as they have been 
estimated by use of corresponding anomalous dimensions. This procedure results in the estimates 





= (38 ± 5) • 


• 10" 


"3 GeV^ 


1 ^cc 1 


= (27 ±3) ■ 


■ 10" 


"3 GeV^ 




= (36 ± 4) • 


• 10" 


"3 GeV^ 


\Z^b 1' 

1 —'be 1 


= (27 ±3) ■ 


■ 10" 


"3 GeV^ 




= (10±1). 


• 10- 


-2 GeV^ 


\Z^bf 


= (70 ± 8) ■ 


• 10" 


GeV^ 



Thus, we have got the reliable estimates for the masses and coupling constants of baryons 
containing two heavy quarks in the framework of NRQCD sum rules. 



2.3. Discussion 

We have considered the NRQCD sum rules for the two-point correlators of baryonic currents with 
two heavy quarks. The nonrelativistic approximation for the heavy quark fields allows us to fix 
the structure of baryonic currents and to take into account the coulomb-like interactions inside the 
doubly heavy diquark. Moreover, we have introduced into the consideration the operators of higher 
dimensions, which are responsible for the quark-gluon condensates in order to reach the convcrgcncy 
of the sum rule method for two scalar correlation functions. To the leading approximation, including 
the pcrturbative term and the contributions of quark and gluon condensates, the correlators of three- 
quarks state and the doubly heavy diquark arc factorized in separate functions, so that the sum 
rules result in the different values of masses and couplings. This fact indicates the divergency of 
approach unless the product of quark and gluon condensates and the mixed condensate are taken 
into account. Then, the interaction of two heavy quarks and light quark destroy the factorization, 
which allows one to get meaningful estimates of masses and couplings. Moreover, we have also 
calculated the binding energy of doubly heavy diquark, which is in a good agreement with the 
estimates in the framework of potential models. In the doubly heavy baryons with the strangeness 
the factorization of diquark correlator is already destroyed in the approximation of quark loop in 
the pcrturbative QCD, so that the better convergency of sum rule approach takes place. So, the 
both correlators have the intervals of stability under the variation of number for the moment of 
spectral density. We have calculated the splitting of masses between the strange baryons Qqq' and 
baryons with the massless light quark Sqq/. We have also estimates the ratio of baryonic coupling 
constants I-^q^^^ P/I^s^^, P for the baryons under consideration. 

Thus, the NRQCD sum rules allow us to improve the analysis of masses and couplings for the 
doubly heavy baryons and to obtain reliable results. 
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Chapter 3. Production processes 

For the production of baryons with two heavy quarks Sqq/ a small ratio of A/mg and, hence, a 
small value of quark-gluon coupling constant in QCD as ~ l/ln{mQ/A) ^ 1 allow us not only to 
consider the production of two pairs of heavy quarks QQ and Q'Q' in the framework of perturbative 
QCD, but also to factorize the contributions caused by the perturbative production of heavy quarks 
and the nonperturbative forming of heavy diquark by the quarks cnctering the Sggz-baryons. So, 
in order to calculate the cross sections for the production of S-wave states at the Z-boson peak 
we have to evaluate the matrix elements for the associated production of bb and cc pairs with the 
anti-triplet color state of be pair having got a definite sum of quark spins {S = 0, 1), while the 
quarks are assigned to the same velocity equal to the velocity of diquark composed by these quarks. 
Then we have to multiply these matrix elements by a nonperturbative factor determined by some 
spectroscopic characteristics of bound state, i.e. by both the wave function of diquark, that gives 
the probability to find the quarks at short distances between them inside the bound state, and 
the quark masses. This approach is justified due to the following: the characteristic virtualities 
of heavy quarks inside the heavy diquark are much less than their masses, since the heavy quarks 
nonrelativistically move in the bound state, while the quark virtualities are about their masses in 
the hard production. Therefore, considering the production of 3^^, we can suppose that the quarks 
b and c are on the mass shell in the diquark, and they are at rest with respect to each other. Thus, 
after the isolation of nonperturbative factor the analysis of Sfec-baryon production is reduced to the 
consideration of matrix elements calculated in the perturbative QCD, if we suppose that the total 
and differential cross sections of baryon repeat the corresponding quantities for the doubly heavy 
diquark. 

Note, first of all, that in the electromagnetic and strong interactions of colliding particles the 
associated production of two heavy quark pairs necessary for the hadronization into Hqq/, can be 
done provided the leading order of perturbative QCD for such processes has the additional factor 
of in comparison with the leading order for the production of single pair of heavy quarks QQ, 
so that (t{Eqq')/(t{QQ) ~ a^l"^ {0)1"^ / rriQ, . This suppression causes a small relative yield of Sqq/ in 
comparison with the inclusive production of heavy mesons Mq. 

So, we have to analyze the leading approximation of perturbative QCD for the production of 
Sqq/, that allows us to get some analytic expressions for the cross sections of Sqq/, particularly, 
for the fragmentation functions of both the heavy quark into the heavy diquark and the diquark 
into the baryon in the scaling limit of M^/s 0. Thus, the production of Sqq/ in the regime of 
fragmentation can be reliably described in terms of analytic expressions, that opens new possibilities 
in the study of QCD dynamics, which is significant in the complete picture of heavy quark physics. 

The mechanism for the production of baryons with two heavy quarks in the hadronic collisions 
involves the consideration of complete set of diagrams in the fourth order of perturbative QCD 
because the fragmentation regime does not dominate in the total cross section determined by other 
nonfragmentational contributions rapidly decreasing at large transverse momenta. We investigate a 
role of these higher twists over the transverse momentum of doubly heavy baryon in the associated 
hadronic production. We numerically determine the limits for the consistent use of factorization 
regime in the hard production of heavy quarks with the consequent fragmentation. 
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3.1. Production of doubly heavy baryons in e^e -annihilation 

The analysis of mechanisms for the production of hadrons with two heavy quarks shows that the 
expected production yield of such hadrons with respect to the number of hadrons with a single heavy 
quark is of the order of 10^'^^^'^-'. For example, at the Z°-boson pole the number of events with 
heavy quarks is about 10^, consequently, the number of hadrons with two heavy quarks is expected 
to be ~ 100 — 1000. Taking into account specific decay modes of hadrons with two heavy quarks 
one could expect the detection of several events with such hadrons, which makes their observation 
problematic at LEP. 

In this section we consider the doubly charmed baryon production (Scc^) under the conditions 
of a 5-factory with high luminosity L = 10^^ cm~^s~^, where the number of Sc*'' is two orders of 
magnitude greater than the yield at the Z°-boson pole. 



3.1.1. Fragmentation mechanism 

In the production cross sections for sic\ and were evaluated in the region of heavy quark 



fragmentation at high energies. These estimations are based on the exact analytical calculations for 
the heavy quarkonium production in the QCD perturbation theory in the limit of small M^/s ratio 
and nonrelativistic potential model. In the cc-diquark momentum spectrum was considered to 
be equal to that of heavy vector quarkonium {cc)f\ 

Dc^Uz) = f al{^ml)F{z), (3.1) 



c 



where 



z(l - zV' 

F{z) = (16 - 32z + 72^2 _ 32z^ + 5z'^) 



(2 

and -Rcc(O) is a radial wave function of the bound diquark at the origin. 

Let us note that identical quarks cc in the color anti-triplet state can have the symmetrical spin 
wave function in the S'-wave, i.e. they must be in the total spin S = 1 state. The normalization of 
the fragmentation function Dc^cc{z) is determined by the model dependent value of Rcc{0). In |SC 



a rather rough approximation with the pure coulomb potential in the system of heavy quarks was 
used. This factor gives a huge uncertainty^ in the estimation of sic^ Moreover, expression (3.1) 
obtained in the scaling limit M^/s 0, is not justified for estimating the si*^ production at the 
-B-factory, where the M'^/s ratio is not small. 

We propose another method to estimate the production of hadrons containing two heavy quarks, 
on the basis of quark-hadron duality. 



3.1.2. Calculations under the quark-hadron duality 

A production cross section of the i?c-nieson S'-wave states at the Z°-boson pole calculated in the 
fragmentation model is in a good agreement with the cross section estimations for the production 
of quark pair (be) in the color singlet state with small invariant masses 

mb + mc< M{bc) < Mth = Mb + Md + AM, (3.2) 
^There is a wrong additional factor of 2 in jscl ]. 

^°The calculation of diquark wave function in the model with the Martin potential accounting for the factor of 1/2 
due to the anti-triplet color state of quarks enforces the corresponding factor by an order of magnitude. 
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Table 3.1. The production cross sections of r]^ and -0 mesons in e+e annihilation at the B-factory. 



meson 


cr, pb 


VcilS) 


0.025 


Vc{2S) 


0.003 




0.055 




0.010 



where AM = 0.5-1 GeV. 

In the same range of duality ( p.2| ) the 6c-diquark production cross section is approximately equal 
to that of 6c-pair. Selecting the color anti-triplet state (6c) and multiplying by the factor of 2/3, we 
obtain the estimate for the the s[*^-baryon production cross section cl^^^J )/a{hh) ~ 6 ■ 10-^ i.e. 6 
times greater than the estimate made in |Q for the IS'-state production. This difference is caused. 



first, by taking into account the contribution by higher excitations of diquark in the framework of 
quark-hadron duality and, second, by the strong suppression due to supposed low value of i?6c(0), 
evidently underestimated in the pure coulomb approximation. 

Let us consider the Scc''-baryon production at the energy of i?-factory {^/s = 10.58 GeV). 
Remember that expression ( |3.1| ) may not be used at the given energy because the power corrections 
over M^/ s are substantial. The method of calculations in the leading order of QCD perturbation 
theory was described in pT|-|5^. 

In the method of quark-hadron duality the cross section for the associated production of quarko- 
nium bound states can be estimated by using the formula 



a(e+e- - inLic-c)j)c-c) = dM^, da{eU- -. {ccU^,,^,cc) ^ 

—J Jm. (^^^^c-c 

where Mi = 2mc is the kinematical threshold of cc-pair production, Mth = ^M^ + AM, AM = 
0.5 — 1 GeV. In Table ^.l| the results for the numerical calculations of the QCD perturbation theory 
diagrams are presented for the production of the bound IS"- and 2S'-levels of charmonium at the 
energy y/s = 10.58 GeV and = 0.2. The values of the radial wave functions at the origin 
Rns{0) have been determined from the experimental data on the lepton decay widths of charmonia 
ip (nS) p3| . As we can conclude from Table p.l| , below the threshold for the decay into the pair of 
DD mesons the sum over the S-wave states of charmonia is equal to 

aC^ri^ip) = 0.093 pb. (3.4) 

Note that the ratio of the vector and pseudoscalar state yields at the energy under consideration is 
equal to ujy/ujp ~ 2.2 in contrast to the value a;y/a;p ~ 1 obtained in the fragmentation mechanism 
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Our estimations of the integral in the r.h.s. of expression (p.3|) give 



acc{^M = 0.5 GeV) = 0.093 pb, (3.5) 
acg(AM=lGeV) = 0.110 pb, (3.6) 



where we set = 1.4 GeV. 
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From (|3.4|) and ( p.5|) , ( ^I6D we see that the relation of quark-hadron duahty (|3.3|) is well satisfied 
for the bound states of (cc). 

Further, the calculations show that the invariant mass distributions for the cc and cc pairs 
coincide with each other in the region of small invariant masses. Hence, in the same duality region, 
the estimates for the production cross sections of cc pair and cc-diquark are approximately equal to 
each other, so that comparing with (|3.5| ) and (|3.6|) we have 



acc{AM = 0.5 GeV) = 0.086 pb, (3.7) 

acc(AM=lGeV) = 0.115 pb. (3.8) 

Selecting the anti-triplet color state, we can obtain the summed total cross section for the production 
of Sec'' baryons 

a(sW) = (70 ± 10) ■ 10"^ pb, (3.9) 
so that the yield fraction of doubly charmed baryons approximately equals 

a(SW)Mcc) = 7- 10-^ (3.10) 



For the luminosity equal to L = 10^'^cm ^s ^ the number of events with the production of si*'' is 
equal to N{Ei*J) = 7 ■ 10^ per year, so it is by two orders of magnitude greater than the yield of 
Ei*^ at LEP. 



In ref. [£5| the distribution over the cc-diquark momentum is presented for the conditions of 
anti-symmetric collider KEK. 

Thus, in this section we have presented the calculations of the doubly charmed Hcc baryon 
production on the basis of the quark-hadron duality in the leading order of QCD perturbation 
theory. We have evaluated the E^ production cross section at the energy of B factory, where the 
fragmentation model ||5^ does not work. 

The main theoretical uncertainty in the estimations for the production cross section of the 
double charmed baryons is related with the description of the process for the heavy cc-diquark 
hadronization. First of all, a considerable fraction of the diquarks (1/3) is produced in the color 
sixtet state and can transmit into both the exotic four quark states (ccqq) and the DD-meson pair. 
As in ref. we assume that color anti-triplet state hadronizes into the Ei*J baryon with 100% 

4 



probability. Thus, at the S-factory one could expect 10 events per year with the production of 
Sec'* at the luminosity L = 10^^ 

3.1.3. Exclusive production of diquark pairs 

Close to the threshold for the production of doubly heavy baryons in e^e~ annihilation a significant 
contribution could be given by the pair production. In order to estimate a yield of such kind events 
the calculation of cross sections was done in [^] for the exclusive production of doubly heavy 
diquark pairs. Supposing the 100% probability for the diquark fragmentation into the baryons we 
can put the pair yield equal to the yield of baryons. Authors of considered both the axial- vector 
states and the scalar ones for the S'-wave diquarks, so that the amplitudes, differential and total 
cross sections were presented for the scalar-scalar, scalar-vector and vector-vecotr pairs. The details 



are given in the original ref. 56 1. For the representativity we show the expression for the total cross 



section of scalar pairs versus the square of total energy s 

a,, = 25671^;^ |vl>,(0)|Ml-^) , (3.11) 



9s3 
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where \l/s(0) is the wave function of diquark at the origin, and the formfactor has the form 



fs 



(3.12) 



so that qi^2 are the heavy quark charges, mi^2 are their masses, and M = mi + 1712- 

Numerical estimates show that the production of axial vector pairs dominates, so that in com- 
parison with the yield of heavy quark pairs the fraction of diquark pairs is about (2 — 6) ■ 10~^. 
This fact implies that, say, at the B factory among the charm events one should try to observe the 
events with the production of doubly charmed baryon pairs yielding a 10% fraction of their inclusive 
production. 



3.2. Perturbative fragmentation of diquark 

In this section we investigate the production of baryons in the fragmentation of vector and scalar 
particles interacting with quarks. From the QCD point of view the doubly heavy diquark of small 
size is a local color-triplet field, therefore, the results of this section can be used for the calculation of 
fragmentation into the baryons with doubly heavy vector and scalar diquarks. In this approach we 
explore the perturbative QCD for the calculation of hard amplitude in the fragmentation, which is 
factorized in front of soft amplitude for the forming of bound state. Sure, this method is quite accu- 
rate, if the hardness is provided by a large value of mass for the quark, which composes the hadronic 
state of baryon by coupling together with the diquark, for example, in the fragmentation of bb into 
bbc. However, the obtained expressions could be used as the QCD-motivated parametrizations for 
the processes with the light quarks, too. 

The fragmentation of scalar color-triplet local field was considered in ^7j. A new problem arising 
in the case of vector particle is a choice of the lagrangian for the vector diquark interaction with 
gluons. Indeed, to the lagrangian of a free vector field — tr [if^^^f/'^'^] , where H^i, = d^Ui, — d^U^, 
Ufj, is the vector complex field with derivatives substituted by covariant ones, we can add the gauge 
invariant term proportional to S^j^tr [G^'^U^Ua], where S^^ = 1/2(5^5^ — 5°5^) is the tensor of spin, 
G'^'^ is the gluon field strength tensor. It leads to the appearance of a parameter in the gluon- 
diquark vertex (the so-called anomalous chromo-magnetic moment). In this section we discuss the 
production of a 1/2-spin bound state containing the heavy vector particle at various values of this 
parameter. 

At high transverse momenta, the dominant production mechanism for the heavy baryon bound 
states is the diquark fragmentation, which can be calculated in perturbative QCD after the 
isolation of soft binding factor extracted from the nonrelativistic potential models . The corre- 
sponding fragmentation function is universal for any high energy process for the direct production 
of baryon. 

In the leading a^-order, the fragmentation function has a scaling form, which is the initial 
condition for the perturbative QCD evolution caused by the emission of hard gluons by the diquark 
before the hadronization. The corresponding splitting function differs from that for the heavy 
quark because of the spin structure of gluon coupling to the diquark, which is the vector or scalar 
color-triplet particle. 

In this work the scaling fragmentation function is calculated in the leading order of perturbation 
theory. The limit of infinitely heavy diquark, m^jg — >^ 00, is obtained from the full QCD consideration 
for the fragmentation. The distribution of bound state over the transverse momentum with respect 
to the axis of fragmentation is calculated in the scaling limit of perturbative QCD. The splitting 
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kernel of the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi-evolution (DGALAP) is derived, while the 
one-loop equations of renormalization group for the moments of fragmentation function are obtained 
and solved. These equations are universal, since they do not depend on whether the diquark will 
bound or free at low virtualities, where the perturbative evolution stops. The integrated probabilities 
of diquark fragmentation into the doubly heavy baryons are evaluated. 



3.2.1. Fragmentation function in the leading order 

The contribution of fragmentation into direct production of heavy baryon has the form 



dalEnip)] = / dz da[diq{p/ z), jj] Ddiq-^Eniz, fJ'), 
Jo 



'0 

where da is the differential cross section for the production of baryon with the 4-momentum p, da 
is that of the hard production of diquark with the scaled momentum p/z, and D{z) is interpreted 
as the fragmentation function depending on the fraction of momentum carried out by the bound 
state z. The value of /i determines the factorization scale. In accordance with the general DGLAP- 
evolution, the /x-dependent fragmentation function satisfies the equation 

dDdiq^^^{z,^) dy 



f dy 

J z y 



where P is the kernel caused by the emission of hard gluons by the diquark before the production of 
heavy quark pair. Therefore, the initial form of fragmentation function is determined by the diagram 
shown in Fig. |3.1| , and, hence, the corresponding initial factorization scale is equal to /i = 2mQ. 
Furthermore, this function can be calculated as an expansion in as(2mQ). The leading order 
contribution is evaluated in this section. 




Figure 3.1. The diagram for the fragmentation of diquark d\c\H into the heavy baryon 'Eh- 



Consider the fragmentation diagram in the system, where the momentum of initial diquark has 
the form q = (go, 0, 0, q^) and the baryon momentum is denoted by p, so that 



g2 = p2 = 
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In the static approximation for the bound state of diquark and heavy quark, the quark mass is 
expressed as rriQ = rM, and the diquark mass equals m = (1 — r)M = fM. The gluon-vector 
diquark vertex has the form 

TaiL'^ = -igst°'[gp,v{q + fp)^ - g^,a{,{l + ae)fp - siq)y - gya{,{l + ?e)q - aefp)^], (3.14) 

where ae is the anamalous chromo-magnetic moment, t'* is the QCD group generator in the funda- 
mental representation. The sum over the vector diquark polarizations with the momentum q (g^ = s) 
depends on the choice of the gauge of free field lagrangian (for example, the Stueckelberg gauge), 
but the fragmentation function is a physical quantity and has not to depend on the gauge parameter 
changing the contribution of longitudinal components of the vector field. So, in a general case, the 
sum over polarization can be taken in the form 

P{.q)^J.u = -g,,u + 

The matrix element of the fragmentation into the baryon with the spin of 1/2 has the form 



^ = - n ,,, Pi<lUP{fp),,T^^l^p^, - M)7VeH Ml (3.15) 

9vM^ rr[s — m^)^ ^ 



where the sum over the gluon polarization is written down in the axial gauge with n = (1, 0, 0, — 1) 



k ■ n 

and k = q — {1 — r)p. The spinors and q correspond to the baryon and heavy quark associated 
to the fragmentation. A4o denotes the matrix element for the hard production of diquark at high 
energy, -R(O) is the radial wave function at the origin. The matrix element squared and summed 
over the helicities of particles in the final state has the following structure: 

In the limit of high energies q ■ n oo W^y behaves like 

W^, = -g^,W + R^,, (3.16) 

where R^y can depend on the gauge parameters. After the separation of Lorentz structures it leads 
to scalar formfactor terms, which are small in comparison with W in the limit of g ■ n — oo. Define 

p ■ n 

z = . 

q ■ n 

The fragmentation function is determined by the expression 



If / m% \ 

^W-IS^/''^<^---T^)«- 

where W is defined in (|3.16 ). The integral in the expression for the fragmentation function diverges 
logarithmically at a constant value of anamalous chromo-magnetic moment if se does not equal 
— 1. We consider two sets for the behaviour of anamalous chromo-magnetic moment. The first is 
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38 = —1. Here we observe that the obtained fragmentation function coincides with that for the 
scalar di quark up to the factor of 1/3 



D(z) 



which tends to 



8al \R{0)\^ z^{l-z 



2437r M3r2f2 [1 - fz] 



D{y) 



■ js + Sr^ - (6 - lOr + 2r^ + 2r^)z + 
+ (3 - lOr + 14r2 - lOr^ + 3r^)-2^)}, 

Sal |i?(0)|My-l)' 



2437r 

at r and y = [1 — [1 — r)z) / {rz). 

D/N 



4y + 3y' 




1000 



0.85 



(3.17) 
(3.18) 



Figure 3.2. The fragmentation function of diquark into the heavy baryon, the A^-factor is determined by 
^ ~ 2437 M'^r'^{i~-r)^ ' fragmentation function at se = — 1 is shown by the dashed line, the fragmentation 
function at 1 + ee = 3M^/(s — m^-^) is given by the solid line (r = 0.02). 

The hmit of D{y) is in agreement with the general consideration of 1/m-expansion for the 
fragmentation function |6^, where 

D{y) = ^aiy) + biy). 

The dependence of a{y) on y has the same form as for the fragmentation of heavy quark into the 
quarkonium [p9| . 

The consideration of fragmentation at ae = — 1 + AM'^/{s — was done in ||6ll . 



The perturbative functions in the leading a^-order are shown in Fig. |3.2| at r = 0.02. We see 
hard distributions, which become softer with the evolution (see ref. |]57[| ). 
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3.2.2. Transverse momentum of baryon 



In the system with an infinite momentum of the fragmentating diquark its invariant mass is expressed 
by the fraction of longitudinal diquark momentum z and transverse momentum with respect to the 
fragmentation axis (see Fig. lOP as 



m + 



M2 



z{\ - z 



■[(1 



[l-r)zf + e] 



where t = px/M. The calculation of diagram in Fig. |3.1| gives the double distribution for the 
fragmentation probability 

(fP _ 

ds dz ' ' 

where at (ae = —1) the function V has the form 

' Slvr r^f^ [1 — r zY [s — 7ii?Y 



|rr^ _l_ ^(^\ + r — z{l + 4r — r^))- 



m 



zil-z) 



m 



2s 2 



(3.19) 



M2 ' 'V M2 

The distribution of baryon over the transverse momentum can be obtained by the integration of 

D{pT)lNt 
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Figure 3.3. The distributions over the transverse momentum with respect to the axis of diquark fragmentation 
into the baryon, A'f-factor is determined by Nt = jvf^r^(i-r) ^ at r = 0.02. The dashed line represents the 
result at SB = —1, the solid line does it at se = — 1 + 3M'^/{s — nidiq). 



(|3.19| ) over z 

2MH 



D{t) = [ dzV{z,s) 
Jo 



zil-zY 

Finally, we get quite a cumbersome expression presented in Appendix II. The characteristic form of 
distribution over the transverse momentum of baryon with respect to the axis of diquark fragmen- 



tation is shown in Fig. |3^ 
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3.2.3. Hard gluon emission 



The one-loop contribution can be calculated in the way described in the previous sections. This 
term does not depend on the part of diquark-gluon vertex with the anomalous chromo-magnetic 
moment, therefore the splitting kernel coincides with that for the scalar diquark. It equals 



diq^diq 



37r 



2x 



X 



where the "plus" denotes the ordinary action: dxf^{x) ■ g{x) = dxf{x) ■ [g{x] 
diquark splitting function can be compared with that of the heavy quark 



(3.20) 

-g{l)]. The 



Pc 



4as(/i) 



1 + X' 



2-, 



X 



which has the same normalization factor at a; — > 1. 

Furthermore, multiplying the evolution equation by z"' and integrating over z, one can get from 
( |3.13| ) the /i-dependence of moments a(„) up to the one-loop accuracy of renormalization group 



dain 



(n) 



1 1 

o + • • • + 7 

2 n + 1 



n > 1. 



(3.21) 



At n = the right hand side of ( p.21| ) equals zero, which means that the integrated probability 
of diquark fragmentation into the heavy baryon does not change during the evolution, and it is 
determined by the initial fragmentation function calculated perturbatively . 
The solution of equation (|3.21|) has the form 



«s(/i) 

as(/io) 



3/30 



(3.22) 



where one has used the one-loop expression for the QCD coupling constant 

2% 



PolniiJ,/ Aqcd) 



where = 11 — 2nf /3, with Uf being the number of quark flavors with nig < fi < mdiq. 

Relation ( |3.22| ) is universal one, since it is independent of whether the diquark is free or bound at 
the virtualities less than /iq. In this paper we take into account the evolution for the fragmentation 
into the heavy baryon. The diquark can lose about 20% of its momentum before the hadronization 



3.2.4. Integrated probabilities of fragmentation 

As has been mentioned above, the evolution conserves the integrated probability of fragmentation 
which can be calculated explicitly 



dz D{z) 



Sal \Rm 



Sivr IQmi 



w(r). 



(3.23) 
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Figure 3.4. The w functions for the diquark fragmentation into the heavy baryon versus the fraction r = mq/M . 
The dashed curve corresponds to the case of ae = —1, while the solid line does to se = —1 + 3M^/(s — m\ 



diq) 



SO that 



w(r) 



16 [(8 + 15r - 60r2 + lOOr^ - 60r^ - 2>r^) + 30r(l 



15(1 



(3.24) 



The w(r) functions are shown in Fig. |3.4| at low r. 

Thus, we have considered the dominant mechanism for the production of bound spin 1/2 states 
composed by a local color-triplet vector field with a heavy anti-quark for high energy processes 
at large transverse momenta, where the fragmentation contributes as the leading term. We have 
investigated two cases for the behaviour of anamalous chromo-magnetic moment. At ae = — 1 we 
observeP^ that the obtained fragmentation function coincides with that for the scalar diquark up to 
a factor. In the infinitely heavy diquark limit, D{z) has the form, which agrees with what expected 
from general consideration of 1/m-expansion for the fragmentation function. The distribution of 
bound state over the transverse momentum with respect to the axis of diquark fragmentation is 
calculated in the scaling limit of perturbative QCD. The hard gluon corrections caused by the 
splitting of vector diquark has been taken into account in the perturbation theory, that has led to 
the corresponding one-loop equations of renormalization group for the moments of fragmentation 
function (see ( ^1211) , (^1221)). 

The numerical estimates show that the probabilities of fragmentation into the bound states 
containing the heavy vector diquark with the mass from 3 to 10 GeV depend on the effective mass 
of quark entering the baryon together with the diquark. In this way the ratio of yields for the baryons 
with the strangeness and without it, is approximately equal to aiVtQQi) / aip^qqi) ^ 0.2. Certainly, 
the introduction of light and strange quarks into the consideration makes the obtained results to be 
not strictly justified, since the calculations suggest that the constituent masses lead to the effective 

^^The expression for the fragmentation function diverges logarithmically at a constant value of anamalous chromo- 
magnetic moment if ae is not equal to —1. 
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and correct description of dominant contributions caused by the infrared dynamics. However, we 
can use the perturbative expressions as models for the fragmentation into the hadrons containing 
hght and strange quarks because in such processes we describe the "fast" valence degrees of freedom 
in the baryon. In the limit of small dispersion these degrees of freedom can be approximated by 
the introduction of ratio for the assigned fractions of longitudinal moment for the parton inside the 
hadron, while we do not provide the consideration of contribution by a soft sea of light quarks and 
gluons. 

Another approach to the fragmentation mechanism for the production of doubly heavy baryons 
was considered in ref. wherein the perturbative formfactors of doubly heavy and heavy-light 



diquarks were calculated and the heavy quark fragmentation function into the baryons was derived 
in the process with the production of vector diquark pairs: Q Hqi^qiq)^^^ + {Q'q)diq- In this way 
the estimates strictly mean the minimal expectations, since they are based on the elastic formfactor 
of diquark. We also believe that the hierarchy of scales rriQ ^ mq -v ^ ^qcd in the relevant strong 
interactions provides that after the hard production of heavy quark (the virtualities about mg) quite 
a fast forming of the doubly heavy diquark occures (the virtualities of the order of mq ■ v), while a 
more slow hadronization of diquark into the baryons (the virtualities sloe to Aqco) takes place in 
the final stage. V.Saleev analyzed the fragmentation into the thriply charmed baryon Qccc due 



to the cascade process of quark into the diquark and diquark into the baryon in comparison with 
the direct fragmentation of quark into the baryon due to the elastic production of vector diquark. 



Unfortumately, the general conclusion of [^] on the significant dominance of direct fragmentation 
is not correct, because, first, the use of elastic formfactor for the vector diquark in the cascade 
process leads to the double counting of suppressing factor due to the charmed quark fusion into the 
diquark in the fragmetation of quark into the diquark as well as in the elastic formfactor, where 
one again took the projection of initial quark state on the bound diquark. This procedure gives 
the overcoming factor of suppression ag|\E'cc(0)P/'m.^ ~ 10~^. Second, the idea on the forming the 
cc-diquark with the consequent hard production of charmed quark off diquark is not correct because 
the time period (the diquark size) for forming the diquark is much greater than the time interval 
(the compton length) for the production of charmed quark. 

3.3. Hadronic production 

Recent years are marked by a rapid increase of charmed particles observed in modern experiments. 
So, the study of about 10^ charmed particles is expected at fixed target FNAL facilities of E831 
and E781. An increase of this value by two orders of magnitude is proposed in experiments of 
next generation. Along with standard problems of CP-violation in the charmed quark sector and a 
measuring of rare decays etc., an investigation of processes with more than one cc-pair production 
becomes actual. The production of additional cc-pair strongly decreases a value of cross section for 
such processes. This fact must be especially taken into account in fixed target experiments, where 
the quark-partonic luminosities are strongly suppressed in the region of heavy mass production. 

An interesting process of the mentioned kind is the doubly charmed baryon production. The 
doubly charmed Scc''-baryon represents an absolutely new type of objects in comparison with the 
ordinary baryons containing light quarks only. The ground state of such baryon is analogous to 
a Q^-meson, which contains a single heavy anti-quark Q and a light quark q. In the doubly 
heavy baryon the role of heavy anti-quark is played by the cc-diquark, which is in anti-triplet color 
state [Q. It has a small size in comparison with the scale of the light quark confinement. The 
study of mechanism for the production of these states are of interest. The ccg-baryon production 
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was discussed in 0,^,^,^|6^. The main problem of calculations is reduced to an evaluation of 
the production cross section for the diquark in the 3-triplet color state. One assumes further that 
the cc-diquark nonperturbatively transforms into the ccg-baryon with a probability close to unit. 

The hadronic production of diquark is subdivided into two parts. The first stage is the hard 
production of two cc-pairs in the processes of gg cccc and qq — >• cccc described by the Feynman 
diagrams in the fourth order of as coupling constant in the perurbative QCD (see Fig. WM - 



gg 



(cch + X 










qq 



(cch+X 







Figure 3.5. The examples of diagrams for the gluon-gluon and quark-antiquark production of cc-diquark. The 
initial quarks are denoted by the thin fermion lines, the final quarks are denoted by the bold fermion lines and 
the gluons are denoted by the helical lines. 

The second step is the nonperturbative fusion of two c-quarks with a small relative momentum 
into the cc-diquark. For the S'-wave states this process is characterized by the radial wave function 
at the origin, R{0). 

The main difference between the existing evaluations of the doubly charmed baryon cross section 
consists in the methods used for the hard subprocess calculation. In a part of diagrams connected 
with the c-fragmentation into the (cc)-diquark is only used instead of the complete set of diagrams. 
As was shown in |^5|] this estimation is not absolutely correct, because it becomes true only at large 
transverse momenta greater than p^*" > 35 GeV, where the fragmentation mechanism is dominant. 
In other kinematical regions the application of fragmentational approximation is not justified and 
it leads to wrong results, especially at \/S being not much greater than p™". 

However, even after taking into account the complete set of diagrams considered in [p5|,|5^, 
essential uncertainties in the estimations of the ccg-baryon production remain. The basic parameters 
determining these uncertainties are the values of as, rric and Rcc{0). In addition, it is not clear, to 
what extent the hypothesis on the hadronization of cc-diquark into the ccg-baryon with the unit 
probability is correct or not. The matter is that the interaction between the diquark and gluons is 
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not suppressed in contrast to the cc-pair production in the color singlet state, when the quarkonium 
dissociation supposes exchange with the quark-gluonic sea by two hard gluons with virtuahties, 
which are greater than the inverse size of quarkonium. 

A decrease of the uncertainty in the ccg-baryon cross section would be possible by means of 
comparing the process of baryon production with the analogous process of J /'^ + DD production]^. 
The associated production of J/\E' is described by practically the same diagrams of fourth order with 
the well-known wave function of J/\E' at the origin]^. In this way of connection to the J/\l/ + DD 
process we could remove the part of uncertainties following from and in the cc-diquark 
production process. 

In the following sections of the paper the joint cross section calculations of these processes in 
7r~p and pp interactions are performed, the description of production models for the ccg-baryons 
and J /'^ + DD is given. We presents the calculation results for the production cross section of 
ccg-baryons and J/\I' + DD in the fixed target experiments E781 and HERA-B as well as at the 
collider energies of Tevatron and LHC. We discuss possibilities of searching for the baryons Sc*'*. 



3.3.1. Production mechanism 

As was mentioned, we suppose that the diquark production can be subdivided by two stages. On 
the first stage the production amplitude of four free quarks is calculated for the following processes 



gg cccc, 
qq — > cccc. 



The calculation technique applied in this work is analogous to that for the hadronic production 
of Be [^,0, but in this case the bound state is composed by two quarks {Qi and Q2) instead of 
the quark and ant i- quark |^,^,^ . 



One assumes that the binding energy in the diquark is much less than the masses of constituent 
quarks and, therefore, these quarks are on the mass shells. So, the quark four-momenta are related 
to the {Q1Q2) diquark momentum -P(QiQ2) in the following way: 

PQi = MQ1Q2) ) PQ2 = T7 P{QiQ2)y (3-25) 

where Mi^q^q^) = mq^ + mq^ is the diquark mass, mg^, mq^ are the quark masses. 

In the given approach the diquark production is described by 36 Feynman diagrams of the 
leading order, corresponding to the production of four free quarks with the combining of two quarks 
into the color anti-triplet diquark with the given quantum numbers over the Lorentz group. The 
latter procedure is performed by means of the projection operators 



2M(Q,Q,) 1 



A/'(0, 0) = W 2mQ, +)^2(PQ., -) - ^i(PQ„ -)«2(PQ„ +)} (3.26) 



^^We calculate the J /ij: + cc production and assume that the cc-pair transforms to DD+ some light hadrons with 
the probability very close to unit. So, we neglect the production of charmed baryons as well as bound states of 
charmonium in the hadronization of associated cc pair. 

^■^The value of |i?;p(0)| is determined by the width of leptonic decay, J/'f l^l^ with taking into account the 
hard gluonic correction, so, numerically, |i?^(0)| — ^JuM/if^, where — 540 MeV. 
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for the scalar state of diquark (the corresponding baryon is denoted by 'E.'q^q_^{J = 1/2)); 



2^(QiQ2) 



2M(Q,Q,) 1 



" ^^ 2m 2m' 7^{^i(PQi' +)"2(pq2, -) + Mi(PQi, -)^i2(pQ2, +)}, 



2^(QiQ2) 



for the vector state of diquark (the baryons are denoted as 'Eq^q^{J = 1/2) and Sq^q^(J = 3/2)). 

To produce the quarks, composing the diquark in the 3c state, one has to introduce the cofor 
wave function as eijk/V^, into the diquark production vertex, so that z = 1,2,3 is the cofor index 
of the first quark, j is that of the second one, and k is the cofor index of diquark. For the identical 
quarks Qi = Q2 possessing the momenta equal to each other, the anti-triplet color state can have 
the summed spin S = 1, only. 

The diquark production amplitude ^4^^^^ is expressed through the amplitude T^^^{pi) for the free 
quark production in kinematics (|3.25| ) under the substitution for the product U1U2 by the projection 



operators as well as under the condition that two heavy quarks in the 3-color state, so that 

Af- = ^^-gL^°^ Tf-(p,.), (3.28) 
v47r 



where -RqiQ2(0) is the diquark radial wave function at the origin, k is the color state of diquark, S 
and Sz are the diquark spin and its projection on the 2;-axis, correspondingly. 

In numerical calculations under discussion we suppose the following values of parameters: 

= 0.2, = 1.7 GeV, = 4.9 GeV, 
i?cc(i5)(0) = 0.601 GeV^/', Rbciis){0) = 0.714 GeV'^', (3.29) 

where the value of -Rcc(O) has been calculated by means of numerical solution of the Schrodinger 
equation with the Martin potential multiplied by the 1/2 factor caused by the color anti-triplet 
state of quarks instead of the singlet one. 

To calculate the production cross section of diquarks composed of two c-quarks, one has to 
account for their identity. One can easily find, that the anti-symmetrization over the identical 
fermions leads to the scalar diquark amplitude equal to zero, and it results in the amplitude of the 
vector cc-diquark production being obtained by the substituting of equal masses in the production 
amplitude of vector diquark composed of two quarks with the different flavors, and taking into 
account the 1/2 factor for the identical quarks and anti-quarks. 

In this approach we suppose that the produced diquark forms the baryon with the unit prob- 
ability by catching up the light quark from the quark- ant iquark sea at small px or having the 
fragmentation into the baryon at large px- 

The typical diagrams of fourth order describing the parton processes are shown in Fig. 



One can subdivide them into two groups. The first group contains the diagrams of fragmentation 
type, wherein the cc-pair emits another one. The second group corresponds to the independent 
dissociation of gluons into the cc-pairs with the consequent fusion into the diquark. The diagrams 
of second group belong to the recombination type. 
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As was mentioned above, the authors of some papers restricted themselves by the consid- 
eration of fragmentation diagrams, only. In this way they reduced the cross-section formulae to 
the cc-pair production cross section multiplied by the fragmentation function of c-quark into the 
cc-diquark. 



Nda^pclpT, pb/GeV 




Figure 3.6. The differential cross section for the associated production of cc diquarkin the gluon-gluon sub- 
process at 100 GeV (histogram) in comparison with the prediction of fragmentation model (solid curve), corre- 
spondingly. 



As we shown in [35|, this approximation is correct only under the two following conditions: 



Mqq -C s and px ^ Mqq. In other kinematical regions, the contribution of recombination dia- 
grams dominates. The typical value of pt, wherefrom the fragmentation begins to dominate in the 
production of -diquark, is px > 35 GeV (see Fig. ^.61 ). It is clear, that at reahstic pt one has to 
take into account all contributions including the recombination. For the first time, the complete set 
of diagrams was taken into account in |6^ and confirmed in [Q. In the both papers the calcula- 
tions are performed only for the gluon-gluon production, which is a rather good approximation at 
collider energies. For the fixed target experiments the value of total energy strongly decreases, and, 
hence, the values of energy in the parton subprocesses decrease too. The contribution of quark- 
antiquark annihilation becomes essential at fixed target energies, especially for the processes with 
initial valence anti-quarks. In the following consideration we take into account the quark-antiquark 
annihilation into four free charmed quarks in the estimation of yield for the doubly charmed baryon. 



3.3.2. Doubly charmed baryon production in fixed target experiments 



In Figs KTl and 



the calculation results for the total cross section of the diquark-production 

subprocesses versus the total energy for the given values of rric and i?cc(0). We have also 
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shown the parametrizations for the dependencies of total cross sections versus the process energy 
for the cc-diquark as given by 



a(-) = 213. f 1 - ^ V'Y^V" Pb, 



99 



= 206. I 1 - 



pb. 



(3.30) 

(3.31) 



We have to stress that the numerical coefficients depend on the model parameters, so that a 
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Figure 3.7. The total cross sections for the gluon-gluon production of cc diquark (a) and J/^ + DD (a) in 



comparison with the approximations of ( p. 30] ) and ( 3.32| ) (solid and dashed lines, correspondingly) 



As was mentioned, the production of J/"^ in the subprocesses of gg J/'^+cc and qq J/^+cc 
is also calculated in this work. The numerical results of such consideration are shown in Figs. |3] 
and p.8|. The parameterization of these results versus the energy \fl are presented below 



a/f = 518. f 1 



0", 
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pb, 
pb. 



(3.32) 
(3.33) 



As well as for the associated B^. + he and Sec + cc, we see the following regularity for the partonic 
production of J/iJ + cc: the fragmentation regime occures at > 25 — 30 GeV. This fact can be 



66 



a- pb 




/s^, GeV 



20 40 60 80 100 120 
Figure 3.8. The total cross sections for the gluon-gluon production of cc diquark (a) and J/'^ + DD (a) in 



comparison with the approximations of ( [3. 311) and (|3.33| ) (solid and dashed lines, correspondingly) 



certainly observed in Fig. |3.9| for the differential cross section oi gg ^ J / Tp + cc y s = 100 GeV. 
Thus, for the associated production of J/ip + cc and Sec + cc the fragmentation works at pt ^ "^c- 
The above parametrizations quite accurately reconstruct the results of precise calculations at 
\/I < 150 GeV, and that is why they can be used for the approximate estimation of total hadronic 
production cross section for the cc-diquark and J/\E' by means of their convolution with the partonic 
distributions 

^ = ^ dxidx2fi/A{xi,fi)fj/B{x2,f^)^, (3.34) 
id 

where fi/A{x, fi) is the distribution of i-kind parton in the A-hadron. The parton distributions used 
for the proton are the CTEQ4 parameterizations and those of used for the 7r~-meson are the 
Hpdf ones [70|. In both cases the virtuality scale is fixed at 10 GeV. As for the choice of fixed 



scale in the structure functions, this is caused by the fact that the cross section of subprocesses 
is integrated in the region of low s close to the fixed scale, so that the account of "running" scale 
weakly changes the estimate of S^c-baryon yield in comparison with the mentioned uncertainty of 
diquark model[^. The cross sections convoluted with the gluon and quark luminosities are presented 
in Fig. p.lO| for both the cc-diquark and J/\I' production in ir'p and pp collisions. 

As we can see in these figures, the cross section of cc-diquark as well as the cross section of 
J/'^ + cc are strongly suppressed at low energies in comparison with the values at the collider 
energies. The ratio for the cc-diquark production and total charm production is dcc/cTcharm ~ 
10~^ — 10"'^ in the collider experiments and ~ 10~^ — 10~^ in the fixed target experiments. The 

^"'We have found the scale-dependent variation to be at the level of 6a/a ~ 10%. 
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Figure 3.9. The differential cross section for the associated production of J/"^ + cc in the gluon-gluon subpro- 
cess at 100 GeV (solid histogram) in comparison with the prediction of fragmentation model (dashed curve), 
correspondingly. 

same situation is observed for the hadronic J/\l/ + DD production. 

The distributions for the ccg-baryon and J/\[^ + DD production are shown in ||71| for the 7r~p- 
interaction at 35 GeV and for the pp-interaction at 40 GeV, correspondingly. The rapidity dis- 
tributions point to the central character of ccg-baryon production as well as J/'^ + DD. The 
PT-distributions of these processes are alike to each other alsoj^. The form of differential distribu- 
tions leads to the conclusion that the process of J/\E' + DD production can be used to normalize 
the estimate of ccg-baryon yield, wherein the following additional uncertainties appear: 

1. the unknown value of |-R(cc)(0)p, 

2. uncertainties related with the hadronization of cc-diquark. 

We see from the given estimates that in the experiments with the number of charmed events at 
the level about 10^ (for example, in the E781 experiment, where a/s = 35 GeV), we have to expect 
about one event with the doubly charmed baryon. The situation is more promising and pleasant for 
the pp-interaction at 800 GeV (HERA-B). The considered processes yield about 10^ sic''-baryons 
and a close number of J/'^ + DD in the experiment specialized for the detection of about 10^ events 
with the 6-quarks. 

^^We assume that at the given energies the cc-diquark has no fragmentational transition into the baryon, but it 
catches up the hght quark from the quark-antiquark pair sea. 
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Figure 3.10. The total cross sections for the production of cc-diquark and J /'^ + DD (solid and dashed curves, 
correspondingly) in ■K~p and pp collisions. 



3.3.3. Production of ccg-baryon at colliders 

As one can see in the previous Section, the observation of Scc'*-baryons presents a rather difficult 
problem in the experiments specialized for the study of charmed particles. As a rule, such exper- 
iments are carried out at fixed targets, so that the effective value of subprocess energy is strongly 
decreased. So, the relative contribution of doubly charmed baryons into the total charm yield is of 
the order of 10^^ — 10^^. The production of ccg-baryons at colliders with large px is more effective. 
In this case the cross section is determined by the region of quark-antiquark and gluon-gluonic 
energy, where the threshold effect becomes negligible and the partonic luminosities are quite large 
at a; ~ M/ -y/i. So, the suppression factor with respect to the single production of cc-pairs is much 
less and it is in the range of 10~^ — 10"'^. 

In [71 1 the pT-distributions for Hcc^ and J/\l' associated with D and D are shown at the energies 
of Tevatron and LHC with the rapidity cut \y\ < 1. 

One can easily understand that the presented si*'' cross sections are the upper estimates for the 
real cross sections because of the possible dissociation of heavy diquark into the DD-pair. 

Further, even if the cc-diquark being the color object, transforms into the baryon with the unit 
probability, one has to introduce the fragmentation function describing the hadronization of diquark 
into the baryon at quite large pt values. The simplest form of this function can be chosen by the 
analogy with that for the heavy quark 



Diz) 



1 



1 



(M2 - 



(3.35) 



where M is the mass of baryon S, 



CC ) 



is the mass of diquark, rriq is the mass of fight quark (we 



suppose it to be equal to 300 MeV). This function practicafiy repeats the form of fragmentation 
function if diquark into the baryon, that has been derived above in this chapter in the framework of 
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perturbative QCD. In ||7T1 the pT-distributions of doubly charmed baryon production are calculated 
with the use of ( ^.351 ). One has to mention, that in the leading order over the inverse heavy 
quark mass, the relative yield of S^c and S*^ is determined by the simple counting rule for the spin 
states, and it equals a{Ecc) '■ cr{E*^) = 1:2. In this approach one does not take into account a 
possible difference between the fragmentation functions for the baryons with the different spins. 
The corresponding difference is observed in the perturbative fragmentation functions for the heavy 



mesons and quarkonia [pl|,[59 



3.3.4. Hadronic production of S 



be 



The total energy dependence of the gluonic production cross sections of E[ 



be 



5) and H^*^ 



be 



baryons 

is shown in Fig. |3.11| . For the sake of comparison, the predictions of the fragmentational mechanism 
for (solid line) and (dashed line) are also presented. One can see from the figure, that the 
fragmentational production mechanism assuming the validity of factorization in the cross section at 
M^/ s <^ 1 and at large transverse momenta in accordance with the formula 



da „, , 



,J*hbc 

bc^ be ' 



dz 



(3.36) 



where z = 2\P\/^/s, does not work at low gluon energies, where it overestimates the cross section, 
because of the incorrect evaluation of the phase space, and it is also not valid at large energies, 
where the predictions of the fragmentational mechanism are essentially less that the exact results. 
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Figure 3.11. The gluonic production cross sections of 



of the fragmentational mechanism for H'j,^ (dashed line) and h[* (solid line) 
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(•) in comparison with the predictions 



70 



So, the fragment at ional values underestimate the Ei'J and cross sections by 10 and 3 times, 
respectively, at = 100 GeV. When the fragmentational predictions give the ratio cr„(.) /cth' ~ 1.4, 

the exact perturbative calculations result in cr„{.)/cr=' ^ 3.5 even at = 100 GeV. 

The agreement with the fragmentational production at \/I = 100 GeV takes place at large 
transverse momenta of the baryon, as one can see from the distributions over px for the H^*"* and 
S'^^ production, shown in Fig. |3.12| . Note, that in contrast to the doubly heavy baryon production, 
the exact perturbative calculations of the gluonic production of Bc{B*) mesons with pt > 35 GeV 
at VS = 100 GeV agree with the fragmentational predictions. For the baryon production, a visible 
deviation is observed up to the largest values of pt, though it is evident that at larger energies of 
gluon-gluon subprocess the region of justified use of fragmentation regime will become wider in the 
direction of large transverse momenta. 

Pb/GeV 



-I 

10 - 





Pt, GeV 

Figure 3.12. The distributions over the transverse momentum in the gluonic production of S'f,j,(s[*^) in com- 
parison with the fragmentation result at the interaction energy of 100 GeV. The solid line corresponds to the 
production of , the dashed curve represents the yield of r!^^, while the complete results and fragmentation 
estimates are given by the histogram and smooth curves, correspondingly. 

The differential cross section da/dpT of the S'^^ and production in pp interactions at ^/s = 1.8 



TeV is presented in [^] in comparison with the fragmentational predictions. So, we can also 
conclude that the fragmentation approach gives rather rough estimate for the yield and momentum 
distribution of Ebc baryons. 

At the chosen values of parameters and with the account for the cuts over the transverse mo- 
mentum and rapidity of the baryons {pt > 5 GeV and \y\ < 1), the production cross section of the 
IS-wave 6cg-baryons and its anti-particles is evaluated as (Xbcg — 1 nb (without cuts the value of 
abcq is about two times greater). After the expected end of Run lb at Tevatron with the integral 
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luminosity 100 -r- 150 pb , one has the yields of 1.0 -j- 1.5 ■ 10^ of 6cg-baryons. 
3.3.5. Pair production 

At the energies of fixed target experiments the luminosity of parton subprocesses with the valence 
quarks is not suppressed in comparison with the luminosity of gluon-gluon collisions in the region 
of large invariant masses. In this way a significant fraction of total cross section for the production 
of baryons with two heavy quarks is given by the pair production of such baryons. The total and 
differential cross sections for the pair production were considered in [^], wherein the contributions 



of scalar and axial- vector diquarks were taken into account. So, the expression for the total cross 
section of scalar pairs versus the square of total energy s has the form 



81s3 ' \ s V 3 20 



^40)r 1 - ^ ^(/W)^ + ^ 1 - ^ Cff}f , (3.37) 



where ^l/s(0) is the diquark wave function at the origin, while the formfactors are determined by the 
functions 



J ss 


= Mai 








m 

J ss 













2M2 /m2 mi 
m^J s \mf m\ 



(3.38) 
(3.39) 

so that mi^2 are the heavy quark masses, and M = mi + 1712. 

Numerical estimates show that the pair production of vector diquarks dominates, so that it gives 
about 10% of single doubly heavy baryon production in the parton subprocess. 

3.3.6. Discussion 

On the basic of perturbative calculations for the hard production of doubly charmed diquark frag- 
mentating in the baryon we have shown that the observation of doubly charmed baryons is a difficult 
problem, because the ratio of a{E!"*c) / (y{charm) for these baryons and charmed particles yields the 
value of 10^^ — lO^^ depending on the process energy. The suppression of doubly charmed baryon 
yield at low energies of fixed target experiments is explained by the threshold effect. As one can 



Table 3.2. The total cross section of doubly charmed baryon production at different facilities. 



facility 


HERA-B 


E781 


Tevatron 


LHC 


total cross section, 
nb / nucleon 


2 ■ 10-3 


4.6 ■ 10-3 


12 


122 



see in Table p.2| , the low value of cross section for the production of doubly charmed baryons in the 
fixed target experiments allows us to expect about 10^ events with the production of these baryons 
at HERA-B. Practically the same number of events at pt > 5 GeV and \y\ < 1 is expected at 
Tevatron with the integrated luminosity of 100 pb~^. The large luminosity and large interaction 
energy allow one to increase the yield of the doubly charmed baryon by 10^ times at LHC. 
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Taking into account the increase of luminosity of collider at FNAL, we can believe that the 
experimental search for the events with the baryons and Sec is the actual task. 

Under conditions of large yield of the doubly charmed baryons, the problem of their registration 
is under challenge. 

First of all it is interesting to estimate the lifetimes of the lightest states of and The 
simple study of quark diagrams shows that in the decay of S+^-baryon the Pauli interference for 
the decay products of charmed quark and valence quark in the initial state takes place as well as in 
the case of D+-meson decay. In the decay of S+ the exchange by the M^-boson between the valence 
quarks plays an important role as well as in the annihilation of D^. Therefore we suppose that the 
mentioned mechanisms]^ give the same ratio for the both baryon and D-meson lifetimes 

r(S+)^0.4-r(H++). 

Further, the presence of two charmed quark in the initial state results in following expressions: 

r(S++) ^ ir(D+) ^ 0.53 ps, 

r(H+) ^ ir(D°) ^ 0.21 ps. 

Following the analogy with the decays of baryon containing a single heavy quark, we can expect 
the decay modes 

BR(S++ ^ irOWs++(*)) ^ BR(S+ i^o{*)(s+(*) + A+)) ^ 
^ BR(Ac ^ ^ 4 ■ 10"^ 

One can observe 4 ■ 10^ events in these decay modes at HERA-B and Tevatron without taking 
into account a detection efficiency. One has to expect the yield of 4-10^ such decays at LHC. Among 
other decay modes, the —>■ 7r+S+ and S+ —>■ 7r"'"S° processes taking place with branching ratios 
of about 1%, can be essential. 

The excited S*^ states always decay into Sec by the emission of 7-quantum, so that the branching 
fraction of transition is equal to 100%, since the emission of vr-meson is impossible in the S*^ decay 
because of the small value of splitting between the ground state and the excited one, in contrast to 
the charmed meson decay. 

In conclusion we mention another possibility to increase the yield of doubly charmed baryons in 



fixed target experiments. In the model of intrinsic charm ||72[ one assumes, that the nonperturbative 
admixture of exotic hybrid state \ccuud) presents in the proton along with the ordinary state \uud) 
including three light valence quarks. The probability Pic of |ccuu(i)-state is suppressed at the level of 
1%. The valence charmed quark from that state can recombine with the charmed quark produced 
in the hard partonic process with the cc-pair production. The energy dependence for such doubly 
charmed baryon production repeats one for the single charmed quark production in the framework 
of pQCD up to the factor of exotic state suppression and the factor of fusion of two charmed quarks 
into the diquark, K ~ 0.1. This mechanism has no threshold of four quark state production in 
contrast to the discussed perturbative one. Therefore at low energies of fixed target experiments, 
where the threshold suppression of perturbative mechanism is strong, the model of intrinsic charm 



^^We present quite naive estimates for the lifetimes, while a strict consideration of operator product expansion for 
the inclusive decays of doubly heavy baryons is addressed in Chapter 4. 
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would yield the dominant contribution into the Hcc production. So, the number of events in this 
model would be increased by three orders of magnitude, and ratio of si*^ and charmed particle yields 
would equal (j(sic^)/cr (charm) ~ 10~^. At high energies the perturbative production is comparable 
with the intrinsic charm contribution. We have to note that the \ccccuud) -state suppressed at 
the level of 3 ■ 10~^, also could increase the doubly charmed baryon production at low energies of 
hadron-hadron collisions. 

Thus, the observation of S*^-baryons in hadronic interactions is a quite realistic problem, whose 
solution opens new possibilities to research the heavy quark interactions. The observation of Sc*''- 
baryons at fixed target experiments |]73[ would allow one to investigate the contributions of different 
mechanisms in the doubly charmed baryon production, as the contribution of the perturbative 
mechanism and that of the intrinsic charm, which strongly increases the yield of these baryons. 
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Chapter 4. Lifetimes and decays of Sqq' baryons 

In the framework of Operator Product Expansion (OPE) over the inverse powers of heavy quark 
mass there is a definite scheme for consistent calculations of QCD effects, that was developed 
for the evaluation of various characteristics in decays of hadrons containing the heavy quarks |P, 
^|7^. The consideration based on the predictions of this approach allows one to extract the 
parameters assigned to the electroweak interactions of heavy quarks in the background of dynamics 
with strong interactions of quarks and gluons, which form the hadrons observed in the experimental 
measurements. The accuracy of QCD description in the sector of heavy quarks has an important 
significance in the search for subtle effects such as the violation of CP-invariance, deviations from the 
predictions of Standard Model as well as in clarifying the mechanisms for contributions of virtual 
corrections caused by a "new" physics operating on characteristic energy scales of TeV range. 
Therefore, the study of OPE over the inverse powers of heavy quark mass is quite an informative 
problem of interest serving a manyfold attention. An importamt challenge is a complex investigation 
of system containing the heavy quarks with the analysis and comparison of various characteristics 
such as l) the convergency of expansion in both the inverse mass and the QCD coupling constant, 
ll) relative and absolute values of various contributions and their dependence on the quark contents 
of systems, ill) qualitative conclusions on the role of some mechanisms |]75[, and iv) uncertainties 
of numerical estimates. 

The efficiency of approach under discussion was convincingly shown in the description of weak 
decays for the hadrons having a single heavy quark, carried out in the framework of the heavy 
quark effective theory (HQET) Q], in the annihilation and radiative decays of heavy quarkonia QQ, 
using the framework of nonrelativistic QCD (NRQCD) [l^, and in the weak decays of long-lived 
heavy quarkonium with mixed heavy flavour^ 5^ . We emphasize that the experimental data 
on the weak decays of doubly heavy hadrons could be able to bring a significant improvement of 
numerical accuracy in the parameters entering the description of systems with the heavy quarks. The 
advantage is caused by the presence of an additional small parameter in the NRQCD expansion in 
contrast to HQET. This parameter is a relative velocity of heavy quark motion v. Moreover, there is 
an essential variability of characteristics assigned to the bound states with the heavy quarks, so that 
these changes of properties allow us to investigate the dependence of OPE on the nonperturbative 
parameters, which can be evaluated, for example, in the potential approach. 

The baryons with two heavy quarks, QQ'q, provide a new insight in the description of systems 
containing the heavy quarks. For these baryons we can apply a method based on the combined 
HQET-NRQCD techniques [§,|l3,|7^, if we use the quark-diquark factorization for the bound states. 



The expansion in the inverse powers of heavy quark mass for the heavy diquark QQ' is a straight- 
forward generalization of these techniques in the mesonic decays [I^jH with the difference that we 
deal with the color anti-triplet system of heavy quarks with the appropriate account of interaction 
with the light quark instead of the color singlet systems. The HQET methods have to reliably work 
for describing the interaction of local diquark with the light quark. 

In this chapter, we present the consistent calculation of lifetimes for the doubly heavy baryons 
S^, S+ and 2++. Taking into account necessary generalizations to the case of hadrons with 
two heavy quarks and other corrections in the description of inclusive decays for the baryons, 
we follow the general consideration of heavy hadron decays in P,|7^, where the decays of the 



hadrons with a single heavy quark and the doubly heavy meson were discussed. The justified 



^^The first experimental observation of Bc-meson has been reported by the CDF Collaboration while a review 
on the theoretical status of Be is given in M . 
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basis for such calculations is the optical theorem for the inclusive decay width combined with the 
Operator Product Expansion (OPE) for the transition currents in accordance with the consequent 
nonrelativistic expansion of hadronic matrix elements derived in OPE. Using OPE at the first 
step, we exploit the fact, that, due to the presence of heavy quarks in the initial state, the energy 
release in the decay of both quarks is large enough in comparison with the binding energy in the 
state. Thus, we can use the expansion over the ratio of these scales. Technically, this step repeats 
an analogous procedure for the inclusive decays of heavy-light mesons as it was reviewed in |76| . 



Exploring the nonrelativistic expansion of hadronic matrix elements at the second step, we use the 
approximation of nonrelativistic QCD |TD[, which allows one to reduce the evaluation of matrix 
elements for the full QCD operators, corresponding to the interaction of heavy quarks inside the 
diquark, to the expansion in powers of where pc = rricVc ~ 1 GeV is a typical momentum of the 
heavy quark inside the baryon. The same procedure for the matrix elements, determined by the 
strong interaction of heavy quarks with the light quark, leads to the expansion in powers of . 

We note, that keeping just the leading term in the OPE, the inclusive widths are determined by 
the mechanism of spectator decays involving free quarks, wherein the corrections in the perturbative 
QCD are taken into account. The introduction of subleading terms in the expansion over the inverse 
heavy quarks mass0 allows one to take into account the corrections due to the quark confinement 
inside the hadron. In this way, an essential role is played by the following nonperturbative charac- 
teristics: the motion of heavy quark inside the hadron and the corresponding time dilation in the 
hadron rest frame with respect to the quark rest frame, and the influence of the chromo-magnetic 
interaction of the quarks. The important ingredient of such corrections in the baryons with two 
heavy quarks is the presence of a compact heavy diquark. The next peculiarity of baryons with 
two heavy quarks is the significant numerical influence on the lifetimes by the quark contents of 
hadrons, since in the third order of inverse heavy quark mass, l/mg, the four-quark correlations 
in the total width are enforced in the effective lagrangian due to the two-particle phase space in 
the intermediate state (see discussion in 0). In this situation, we have to add the effects of the 
Pauli interference between the products of heavy quark decays and the quarks in the initial state 
as well as the weak scattering involving the quarks composing the hadron. Through such terms 
we introduce the corrections involving the masses of light and strange quarks in the framework of 
nonrelativistic models with the constituent quarks. We include the corrections to the effective weak 
lagrangian due to the evolution of Wilson coefficients from the scale of the order of heavy quark 
mass to the energy, characterizing the binding of quarks inside the hadron. 

Following the picture given above, we describe the general scheme for the construction of OPE for 
the total widths of baryons with two heavy quarks with account of the corrections to the spectator 
widths in Section 4.1. In Section 4.2, the procedure for the estimation of nonperturbative matrix 
elements over the states of doubly heavy baryons is considered for the operators of nonrelativistic 
heavy quarks. Section 4.3 is devoted to the numerical evaluation of lifetimes for the baryons with two 
heavy quarks and their partial decay rates, as well as to the discussion of underlying uncertainties. 
We conclude the chapter by summarizing our results. 

4.1. Operator Product Expansion for the heavy baryons 

In this section we describe the general scheme for the construction of OPE for the total widths of 
baryons containing two heavy quarks. 



^^It was shown in | |7^ that the leading order correction in l/mg is absent and the corrections begin with I/ttiq. 
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Now let us start the description of our approach for the calculation of lifetimes for the doubly 
charmed baryons. The optical theorem, taking into account the integral quark-hadron duality, 
allows us to relate the total decay width of the heavy hadron F with the imaginary part of its 
forward scattering amplitude. This relationship, applied to the H^^-baryon total decay width r=o^, 
can be written down as 

r^.- = ^{Kc\T\Kc). (4.1) 



where the state in (|4.1|) has the ordinary relativistic normalization 



{Kc\Kc) = ^Ev, 

and the transition operator T is determined by the expression 



T 



J d^x {THeffix)Heff 



(0)}, 



(4.2) 



where H^jf is the standard effective hamiltonian describing the low energy weak interactions of 
initial quarks with the decay products. For the transition of c-quark into the w-quark and the 
quarks qi^2 with the charge —1/3, the lagrangian has the form 



H, 



eff 



-:Vu,,v:\cMo^ + C4fi)0^] + h.c. 



where V is the matrix of mixing between the charged currents, and 
SO tat the indices a,P denote the color states of quarks, and 



(4.3) 



a,(M, 



w) 



n/ is the number of flavors. 

Assuming that the energy release in the heavy quark decay is large, we can perform the operator 
product expansion for the transition operator T in ( |4.3| ). In this way we find a series of local 
operators with increasing dimensions over the energy scale, wherein the contributions to F are 
suppressed by the increasing inverse powers of the heavy quark masses. This formalism was already 
applied to calculate the total decay rates for the hadrons, containing a single heavy quark (see 
refs. and |]§,|T^). Here we would like to stress that the expansion, applied in this paper, is 
simultaneously in the powers of the inverse heavy quark mass and the relative velocity of heavy 
quarks inside the hadron. Thus, this fact points to the difference from the description of both 
the heavy-light mesons (the expansion in powers of -^p^) and the heavy-heavy mesons [73|| (the 



expansion in powers of relative velocity of heavy quarks inside the hadron, where one can apply the 
scaling rules of nonrelativistic QCD [p!0|j). 

In this work we will extend this approach to the treatment of baryons containing two heavy 
quarks. The operator product expansion applied has the following form: 



mi mj. 



(4.4) 
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The leading contribution in OPE is determined by the operator cc, corresponding to the spectator 
decays of c-quarks. The use of the equation of motion for the heavy quark fields allows one to 
eliminate some redundant operators, so that no operators of dimension four contribute. There 
is a single operator of dimension five, Qgq = Qd'^fiuG'^'^Q- As we will show below, significant 
contributions come from the operators of dimension six Q2Q2q = QTqqT Q, representing the effects 
of Pauli interference and weak scattering for S^"*" and correspondingly, which are enforced 
by the two-particle phase space. Furthermore, there are also other operators of dimension six 
QeiQ = Q'^nvliD^G^^Q and Qq2Q = QD^G^^VyQ. In what follows, we neglect the corresponding 
contributions for the latter two operators, since they are suppressed by the mentioned smallness of 
three-particle phase space, so that the expansion is certainly complete up to the second order of ^, 
only. The reason for the restriction of dimension-six operators is twofold. First, the operators Qqiq 
and Q&2Q do not contribute to the lifetime difference between the doubly charmed baryons under 
consideration, since they are independent of the quark contents of hadrons. Second, the four-quark 
operators are enhanced in comparison with the quark-gluon terms with the same dimension because 
the two-particle phase space integrated in the calculation of coefficients in front of Pauli interference 
and weak scattering has the additional factor of IGvr^ in contrast to the the three-particle phase 
space expressed in the units of heavy quark mass as it occurs for the coefficients of operators Qqiq 
and Qq2Q, so that they are suppressed. 

Further, the different contributions to OPE are given by the following: 

^++ = %,5c + '^,P/, 
^+ = ^5c + %,WS^ 

where the first terms account for the operators of dimension three O^q and five Ogq-, the second 
terms correspond to the effects of Pauli interference and weak scattering. The explicit formulae for 
these contributions have the following form: 

Tssc = 2(r,,,peccc - %[(2 + iroc)Pi + K2cP2]Ogc), (4.5) 
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where Fqc = and Kq^ = + 2C^, K2C = 2(C^ — Ci). This expression has been derived 

in |78| (see also fT^), and it is also discussed in |76|. The phase space factors Pi [|76|,Q have the 
form 

Pi = (l-#, P2 = il-yf 

2 

where y = Fc^spec denotes the contribution to the total decay width of the free decay for one of 



the two c-quarks, which is explicitly expressed in |81[] 



Quite cumbersome expressions for the contributions of Pauli interference and weak scattering 
terms in the inclusive widths are given in Appendix III. 

In the numerical estimates for the evolution of coefficients (7+ and C_, we have taken into 
account the threshold effects, connected to the 6-quark, as well as the threshold effects, related to 
the c-quark mass in the Pauli interference and weak scattering. 

In the expressions for C± the scale /i is approximately equal to rric. For the Pauli interference and 
weak scattering, this scale in the factor k is chosen in the way to obtain an agreement between the 
experimental differences of lifetimes for the A^, S+ and S[?-baryons and the theoretical predictions 
based on the effects mentioned above. This problem is discussed below. Anyway, the choice of these 
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d, u 

Figure 4.1. The spectator contribution into the total widths of doubly charmed baryons. 



scales allows some variations, and a complete answer to this question requires calculations in the 
next order of perturbative theory. 

The contribution of the leading operator cc corresponds to the imaginary part of the diagram 
in Fig. [4.1| , as it stands in expression ([4.5|) . The coefficient in front of cc can be obtained in 



the usual way by projection of contribution due to the diagram in Fig. to the operator cc. 
This coefficient is equivalent to the free quark decay rate, and it is known in the next-to-leading 
logarithmic approximation of QCD |]82|-p6[], including the strange quark mass effects in the final state 
To calculate the next-to-leading logarithmic effects, the Wilson coefficients in the effective weak 



lagrangian are required at the next-to-leading accuracy, and the single gluon exchange corrections 
to the diagram in Fig. ff]T| must be considered. In our numerical estimates we use the expression 



for Tspec, including the next-to-leading order corrections, s-quark mass effects in the final state, but 
we neglect the Cabibbo-suppressed decay channels for the c-quark. The bulky explicit expression 
for the spectator c-quark decay is placed in the Appendix of |^ . 



Similarly, the contribution by Ogq is obtained, when an external gluon line is attached to the 
inner quark lines in Fig. |4.1| in all possible ways. The corresponding coefficients are known in 
the leading logarithmic approximation. Finally, the dimension six operators and their coefficients 
arise due to those contributions, wherein one of the internal u or d quark line is "cut". The 
resulting graphs are depicted in Figs. |4.2| and [4.3| . These contributions correspond to the effects of 
Pauli interference and weak scattering. We have calculated the expressions for these effects with 
account for both the s-quark mass in the final state and the logarithmic renormalization of effective 
electroweak lagrangian at low energies. This effective lagrangian is shown in Appendix III. 




d 

Figure 4.2. The Pauli interference of c-quark decay products with the valence quark in the initial state for the 
H++-baryon. 



To calculate the contribution of semileptonic modes to the total decay width of S*^-baryons (we 
have taken into account the electron and muon decay modes only) we use the following expressions 
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Figure 4.3. The weak scattering of the valence quarks in the initial state for the H+-baryon. 



9| (see also 



si 



4rc({l - 8p + 8p^ - - 12p2 In p} + 
Ee{5 - 24p + 24p2 _ 8p3 + 3p^ - 12p2 Inp} + 
Kc{-6 + 32p - 24p2 - 2p^ + 24p2 Inp} + 
Gc{-2 + 16p- 16p3 + 2p^ + 24p2lnp}), 



(4.6) 



where Tr 



IT/ 12^2 



= -ft'c + Cc, while Kc and Gc are given by the expressions 



2m? 



0/3 



4m^ 



(4.7) 



where the spinor field Cy in the effective heavy quark theory is defined by the form 



c{x) 



1 + 



2mr 



(4.8) 



where denotes the partial derivative over coordinates. 

The analogous scheme for the calculation of inclusive widths T^o for the baryon S^^, where o 



denotes the electric charge of the system, was developed in ||8^. Then, for the total widths we have 
the expressions 



^bc 
^bc 



%5b + T35C + T^^pi + '^%%s-> 



where two initial terms denote the contributions into the decays of quark Q by the operators with the 
dimensions 3 and 5, and the forthcoming terms are the interference and scattering of constituents. 
Various contributions were explicitly presented in 

The calculation of both the Pauli interference effect for the products of heavy quark decays with 
the quarks in the initial state and the weak scattering of quarks composing the hadron, results in 
summming up the various channels of decays, so that 



-^6,PI 
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'^Q PI — ^PI,sc + ^PI,du + '^P^I,du + ^ '^PI,lPi'> ^"^"^^ 



6,WS ~ -^WS,bu -r -iWS,bc, 



and 



'^PI,du — '^PI,sc (^- 0)' 

^P/,eP. = ^P/,MP„ = ^P/,rP. (^r - 0). (4.10) 



Further, we have taken into account the known ttg-corrections to the semileptonic widths of quarks. 

Thus, the calculation of lifetimes for the baryons with two heavy quarks is reduced to the problem 
of evaluating the matrix elements of operators, that is the subject of next section. 



4.2. Hadronic matrix elements 

According to the equation of motion for the heavy quarks, the matrix element of local operator 
Q^Q^ can be expanded in the following series over the powers of l/niQj, so that 

2M 4Mml^ Mm^Q/' ^ ' ' 

Thus, we have to make numerical estimates for the following list of operators: 

0'7a75<5'?7°(l - 75)?, Q'lalsQ'Q^ri^ - 75)(5^ (4.12) 

Q^7a(i-75)Q^QV(1-75)Q'- 

The first operator corresponds to the time dilation connected to the motion of heavy quarks inside 
the hadron. The second is related to the spin interaction of heavy quarks with the chromo-magnetic 
field of light quark and the other heavy quark. The further operators are the four-quark operators 
representing the effects of Pauli interference and weak scattering. 

Next, following the general methods of effective theories, we introduce the effective filed ^q, 
which represents the nonrelativistic spinor of heavy quark, so that we integrate out the contributions 
with the virtualities /i in the range mg > n > mqvg in the framework of perturbative QCD, 
while the nonperturbative effects in the matrix elements could be expressed in terms of effective 
nonrelativistic fields. So, we have 

QQ = *^*Q-^vI;^(iD)^*c? + ^*^(zD)^*Q- 

^^^^,gaB^,-^^^liDgm, + ... (4.13) 

Qga.^G^-'Q = -2*^^o-S^Q-— ^^(£>5i;)^Q + ... (4.14) 
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In these expressions we have omitted the term "^gcr^gE x Z))\1/q, corresponding to the spin-orbital 
interaction, because it vanishes in the ground states of baryons under consideration. For the nor- 
mahzation we take 

j ci^x^^^Q = j £xQ^Q. (4.15) 

Then for Q determined by 

Q = e-'™* (^^^ (4.16) 

we have 

Let us emphasize the difference in the descriptions of interactions between the heavy quarks 
and hght quark and interactions between two heavy quarks. As we have mentioned, there is an 
additional small parameter in the heavy subsystem. It is the relative velocity of quark motion, that 
introduces the energetic scale rnqv. Therefore, for example, the Darvin term (DE) in the heavy 
subsystem turns out to be the same order of magnitude in comparison with the chromo-magnetic 
term {crB) in the expansion over the inverse power of heavy quark mass, since it has the same 
degree of smallness in v. This fact becomes clear if we explore the scaling relations derived in the 



nonrelativistic QCD |]10 



1 



"^Q^irriQVQ)^, \D\^mQVQ, g\E\ r-. uiqVq, g\B\ ^ mgVQ, c/ ~ v^. 

For the interaction of heavy quark with the light one there is no small parameter defined by the 
relative velocity, so that the Darvin term is suppressed by the factor of k/mq ~ KqcdIi^q- 



4.2.1. Baryons S+ and 

Let us now start the calculation of matrix elements with the use of potential models for the bound 
states of hadrons. While estimating the matrix element value of the kinetic energy, we note, that 
the heavy quark kinetic energy consists of two parts: the kinetic energy of the heavy quark motion 
inside the diquark and the kinetic energy related to the diquark motion inside the hadron. According 
to the phenomenology of meson potential models, in the range of average distances between the 
quarks: 0.1 fm < r < 1 fm, the average kinetic energy of quarks is constant and independent of 
both the quark flavors, composing the meson, and the quantum numbers, describing the excitations 
of the ground state. Therefore, we determine T = mdv'^/2 + mivf /2 as the average kinetic energy 
of diquark and hght quark with nii = niq, and T/2 = mciv'^i/2 + mc2fc2/2 as the average kinetic 
energy of heavy quarks inside the diquark (the coefficient 1/2 takes into account the anti-symmetry 
of color wave function for the diquark). Finally, we have the following expression for the matrix 
element of the heavy quark kinetic energy: 

2M=o^m^ 2m1 + mcmi 2 

We use the value T ~ 0.4 GeV, which results in = 0.146, where the dominant contribution comes 
from the motion of heavy quarks inside the diquark. 

Now we would like to estimate the matrix element of chromo-magnetic operator, corresponding 
to the interaction of heavy quarks with the chromo-magnetic field of the light quark. For this 
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purpose, we will use the following definitions: Omag = Yl'i=i -^(^ tivG^^ c!' and Omag ~ + 1) — 
Sdi^d + 1) ~ ^li^i + 1))) where is the diquark spin (as was noticed by the authors earlier, there 
is only the vector state of the cc-diquark in the ground state of such baryons), Si is the light quark 
spin and j is the total spin of the baryon. Since both c-quarks additivcly contribute to the total 
decay width of baryons, we can use the diquark picture and substitute for the sum of c-quark spins 
the diquark spin. This leads to the parametrization for Omag, as it is given above, and, moreover, it 
allows us to relate the value of the matrix element for this operator to the mass difference between 
the excited and ground state of baryons 

Omag = "^(Ms^o - M^oJ. (4.19) 

The account for the interaction of heavy quarks inside the diquark leads to the following expressions 
for the chromo-magnetic and Darwin terms: 

2M=o "9^ me ' ^ 

'— 'CC ^ 

(iyii(£_Ma.2^.|*.(0)p, (4.21, 

■— 'CC 

where ^'^(O) is the diquark wave function at the origin. Collecting the results given above, we find 
the matrix elements of dominant operators determining the spectator decay 



^\0)\' - —rg'\^\0)\' + ... 



^cc ^ C C 

^ 1 - 0.074- 0.004- 0.003- 0.005 (4.22) 

We can see that the largest contribution to the decrease of the decay width comes from the time 
dilation, connected to the motion of heavy quarks inside the baryon. For the matrix element of the 
operator cgGnj^G^^c, we get 

{^l.\c9<^,.G^'c\^oj ^ _4 (M.,, - M.,J _ V ,^..Q.|2 _ ^|^.(o)|2 .4 23) 

2M^.ml 3 9m3 ' ^ '^ 3m3 ' ^ ^' ' ^ > 

Now let us continue with the calculation of the matrix elements for the four-quark operators, 
corresponding to the effects of Pauli interference and weak scattering. The straightforward calcu- 
lation in the framework of nonrelativistic QCD gives the following expressions for the operators 
(cTm(1 - lb)c ■ q^y'^il - 75)9) and (07^750 • ^7^(1 - 75)?): 

(C7^(l-75)c-g7ni-75)g) = 2m,V-\l-4S,Sg), (4.24) 
(c7M75C-g7'^(l-75)g) = -iS,Sg-2m,V-\ (4.25) 

where V"^ — \^i{0)\'^, is the hght quark wave function at the origin of two c-quarks. We 

suppose that, for estimates, |^z(0)| has the same value as that in the D-meson. So, we find 

l*KO)r«^. (4.26) 
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We would like to note that in the numerical estimate the supposed value of light quark wave function 
obtained from the leptonic constant of charmed meson fo 200 MeV, is about two times less than 
the wave function calculated in Chapter 1 in the approximation of quark-diquark factorization. 
This is related with the fact that the leptonic constants of charmed hadrons have large corrections 
of both logarithmic form and powers of inverse charmed quark mass. So, the leptonic constant of 
D meson calculated in the potential models is two times greater than the value obtained in the 
QCD sum rules with account for the mentioned corrections. Therefore, we think that the supposed 
approximation for the wave function of light quark in the baryon is justified enough. 

Then, again remembering that both c-quarks additively contribute to the total decay width and 
using the diquark picture, we can substitute for Sc^ + Sc^ by Sd, where Sd is the diquark spin. Thus, 
for the matrix elements of operators under study we have the following expressions: 

(H:,|(c7^(1 - j,)c){qr{l - 75)g)|5L) = 12^^ • I^^O)!', (4.27) 

(S:j(c7M75c)(g7''(l - 75)g)|S:,) = 8m, ■ |vl/,(0)|l (4.28) 

The color anti-symmetry of the baryon wave function results in relations between the matrix ele- 
ments of operators with the different sums over the color indices: (S*^|(ciT^Cfc)(gfc7''(l — 75)5,12*^) = 
-(S^^|(c7)(C)(^7^(l - 75)g'|S*^), where is any spinor structure. 

Thus, we have formally constructed the procedure for the evaluation of matrix elements obtained 
in the OPE of T for the baryon with two identical heavy quarks. 



4.2.2. Baryons r.'^^ and 

Considering the baryons with two heavy quarks of different flavors, we point out the modiflcations, 
which have to be introduced in the estimation of hadronic matrix elements for the quark operators 
determining the inclusive decay widths. 

First of all, according to the quar-diquark factorization we get the following expressions for the 
kinetic terms: 

2Mwo^ml {mi + rrib + mc){mb + rric) mc{mc + mf,) 

2Ms«^m^ ^ (mi + nib + mcjirrih + rric) ^,5(772, + mt) 

Numerically, we assign T ~ 0.4 GeV, which results in v I = 0.195 and = 0.024, where the 
dominant contribution is provided by the motion inside the diquark. 
Let us define 

Omag = -^ca^.G^''' c + ^ba^^C^^b, (4.31) 

(Omag) = —{Sa{Sa + l)-S,{Sc+l)-SiiSi + l))+ (4.32) 
rric 

—{Su{Su + 1) - SbiSb + 1) - Si{Si + 1)), 
nib 
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where Su = Sb + Si, Sd = Sc + Si , Sb is the 6-quark spin, Sc is that of c-quark, and 5*^ is the hght 
quark spin. The operator under study is related to the hyperfine sphtting in the baryon system 

3 3 1 1 

{Sbc = 1,5'= —\Omag\Sbc ~ '^^ S = —) — {Sbc =1,5 = —\Omag\Sbc ~ 1) 5 = — ) 

= {Sbc =l,S = ^\VHf\Sbc=l,S = ^)- {Sbc =l,S = ^\VHf\Sbc = l,S = ^), (4.33) 

so that S denotes the total spin of the system, and Sbc is the diquark spin. Further, the perturbative 
term, depending on the spins, is equal to 

Vhj = la^^SiSclR'^m' + ^a.^SiSblR'^m', (4-34) 
9 minic 9 minib 

where W^^O) is the radial wave function at the origin of quark-diquark system. In contrast to the 
diquark system with the identical quarks, this operator is not diagonal in the basis of S and Sbc- 
To proceed, we use the change of basises 

\S- Sbc) = 5^(-l)(^+^^+^^+^^) v/(25,, + 1)(25,, + 1) I I I \S- Sbi) (4.35) 

and 

\S; Sbc) = ^(-i){s+^^+^^+^^)^(25,, + 1)(25,, + 1) I I I 1^ I \S; Sd) (4.36) 
The result of substitutions gives 

A = &>^. (4.37) 

9m, ' ^ ^ 

however, for the state with the zero spin of heavy diquark, which is considered in what follows, we 
have 

\^bc\^mag\^bc/ _ g ooN 

"be 

The account of Darwin and chromo-magnetic terms results in 

{EimgcT-B^cl^tc) ^ _2 .4 39^ 

^bc " 

{ZlM{D-gE)^>c\'^i:, 2 ,,^,,2 



2M= 



^'|^''(0)|^ (4.40) 



where \E'''(0) is the wave function at the origin of diquark. The analogous matrix elements^ for the 
operators of beauty quarks can be written down by the permutation of heavy quark masses. 
Combining the results above, we find 

= l-k + Tr^l*'(0)P-^/|^^(0)P + ... 

^ 1 - 0.097 + 0.004- 0.007... (4.41) 



-•^^The obtained expressions differ from those for the Be meson ||7J] because of the color structure of the state, 
providing the factor of i . 
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The dominant role in the corrections is played by the term, connected to the time dilation because 
of the quark motion inside the diquark. Next, for the operator cga^uG^'^c we have 



{El\cga,,G'^''c\El) 
2M=o m2 

^bc 



9^ 



3ml 



|^"'(0)|^ ^ 0.002. 



(4.42) 



The permutation of quark masses lead to the required expressions for the operators of bb and 
bga^uGf^'^b. 

Making use of ( [4.35|) and (|4.36|) for the ground states of baryons results in 

(4.43) 
(4.44) 
(4.45) 
(4.46) 
(4.47) 
(4.48) 

Thus, we can make numerucal estimates of inclusive decay widths. 



(S:J (67,(1 -75)&)(C7"(1 - 


- 75)c) 


^bc) = 


= 8(mfe 


+ 


rric) 






- 75)c) 


^bc) = 


= 6(mfe 


+ 


rric) 


■ |^''(0)r, 




- 75)9) 


^bc) ~ 


= 2(mfc 


+ 


m) 


|^'^'(o)p. 


(s:j(67,75&)(gV(i- 


- 75)9) 


^bc) = 


= 0, 








(S^J(c7,(l-75)c)(g7^(l- 


- 75)g) 


^bc) ~ 


= 2(mc 


+ 


mi) 


|^'^^(o)p. 




- 75)?) 


^bc) = 


= 0. 









4.3. Numerical estimates 

Now we are ready to collect the contributions, described above, and to estimate the lifetimes of 
baryons and H+. For the beginning, we list the values of parameters, which we have used in 
our calculations, and give some comments on their choice. 

m^ = 1.6 GeV, = 0.45 GeV, |Ks| = 0.9745, 

M^++ = 3.56 GeV, M=+ = 3.56 GeV, Ms- - M^o = 0.1 GeV, (4.49) 
T = 0.4 GeV, 1^(0)1 = 0.17 GeV^, = 0.30 GeV. 



For the parameters Mc++ , and M=.o — M=o we use the mean values, given in the literature. 
Their evaluation has been also performed by the authors in the potential model for the doubly 
charmed baryons with the Buchmiiller-Tye potential and also in refs. For fo 



we use the value, given in and for T we take it from |^8[. The mass mc corresponds to the 



pole mass of the c-quark. For its determination we have used a fit of theoretical predictions for 
the lifetimes and semileptonic width of the D°-meson from the experimental data. This choice of 
c-quark mass seems to effectively include unknown contributions of higher orders in perturbative 
QCD to the total decay width of baryons under consideration. 

The renormalization scale /z is chosen in the following way: /ii = in the estimate of Wilson 
coefficients C for the effective four-fermion weak lagrangian with the c-quarks at low energies and 
II2 = 1-2 GeV for the Pauli interference and weak scattering. The latter value of renormalization 
scale has been obtained from the fit of theoretical predictions for the lifetimes differences of baryons 
Ac, H^, H° over the experimental data. Here we would like to note, that the theoretical approxi- 
mations used in include the effect of logarithmic renormalization and do not take into account 
the mass effects, related to the s-quark in the final state. For the corresponding contributions to 
the decay widths of baryons with the different quark contents we have 

Ar„;(Ac) = Cd{Od)k, + Cu{Ou)Ac, 

Ar„,(S+) = c,(0,)s++c„(0„)s+, (4.50) 
Ar„KS°) = Cd{Od)E0 + Cs{Os)E0, 
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where {Oq)xc = {^c\Oq\Xc), Oq = {c~^fic){q'-y^q) with q = u,d or s. The coefficients Cq{fi) are equal 
to 



cu = -^^[(C+ + C_)^ + ^(1-4A;^)(5C^ + C2 -6C+C_)], (4.51) 

cs = + i(l-4fc^)(5CUC^+6C+C_)]. 

lovr 6 

We use the spin averaged value of the D-meson mass for the estimation of the effective hght quark 
mass = A as it stands below 

mD = mc + A + -^ = mc + mi + = 1.98 GeV. (4.52) 

znic Trie + TUi 

The s-quark mass can be written down in terms of mi introduced above, so that 

= + 0.15 GeV. (4.53) 
As we have already mentioned, the spectator decay width of c-quark Tc,spec is known in the next- 



to-leading order of perturbative QCD [p4-|86|. The most complete calculation, including the mass 
effects, connected to the s-quark in the final state, is given in [^. In the present work we have used 
these results for the calculation of the spectator contribution to the total decay width of doubly 
charmed baryons. In the calculation of the semileptonic decay width, we neglect the electron and 
muon masses in the final state. Moreover, we neglect the r-lepton mode. 

Now, let us proceed with the numerical analysis of contributions by the different decay modes 



into the total decay width. In Table |4.1| we have listed the results for the fixed values of parameters, 
described above. From this table one can see the significance of effects caused by both the Pauli 
interference and the weak scattering in the decays of doubly charmed baryons. The Pauli interference 
gives the negative correction about 63% for the S+^-baryons, and the weak scattering increases the 
total width by 61% for S+. As we have already noted above, these effects take place separately in 
these baryons, and, thus, they enhance the difference of lifetimes. 



Table 4.1. The contributions of different modes to the total decay width of doubly charmed baryons. 



Mode or decay 
mechanism 


Width, ps ^ 


Contribution, % 


Contribution, % 


c sdu 


2.648 


127 


31 


c se'^u 


0.380 


18 


4.2 


PI 


-1.317 


-63 




WS 


5.254 




60.6 




2.089 


100 






8.660 




100 



It is worth here to recall that the lifetime difference of and D -mesons is generally explained 
by the Pauli interference of c-quark decay products with the anti-quark in the initial state, while 
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in the current consideration, we see the dominant contribution of weak scattering. This could not 
be surprising, because under a more detailed consideration we will find, that the formula for the 
Pauli interference operator for the Z)-meson coincides with that for the weak scattering in the case 
of baryons, containing, at least, a single c-quark. 

Finally, collecting the different contributions for the total lifetimes of doubly charmed baryons, 
we obtain the following values: 

r^++ = 0.48 ps, r=+ = 0.12 ps. 

Note that the supposed exploration of fitting the data on the semileptonic decay width of D mesons, 
the difference between the widths of baryons containing the charmed quark as well the spectroscopic 
characteristics allow us to significantly decrease the variations of model parameters, i.e. the quark 
masses, the scale for the normalization of Wilson coefficients, and the wave function of light quark in 
the nonrelativistic model. In this way we decrease the uncertainties of theoretical predictions. The 
variation of charmed quark mass in the limits of 1.6 — 1.65 GeV and difference between the masses of 
strange quark and light quark in the range 0.15 — 0.2 GeV results in the following uncertainties for 
the lifetimes of baryons under consideration: 6t^++ = 0.1 ps, (5r=+ = 0.01 ps, whereas we see that 
the abolute values of uncertainties in the widths are equal to (5rc++ = 0.4 ps~^, (5r^+ = 0.9 ps~^. 

J. <— 'CC ^ ' ^cc 

However, since the total width of S+ is significantly enhanced by the contribution of weak scattering 
for the constituent quarks, the relative uncertainty in the estimate of lifetimes for this baryon is 
much less, i.e. it is 10% in comparison with 20% for 

Calculating the inclusive widths of decays for the and we have used 

mb = m^ + 3.5 GeV, (4.54) 

in addition to the parameters given above. The baryon mass has been put to 7 GeV. For the wave 
function in the subsystem of diquark we have used the results of calculations in the nonrelativistic 
model with the potential by Buchmiiller and Tye [T^\ corrected for the color structure of diquark, 
so that 

^'^(0) = 0.193 GeV^/^ 

Further, it is quite evident that the estimates of spectator widths of free heavy quarks do not 
depend on the system, wherein they are bound, so that we can exploit the results of the calculations 
performed earlier. We have chosen the quark masses to be the same as in |]7^ , and we have put the 
corresponding values, presented in Table 



Table 4.2. The widths of inclusive spectator decays for the b and c quarks, in ps 



mode 


b cud 


b — > CCS 


b — > ce+z/ 


r 


0.310 


0.137 


0.075 


mode 


b CT^V 


c — s> sdu 


c — ^ se'v 


r 


0.018 


0.905 


0.162 



Then the procedure, described above with the shown parameters, leads to the lifetimes of the 
1^ and baryons equal to 

= 0.33 ps, r^o^ = 0.28 ps. (4.55) 
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We can clearly see that the difference in the lifetimes caused by the decay processes with the Pauli 
interference and weak scattering is about 20 %. The relative contributions by various terms in the 
total widths of the baryons under consideration are presented in Table |4.3| . 

Table 4.3. The branching fractions (in %) of various modes in the inclusive decays of and baryons. 



mode 


Tb 


Tc 


Tpi 






20 


37 


23 


20 


"6c 


17 


31 


21 


31 



Note that the contributions by the Pauli interference and weak scattering depending on the 
baryon contents, can be quite significant, up to 40 — 50%. The corrections due to the quark-gluon 
operators of dimension 5 are numerically very small. The most important corrections are those due 
to the operator of dimension 3, where the role of time dilation is essential for the heavy quarks in 
the hadron rest frame. 

For the semileptonic decays, whose relative fractions are presented in Table |4.4| , the largest cor- 
rections appear in the decays of 6-quark because of the Pauli interference, so that the corresponding 
widths practically increase twice. This leads to the result that the semileptonic widths of b and c 
quarks in the electron mode are equal to each other, whereas for the spectator decays, the width of 
the charmed quark is twice that of b. 



Table 4.4. The branching ratios for the inclusive semileptonic widths of Et and E^^, in %. 



mode 




b 


^ b 




5.0 


4.9 


2.3 


^bc 


4.2 


4.1 


1.9 



As for the sign of terms, caused by the Pauli interference, it is basically determined by the 
leading contribution, coming from the interference for the charmed quark of the initial state with 
the charmed quark from the 6-quark decay. In this way, the anti-symmetric color structure of baryon 
wave function leads to the positive sign for the Pauli interference. 

Finally, concerning the uncertainties of the estimates presented, we note that they are mainly 
related to the following: 

1) the spectator width of charmed quark, where the error can reach 50 %, reflecting the agreement 
of theoretical evaluation with the lifetimes of charmed hadrons, so that for the baryons under 
consideration this term produces an uncertainty of 6r /T ^ 10%, 

2) the effects of Pauli interference in the decays of beauty quark and in its weak scattering off 
the charmed quark from the initial state, wherein we use the nonrelativistic wave function, which, 
being model-dependent, can lead to an error estimated close to 30%, producing an uncertainty of 
ST/T ^ 15% in the total widths. 

Thus, we estimate, that the uncertainty in the predictions of total widths for the and 
baryons is about 20%. 
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4.3.1. Parametric dependence of results 

Despite of tiiat we have already pointed out what is the expected accuracy of predictions on the 
inclusive widths and lifetimes of baryons with two heavy quarks, let us consider this important 
problem in a more detail way. 

First of all, we have to investigate the dependence of widths on the masses of heavy quarks 
composing the diquarks in the baryons. Indeed, the spectator decay widths of heavy quarks are 
determined by the fifth power of masses, while the dominant corrections due to the Pauli interference 
and weak scattering of constituents depend as the third degree of the heavy quark masses. In this 
way, a natural challenge is the justification of quark-hadron duality with the OPE of quark operators, 
that leads to zero term of correction in the first inverse power of heavy quark mass, I/uiq. Such 
the nullification follows from the theorem by Adomolo-Gato claiming that the introduction of term, 
which is proportional to A and breaking a symmetry of lagrangian, leads to corrections of second 
order, i.e. A^, in observable quantities, so that under the introduction of I/ttt-q in the lagrangian 
of heavy quarks, the decay widths of heavy hadrons do not contain contributions linear in l/rriQ, 
if we suppose the quark-hadron duality. In the connection with the problem on the lifetime of 
baryon, whose total width measured experimentally is about 20% greater than the widths of B 
mesons, that is in a contradiction with the predictions of heavy quark theory ||90|, authors of |[91[ 



suggested the hypothesis on a strong violation of quark-hadron duality, i.e. on the possibility of 
significant contribution by terms linear in 1/ mq into the inclusive decay widths of heavy hadrons. 
This assumption, in fact, implies that the heavy quark effective mass determining the contribution 
of leading term, can vary with the mass and contents of hadron. So, considering a more large 
system of A;, baryon with the light diquark, wherein the string tension is twice less than the tension 
in the meson, we have to introduce a larger effective mass of heavy quark, since it is determined at a 
smaller energetic scaleQ. Then, the total width of A;, increases. Such the approach is not acceptable 
in the OPE with the quark-hadron duality, which operates with the heavy quark mass being the 
same for all kinds of hadrons, otherwise the corrections linear in l/mq could appear in the widths. 
However, the hypothesis of the strong violation of quark-hadron duality is practically removed 
by the experimental measurement of meson lifetime. Indeed, the author of |Q following the 



ideology of predicted the value of lifetime t[Bc] ~ 1.3 — 1.5 ps, since the quarks inside the heavy 



quarkonium are strongly bound, and their effective masses as well as the admissible phase space of 
the final state in decays decrease. This fact leads to a significant suppression of decay widths for 
both b and c quarks in the Be meson. The experimental measurement yields t[Bc] = 0.48 ± 0.19 
ps, which is in a good agreement with the estimates in the framework of OPE |7^,^, QCD sum 
rules [0 and potential models Thus, at present we can claim that the OPE with the quark- 
hadron duality is a correct tool in the calculations of inclusive decay widths and lifetimes for the 
hadrons with heavy quarks. 

As we have mentioned, the estimates presented in the above sections have been done under the 
assumption that the chosen value of charmed quark mass results in quite the precise agreement of 
theoretical prediction on the semileptonic widths of D mesons with the experimental values. Behind 
our paper |8l[| the analogous calculations were performed in ref. [^], whose authors supposed 



essentially lower value for the charmed quark mass iric = 1.35 GeV, which certainly resulted in 
that the semileptonic widths of charmed mesons could not be satisfactorily described in the same 
approach. Such the preference was probably caused by the following: first, a low value of current 
mass of c quark is usually obtained in the QCD sum rules for the charmonium, and, second, the 
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The cloud of virtual gluons and quarks has a larger size in the heavy baryon than in the heavy-light meson. 
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description of inclusive decay widths of D mesons is generally considered to be quite qualitative, 
but quantitative, because the charmed quark mass is not large enough, therefore, the convergency 
of series in l/niq could be slowQ. Such the suggestions do lead to just a quahtative predictions 
on the lifetimes of baryons with two heavy quarks in |95], so, their estimates are twice or thrice 
different from the values obtained in previous sections. Indeed, the significant decrease of leading 
term because of the variation of charmed quark mass leads to that the negative contribution of 
Pauli interference strongly diminishes the total widths of so that the lifetime is essentially 

overestimated. 

Another important source of uncertainty in the estimates of widths for the decays of charmed 
quark is the mass of strange quark. Indeed, since the charmed quark has the mass about 1.5 GeV, 
the phase space of its decay strongly depends on the chosen mass of strange quark: the current mass 
about 150 — 200 MeV or the constituent mass close to the mass of K meson. In previous estimates 
we have supposed that the suppression of phase space is determined by the constituent mass. The 
problem on the dependence of inclusive decay widths on the masses of heavy and strange quarks 



was studied in detail in ref. |]96| under the consideration of Be meson lifetime, where the uncertainty 
caused by the modelling the wave function is lowQ and, moreover, the corresponding contribution 
of weak annihilation (corrections in the second order of I/tuq) is small (about 10%). 

Accepting ( [4.54| ), that follows from the analysis of data on the decays of B mesons, authors 
of [|96| get the estimates presented in Table [4.5| . In this table we see that, first, the low value of 
charmed quark mass supposed in [^| gives quite an overestimating value for the lifetime of -Be- 
Second, the choice of current mass for the strange quark is slightly preferable, since it leads to the 
value of Be lifetime, which well agrees with the central value of experimental inteval, though the 
uncertainty of data allows the description with the constituent mass of strange quark, too. It is 
worth to note that this analysis yields the quark mass values consistent with the choice of charmed 
quark mass under the semileptonic decay widths of D mesons. 



Table 4.5. The lifetime of Be meson and contributions of spectator widths as well as of corrections caused by 
the Pauli interference (PI) and weak annihilation (WA) at different values of quark masses. 



Parameters, GeV 


b ^ c, ps ^ 


c — s, ps ^ 


PI, ps-i 


WA, ps-i 


tb,, ps 


rrtfe = 5.0, rric = 1.5, = 0.20 


0.694 


1.148 


-0.115 


0.193 


0.54 


ruh = 4.8, rric = 1-35, = 0.15 


0.576 


0.725 


-0.132 


0.168 


0.75 


nib = 5.1, rric = 1-6, = 0.45 


0.635 


1.033 


-0.101 


0.210 


0.55 


rrib = 5.1, rric = 1-6, rris = 0.20 


0.626 


1.605 


-0.101 


0.210 


0.43 


nib = 5.05, mc = 1.55, = 0.20 


0.623 


1.323 


-0.107 


0.201 


0.48 


rrib = 5.0, TTT-c = 1.5, ms = 0.15 


0.620 


1.204 


-0.114 


0.193 


0.53 



Along with the quark masses, which, to the moment, could be generally considered as quantities 
with quite low uncertainties, the variation of light quark wave function plays a significant role in the 
calculations of inclusive decay widths for the baryons with two heavy quarks. As we have explained 
above, this quantity was estimated under the assumption that it is analogous to the wave functions 
of D mesons, i.e. we suggest that the corrections to the value obtained in the framework of potential 

Again, one does not take into account the data on Be- 
^^The heavy quarkonium is quite precisely described due to lots of data on the charnionium and bottomonium. 
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models are similar in the mesons and baryons. In the analysis of |9^,^ the following relation for 
the wave function of light quark in the doubly heavy baryon was supposed: 



4 



1^,(0)1= ^ (4.56) 



4 



where = 170 MeV, and the factor of A;~9 is caused by the so-called "hybrid" logarithms for the 
nonrelativistic heavy quarks. Such the expression is derived if we suppose the scaling of hyperfine 
spin-spin splitting in the charmed mesons and baryons and take into account corresponding spin 
factors as well as the double mass of diquark composed by two heavy quarks. The assumption 
on the independence of splitting on the mesonic or baryonic state looks quite speculative, while, 
nevertheless, if we leave physical motivations, the numerical effect is the reduction of wave function 
of light quark by two or three times. In contrast, the calculations in the framework of potential 
models result in the twice enhancement of wave function factor. Thus, the numerical value accepted 
in previous sections gives the central value for the widths under the variation of wave function of 
light quark. 



Table 4.6. The lifetime and inclusive widths of baryon. 



Parameters, GeV 


Y,c^ s, ps-i 


PI, ps-i 


7"=++, ps 


rric = 1.35, = 0.15 


1.638 


-0.616 


0.99 


TUc = 1.6, TUs = 0.45 


2.397 


-0.560 


0.56 


rric = 1.55, = 0.2 


3.104 


-0.874 


0.45 



Table 4.7. The lifetime and inclusive widths of baryon. 



Parameters, GeV 


x;c ^ s, ps ^ 


WS, ps-i 


r=+ , ps 


rric = 1-35, = 0.15 


1.638 


1.297 


0.34 


rric = 1-6, rrig = 0.45 


2.397 


2.563 


0.20 


rric = 1-55, = 0.2 


3.104 


1.776 


0.20 



Table 4.8. The lifetime and inclusive widths of 17+ baryon. 



Parameters, GeV 




PI, ps-i 




rric = 1.35, = 0.15 


1.638 


1.780 


0.30 


rric = 1-6, rris = 0.45 


2.397 


0.506 


0.34 


rUc = 1.55, rris = 0.2 


3.104 


1.077 


0.24 



For the sake of presentation on the degree of variations in the theoretical predictions for the 
inclusive widths of doubly heavy baryons we show the estimates from |9^ exploring the underesti- 



mated value of wave function of light quark, in Tables [4.6[-[4.8|, which should be compared with the 



results in Table 4.1 



92 



Remember that the estimates with the low value of charmed quark mass give the illustrations, 
only, and they cannot be correct because of the contradiction with the data on the lifetime of B^. 

Summimg up the analyzed uncertainties caused by the masses and wave function of light quark 
in the baryon, we present our final estimates in Table 4.9. 



Table 4.9. The lifetimes of doubly heavy baryons. 



baryon 


r, ps 


baryon 


T, ps 


baryon 


r, ps 


■=:++ 

' — 'cc 


0.46 ±0.05 


"be 


0.30 ±0.04 


"66 


0.79 ±0.02 


' — 'cc 


0.16 ±0.05 


^0 

"be 


0.27 ±0.03 


"66 


0.80 ±0.02 


^tc 


0.27 ±0.06 




0.22 ±0.04 


^6^6 


0.80 ±0.02 



In ref. |93l the analysis comparing the structures of OPE for the heavy hadrons was done on 
the basis of symmetry in the hadronic matrix elements determining the contributions by the Pauli 
interference and weak scattering of constituents^, so that up to both corrections in the inverse 
powers of heavy quark mass and logarithmic terms, which are given by anomalous dimensions of 
corresponding operators, the following scaling relations were derived: 



r[B-] - _ r[E-] - r[5°] _ r[E,,] - r[H°,] _ \v^,\ 
T[D+] - r[DO] T[Et] - r[HO] r[s++] - r[s+] 



(4.57) 



The accuracy of such relations should be of the order of 50%, since, for instance, according to the 
consideration of leptonic constants for the heavy mesons, the hadronic matrix elements of quark 
currents with the charmed and light quarks acquire large corrections about 50-90% certainly due 
to both the 1 / mq terms and the logarithmic renormalization. Since the cc-diquark is twice heavier 
than the charmed quark, in the leading order we can put that the mentioned corrections could be 
twice smaller for the doubly heavy baryons. We have explored the data of Table [4.9| in order to 
test the last equality in ( [4.571 ). First of all, as we see, the accuracy of theoretical predictions does 
not allow us to make some convincing quantitative conclusions on the difference of lifetimes for 
the baryons with two b quarks. If we restrict ourselves by the difference of central values, then 
the studied part of relation ( [4.57|) is really satisfied with the accuarcy of 50%, which points to a 



qualitative consistency of relation, whose quantitative accuracy is sadly low. 
4.4. Exclusive decays in sum rules of NRQCD 

In this section we describe the calculations of exclusive semileptonic cascade decays of doubly heavy 
baryons as well as two-paricle hadronic decays in the approximation of factorization for the weak 
trasition current of quarks [p3]. 



In the framework of NRQCD sum rules, the following form of baryonic current was chosen in |^ 



Jho^ = 5"^^ : {QlC^,qp)Q, :, (4.58) 

which leads to the necessary anti-symmetrization in corresponding diagrams, since two identical 
heavy quarks can enter the baryon. In this way, while taking the matrix element of chosen current 



^•^In this way, we deal with the wave function of Ught quark independent on the flavor of infinitely heavy quark 
being the source of gluon field determining the motion of light quark. 
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over the physical baryonic states and the vacuum there is a component giving zero unphysical 
contribution, which is not essential In this approach, the baryonic couphng constants are generally 
different from the values calculated above. Therefore, we have to analyze the two-point correlation 
functions with the new choice of current structure, that was done in However, the new choice 
has a formal advantage in the consideration of three-point correlators determining, for instance, the 
formfactors of semileptonic decays (see Fig. 



Pi 




Pf 



Figure 4.4. The quark loop for the thre-point correlator in the decay of baryon. 



Consider the correlator 



= jd^xdMo\T{JHAy)Uo)JHAx)}\o)e'p^-''e-'p'-y, (4.59) 

where is the weak current of quark decay. The theoretical part of sum rules can be expressed in 
the form of dispersion relation 

where dots denote possible subtractions providing the convergency of integrals, and the spectral 



densities were calculated in in the limit of spin symmetry in the effective lagrangian of heavy 
nonrelativistic quarks under account of quark loop and corrections due to the condensate of light 
quarks. In the approximation of symmetry the only scalar correlator should be calculated. Indeed, 
the hadronic part od sum rules has the form 

nr-'fe.^..^^) = E '"'"^Jm'"'" {h.MMH,M} (4.61) 

spins 2 Hp 1 Hj 

where the formfactor of decay for the baryon with the spin | into the baryon with the spin | is 
expressed in the general form 

{Hf{pf)\J,\Hj{pj)) = u{pp) H,Gr + + v^G^ + 

l^il.Gf + v^Gi + vlGi)} u{pi). (4.62) 
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In general, all of six formfactors are independent. However, in the leading order the lagrangian 
of NQRCD possesses the symmetry, so that at small recoil momenta restricting the virtualities of 
gluonic exchanges in the hadronic state we can derive relations connecting the formfactors with 
each other, if they give nonzero contributions. So, in this limit the 4-velocities of baryons in the 
initial and final states slightly deviate from each other vj 7^ vp-, while their scalar product is close 
to unit w = {vi ■ vp) — ^ 1. Then, the correlation function for the decay of heavy quark into the 
heavy quark has the form 



where 



Vi 

vf 



lull 
^3 

Vf + - — [vf - Vi). 
2m2 



(4.63) 

(4.64) 
(4.65) 



We see that the correlation function is determined by the only formfactor at the minimal recoil 
momentum. However, this formfactor is not universal, i.e. it depends on the quark contents of 
baryons in the process of decay. 




15 5^, GeV' 



Bi, GeV 



Figure 4.5. The forfactor ^(1) for the transition of —>■ SL in the Borel schenne of sum rules 



For the decay of heavy quark into the light one we find 
wherefrom we can get the spin-symmtry relation 



(Cr + + G^) 



and the following connection between the functions C,i{w) takes place: 



(4.66) 



(4.67) 
(4.68) 
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Exploring the presence of heavy hadrons in both the initial and final states, we draw the conclusion 
that two form factors GX = = C,{w) are not suppressed by the heavy quark mass, only. 
The formfactor ^ at nozero recoils is determined in accordance with 



(27r)2 SMiMpZiZp J^mi+msy i(mi+m2)2 ' 

, si-Mh , sf-MI, 
X exp(-^-^) exp(-^-^), (4.69) 



where Bi and Bp are the parameters of Borel transform over the invariant squares of masses in the 
initial and final states of decay, while Zi and Zp are the coupling constants of baryons with the 
appropriate currents of quarks. 

4.4.1. Numerical estimates 

The estimates obtained in the sum rules for the function analogous to the Isgur-Wise function are 



presented in Table [4.10| for the formfactors of semileptonic decays of doubly heavy baryons into the 
baryonic states with the spin |, in comparison with the values calculated in the potential model. 
The deviation between the values of ^(1) in these two approaches does not lead to any discrepancy, 
since the corresponding systematic uncertainty is about 10%. 



Table 4.10. The normalization of Isgur-Wise formfactor at zero recoil in the sum rules and potential model. 



Mode 


^(1), sum rules 


^(1), pot. model 


^bb ^bc 


0.85 


0.91 


^bc ^ ' — 'CC 


0.91 


0.99 


^bc ^ ' — 'bs 


0.9 


0.99 


^cc ^ ' — 'CS 


0.99 


1. 



The result on the normaliztion of Isgur-Wise function in the transition of H^f, S^c is presented 



in Fig. [4.5| obtained in the Borel scheme of NRQCD sum rules. 

In order to calculate the exclusive widths we suggest that the dependence of formfactors on the 
transfer momentum has the form of pole 

e'^iw) = eo-^-^, (4.70) 



pole 



with 



^poie = 6.3 GeV for b ^ c transitions, 
^poie = 1-85 GeV for c — * s transitions. 

The results of calculation for the exclusive decay widths of doubly heavy baryons in the framework 
of NRQCD sum rules are given in Table [4.11| , where the total widths have been supposed equal 
to the mean values presented in the previous section. The contribution of cascade decays into the 
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Table 4.11. The branching ratios (Br) for the exclusive decays of baryons with two heavy quarks. 



Mode 


Br (%) 


Mode 


Br (%) 


00 be 1' 


14.9 


6c cc ' 


4.9 


so ^E+Wi 

DC CC <• 


4.6 


-+ "0/,, 

be "b i 


4.4 


be b ' 


4.1 


cc c ' 


16.8 


cc c t 


7.5 


' — ^66 ^ ' — ^6c^ 


2.2 


00 OCr 


5.7 


' ^6c ^ ' 'cc ^ 


0.7 


' — 'be ^ ' — 'cc'^ 


0.7 


' — ^6c ^ ' — 'cc P 


1.9 


"be ^ "ceP 


1.7 


"6c ^ "6^ 


7.7 


^bc ~^ "6 


7.1 


"6c ^ "6P 


21.7 


"6c ~^ ^b P 


20.1 


cc c 


15.7 


"cc "c^ 


11.2 


' — 'cc ^ ' — 'c P 


46.8 


' — 'cc ^ ' — cP 


33.6 







baryons with the spin | is also taken into account in the table. In this estimate the results of 
on decays '^cc + were used. In the decays with S^^ and S*^ the factor caused by the 

anti-symmetrization of identical quarks was taken into account. For the decays of — >■ S+X the 
correction factor of 0.62 was introduced because of negative interference of Pauli. The author of |9^ 
claims that the values obtained in these sum rules agree with the results in the potential approach 
of 1 99 1 as well as with the above estimates of inclusive decay widths if we sum up the corresponding 
exclusive widths calculated in the sum rules. 



4.5. Discussion 

In this chapter we have calculated the lifetimes of baryons with two heavy quarks in the framework 
of consistent consideration of OPE in the inverse powers of heavy quark masses. In this expansion 
the leading constribution is determined by the spectator widths of heavy quarks in the inclusive 
decays, while the significant corrections appear under the taking into account the effects of Pauli 
interference and weak scattering, which contribute about 20-30% for the baryons Sec and Hbc- The 
measurement of lifetimes for the doubly heavy baryons allows us to make the comparative analysis 
of deacy mechanisms for the hadrons with heavy quarks, that is especially actual in the light of 
searching for the fine effects with the violation of combined -parity in the sector of heavy quarks, 
because the characteristics of quark inteactions enter the measured quantities with the factors 
composed by the hadronic matrix elements of quark currents. A reliable knowledge of properties for 
such matrix elements could be essentially enriched by the study of decays and lifetimes of baryons 
with two heavy quarks. Such the investigations allow us to numerically analyze possible effects 
with the violation of quark-hadron duality (these effects should be small, as we have stressed). An 
actual challenge is studying the dependence of heavy quark decay widths on the contents of hadron, 
which could be essential in the clarification of reasons causing the large deviation from unit for the 
ration of Af, and B lifetimes. The measurement of doubly heavy baryon lifetimes allows us also to 
investigate the characteristics of confinement for the heavy quarks in various systems. 

An open field of activity is a calculation of exclusive decay widths for the baryons with two heavy 
quarks. The presented estimates in the framework of NQRCD sum rules and potential models are 
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Figure 4.6. The decays of Ehc — > D^*^K^*^p in the process with the weak scattering of constituent b and c 
quarks. 

quite preliminary, of course, since the challenge here is a role of corrections in the inverse powers of 
heavy quark mass, which could be very significant for the hadrons containing the charmed quark. 

Moreover, the studies of exclusive hadronic decays for the doubly heavy baryons are of interst for 
the experimental practice. Some channels have to be distiguished among the hadronic decays. In 
contrast to the cascade decays of two heavy quarks, for example, S^c — > '^cc + X Ec + X — > E + X, 
which require a reconstruction of three "secondary" vertices, there are the processes with the weak 
scattering of constituents, in part, the heavy quarks, that will lead to the only heavy quark in the 
final state as it happens in Fig. for S^c, so that the only heavy hadron should be detected. The 
contribution of weak scattering into the total width is large enough (about 20%). Simple estimates 
of suppression factors show that the branching fraction of such decay is about several tenth per 
cent. 

We believe that the experimental investigations of doubly heavy baryons are quite actual prob- 
lem, first of all, at the hadron colliders, and the measurements of their decay characteristics could 
significantly enrich the knowledges on the mechanisms of heavy quark decays. 
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Conclusion 

In this review we have presented the basic physical characteristics of baryons containing two heavy 
quarks. The description of such hadrons is based on the hierarchy of scales in the interactions 
determining the time intervals and distances ascribed for the strong interactions in the subsystems 
composing these baryons. So, because of the nonrelativistic motion of heavy quarks with respect 
to each other the time interval for forming the system of two heavy quarks is greater than the 
time period for the hard production of heavy quarks or "dressing" the quarks by hard gluons. On 
the other hand, this time for the forming the heavy diquark is less than the characteristic time 
interval for the interactions providing the confinement of light quarks involving the low-frequency 
strong-interacting fields. Due to the mentioned hierarchy of strong interactions in the quark systems 
under consideration we have developed the methods of heavy quark effective theory for the baryons 
with two heavy quarks. In the effective theory we have isolated the leading approximation and 
constructed the systematic way for the calculation of corrections to it. 

Having started from the general approaches in the description of hadronic systems with the 
heavy quarks, we have developed the following particular methods for the consideration of baryons 
with two heavy quarks: 

- the potential of static quarks and the potential model for the doubly heavy baryons, the separation 

between the heavy quark motion inside the diquark and the light quark motion in the field of 
diquark; 

- the formulation of two-point sum rules in NRQCD for the quark currents corresponding to the 

baryons with two heavy quarks, and the calculations of both ground state masses for such the 
baryons and their coupling constants with the currents, taking into account the corrections 
to the local condensates of light and strange quarks; 

- the calculation of anomalous dimensions for the baryonic currents with two heavy quarks in 

NRQCD; 

- the scaling functions of diquark fragmentation; 

- the numerical calculations for the complete set of diagrams in the fourth order of QCD and the 

analysis of higher twists over the transverse momentum; 

- the generalization of OPE for the inclusive widths of baryons with two heavy quarks and a single 

light quark; 

- the formulation of three-point sum rules in NQRCD for the exclusive semileptonic decays and 

hadronic decays in the approximation of factorization for the transition current. 

The most bright physical effects in the baryons containing two heavy quarks are the following: 

- there are the families of doubly heavy baryons, i.e. the systems of quasi-stable excited levels for 

the baryons with identical heavy quarks due to a suppression of operators determining the 
transitions into the low-lying and ground states, since the quantum numbers of diquark have 
to change in such transitions, so that for some states the operators, which are not suppressed 
by the heavy quark mass or small size of diquark, have to be equal to zero because of the 
Pauli principle; 
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- the cascade processes of fragmentation at high transverse momenta, wherein we can get analytical 

results for the universal fragmentation functions in the perturbative QCD for the heavy quark 
into the heavy diquark and for the heavy diquark into the doubly heavy baryon, while the 
evolution of fragmentation function because of emission of hard gluons can be described in 
the framework of renornalization group in QCD; 

- the separation of regimes for the fragmentation and recombination in hadronic processes by taking 

into account the higher twists over the transverse momentum, which can be described in the 
framework of peturbative QCD by calculating the complete set of diagrams in the given order 
of coupling constant; 

- large, about 20-50%, contributions of nonspectator terms into the lifetimes of doubly heavy 

baryons depend on the valence quark contents of baryons and are given by the effects of 
Pauli interference and weak scattering, especially in the presense of charmed quark, that leads 
to a strong splitting of lifetimes: 

r[Si+] > r[Qtc] > r[S+], 

^ T[n-,] > t[E%], 

- cascade mechanisms of decays for the baryons with two heavy quarks, while there are peculiar 

modes due to the weak scattering, which have sizable branching ratios. 

Sure, direct measuring the masses of ground states and excited levels allows us to essentially 
move forward more deep understanding the dynamics of forming the bound states with the heavy 
quarks. The theory for the radiative transitions between the quasi-stable levels of doubly heavy 
baryons due to either the electromagnetic or strong forces should be developed, wherein the method 
of chiral lagrangian should be generalized for the soft emission of Goldstone's mesons, particularly, 
pions. 

A quite complete picture of production mechanisms has been constructed for the doubly heavy 
baryons, so that the search for these baryons at hadron colliders with high luminosities could be 
prospectively proposed. The distributions over the transverse momentum of baryons could give the 
most important information on the regimes of production and probably make a contribution into 
the understanding of reasons resposible for the discrepance between the theory and the experimental 
data on the yield of 6-hadrons in hadron collisions. 

The experimental data on the lifetimes of doubly heavy baryons could represent the most inter- 
esting information, since they are strictly connected to the whole system for describing the inclusive 
decays of heavy hadrons. The branching ratios of semileptonic decays could be important for de- 
termining the role of gluon corrections to the nonleptonic lagrangian for the weak charged currents 
of quarks. The lifetimes would allow us to enrich the knowledge on the masses of heavy quarks 
and the relative contributions of various decay mechanisms, since the total widths are extremely 
sensitive to the above physical characteristics. 

We could suppose that the description of exclusive hadronic decays would require extensive 
theoretical efforts for the baryons with two heavy quarks. 

Thus, the physics of baryons containing two heavy quarks is sufficiently rich and informative, so 
that it will justly occupy a honour place in the experimental investigations, in the light of optimistic 
occurancies and tendancies depicted by the modern theory, which deals with this field entering the 
time of mature progress. 
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Appendix I. The coefficients of spectral densities 



The spectral densities in ( p.l9|) have the coefficients 

r/i,i(cu) = 104uJ{231Mlg + 297MlgUJ + 121Md^gUJ^ + 15uj'^), 



'dig' 

The coefficients of spectral densities in (|2.2CI| ) have the form 



r/2,0 = A2u{Ml^^ + A%Mdiqu; + Uu''), 

^2,2 = -r-^— -^(105A<|, + 315^^ 27g_^ 2^77^3)^ (A.I.2) 

The spectral densities taking into account the coulomb corrections in ( p.24 ) have the coefficients 
r^i^o = (27Wd,g + cj) V, 

= ^Md^^^ + ^)2 + 24A<i.^ + 32.M^,,u;^ + 20A^,,,u;^ + 5a;^). (A.I.3) 



The coefficients of spectral densities taking into account the coulomb corrections in (|2.25| ) have the 
form 



T]^Q = {2Mdiq + UJ)UJ 

2 

in A / 

^diq 



J^diq + ^ 

"^So = T-TT^— ^(2A<^,, + 2A^,,,a; + a;2). (A.I.4) 
[Mdiq + ujy 



For the spectral densities with the gluon condensate in ( |2.3CI| ) we find 



„2 2(210A<3 + 70A<2 + 2lMd^qUJ^ + 3cu3) 
Via = 



Mdiq + ^ 
{Mdiq + 

Appendix II. The distribution over the transverse momentum 



_ 2{2lQMl^ + 7QMl^UJ + 2lMd^qUJ' + ^UJ'') 

^^1,2 — TT~A \ — Tv) • (^A.i.bj 



For the fragmentation of vector diquark into the baryon state in the scaling limit we get the following 
distribution over t = px/M with respect to the fragmentation axis: 

64^2 |i?(o)|2 1 



D{t) 



SItt 3(1 - r)5M3 t6 
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|t(30r3 - 30r^ - 61tV + 45r¥ + sSrH^ _ 

-UrH'^ + 3t^ - 9rt^ + IbrH^ - Or^t^) - 
(30r^ - mrh'^ - 54r^t2 + 27rH'^ + 9t^ + 18rt^ - Gr^t^ + 

(1 - r)t 

'r2(l + t2)- 



+18r^t^ + 3r V + 3t^ - 6rt^ + 9r2t^)arctg f - 

24(2r3t + rt3 + r¥) ln(^^^^-^)}. (A.II.l) 



Appendix III. The spectator effects in the baryon decays 

For the effects of Pauli interference and electroweak scattering in the baryons Sec we have the 
following expresions: 

47r iriu 

([(^^^ - ^^^4^)(Q7«(1 - 75)c.)(g,r (1 - l,)q,) + 

(^^^ - ^^^^)(Q7.75C.)(g,7"(l - 75)g,)][(C+ + C^f + (A.III.l) 
J(1-A;^)(5C^ + C2 -6C_C+)] + 

[(^^^ - ^^^)(Q7«(1 - 75)c,)fe7"(l - 75)g.) + 

(^^^ - ^^^)(Q7«75C,)(g.7°(l - 75)g.)]fc^(5C^ + - 6C_C+)), 

^(1 - k^){Cl - C!))(Q7a(l - 75)Q)fe7"(l - 75)g,) + (A.III.2) 
k^{Cl - C7!)(Q7a(l -75)9)fe7"(l -75)g.)], 

2 

where p+ = Pc + Pq, P- = Pc — Pq z± = Further, for p+ and p_ we use their threshold values 

rric rric 

These expressions were derived with account of low-energy renormalization of nonleptonic lagrangian 
for the weak interaction of nonrealtivistic heavy quarks. This lagrangian has the form 

^{Cl - Cl)kHcT,d)idT'^c) + ^{Cl - Cl)k'^{k-l - l){cT,t-c)j; - (A.III.3) 
+ C^f + i(l - kl){hCl + Cl- 6C+C_)](cr^c + h^,^,c){uV^u) - 

o o 
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o o 

^-k'^k-i - l){bCl + C!)(cVc + h^^^,t''c)r}, 

where = 7^(1 — 75), k = {as{fJ.) / as{mc)) and = U'y^t"'u + d'^^t"'d + sj^t^-s is the color current 
of hght quarks (i" = A"/2 are the color generators). 

Here we would like to make a note on the effective lagrangian terms containing the colored 
current of light quarks, we have neglected those terms, since they enter into the lagrangian with 

2 

the factor of k~9 — 1 numerically about 0.054. 
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